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RECENT APPLICATIONS OF CONVEX FUNCTIONS* 
J. W. GREEN, University of California at Los Angeles 


1. Introduction. A comprehensive account of the history, properties, and ap- 
plications of convex functions up to 1946 has been given by E. F. Beckenbach 
[1]. However, even in the few years since that article appeared, convex func- 
tions have been so consistently of value in analysis, geometry, and other 
branches of mathematics, notably mathematical economics, that a considerable 
wealth of material exists for a supplementary article. In the following I will re- 
port on some of this material. Coverage of the field will not be exhaustive; the 
choice of subject matter was determined by my own fields of interest; namely, 
convex functions per se, and their generalizations, geometry of convex figures, 
function theory, and (through various logistics and programming projects) lin- 
ear and non-linear programming. 


2. Convex functions and generalizations. We recall that a convex function f 
is one which satisfies the inequality 


fQx + (1 — Ay) S + (1 — MSY) 


for all \ such that 0S 31. The variables x and y may be regarded as vectors 
in any linear vector space, which in this article will be taken to be Euclidean. 
This inequality states that between any two points of the graph of f, the graph 
lies nowhere above the chord joining those two points. For elementary proper- 
ties of convex functions, see [1] and [2]. 

Now, as far as the study of properties of convex functions, regarded as func- 
tions of one or more real variables, is concerned, there has not been a great deal 
of activity of late. This is due perhaps in one dimension to a shortage of unsolved 
problems and in more dimensions to a shortage of solvable ones. Here we shall 
describe several generalizations of convex functions which seem interesting or 
useful. 

First, let us stay in one dimension and suppose that there is a family F of 
continuous functions F such that through any two points (x, yi) and (x2, ye) 
with x; x2 we may interpolate exactly one function F of F. A given function f 
is said to be convex relative to 7, or sub-F, provided it agrees with F at x, and 
%2, it lies nowhere above F for x between x; and x. It is clear that the functions 
F play the role which the linear functions play in the ordinary convexity. The 
sub-F functions were studied extensively by Beckenbach and Bing [3] and 
were reported on in [1]. However in the recent past, a number of papers [4, 5, 
6] have appeared concerning these functions, particularly in connection with 


* This article was originally prepared to be presented at the Christmas, 1952, meeting of the 
Association in St. Louis. However, an overly convex appendix on Christmas eve made its presenta- 
tion in person impossible. 
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their smoothness and convergence properties. For example if the F are sufficient- 
ly well-behaved, the sub-F functions have very nearly the same differentiability 
properties as the convex functions—unilateral derivatives exist everywhere, 
derivatives everywhere with countable exceptions, second derivatives almost 
everywhere. A convergent sequence of sub-F functions converges uniformly, 
and the sequence of derivatives converges boundedly. 

Now the sub-F functions, instead of being only generalizations for generali- 
zation’s sake, have proved to be applicable to a number of problems. In fact, 
a particular case of them, the sub-trigonometric functions for which 7 consists 
of the functions A cos x+B sin x, were introduced by Phragmén and Lindeléf 
and used by Pélya in connection with problems in function theory. Pélya [7] 
proved that a function p(@) of period 27 is the support function of a convex area 
in the plane if and only if it is sub-trigonometric. The present author rediscov- 
ered this fact [8] and made use of it in several geometrical problems (see §3 
below). In [8], proper acknowledgement is not made to Pélya, and the author 
wishes to make it here. In [6], by use of the differentiability and convergence 
properties mentioned above, a simple proof is given of the known but appar- 
ently not easily accessible fact that the area of the general convex curve is 
given by (p?— p’*)d@ without any smoothness assumptions. 

Another and novel generalization of convex functions has been given by 
Hyers and Ulam [9]. They call a function approximately convex, or €-convex 
provided it satisfies the inequality 


fQx + (1 — Ny) S + (1 — +, 


where e>0 and x, y are vectors in m-space. The chief theorem about an e-convex 
function is that it is approximately equal to a convex function. More precisely, 
there is a convex function g such that gSf<g+k,e, where k, depends only on 
the dimension n. The best value for k, has not been determined for all . For 
n=1, it is 1, and it is easily seen that in one dimension, the most general e- 
convex function is obtained by adding to an arbitrary convex function an 
arbitrary function whose values lie between 0 and «. In [9], the value kn= 
(n?+3n)/(2n+2) is given, which is the best value for n=1, 2, but not beyond. 
The present author [10] has obtained an improved value of k, asymptotic to 
log: ”; this value is the best possible through m =3 but beyond that the question 
is open. 

A third type of function akin to convex and studied of late is the sub-addi- 
tive function, satisfying the inequality f(x+y) Sf(x)+f(y). Rosenbaum [11] 
has studied these functions systematically. If f is sub-additive and homogene- 
ous, it is convex; if f is convex for x20 and f(0)20, it is sub-additive. These 
functions have many properties analogous to those of convex functions, but 
behave quite differently as far as continuity properties are concerned. For ex- 


ample, there exist measurable everywhere discontinuous sub-additive func- 
tions. 
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3. Geometry. Each year dozens of papers on the geometry of convex bodies 
are written, and these will usually involve convex functions, directly or in- 
directly. However here are a few problems where the analysis of functions is 
used to produce geometric results. 

In [8] the following problem was considered: Let K be a convex curve lying 
interior to a circle C. What can be said about K if it subtends at each point of C 
the same angle a? It was shown that if r—a is an irrational multiple of z or is 
of the form (m/n)mr with m even and prime to m, then K must be a circle con- 
centric to C; in other cases—for example a =7/2—there is a considerable variety 
of curves K.* The author used the sub-trigonometric characterization of the 
supporting function (6) of K to solve several external problems involving pe- 
rimeter, area, diameter, and width of K for those a where there is any prob- 
lem. Since, for example, the length is given by />* p(6)d#, it is not surprising that 
this was possible. A typical result is that the maximum length of K is attained 
by a figure consisting of elliptical arcs pieced together with m corners of angle 

The same characterization also proved useful in connection with another 
geometrical problem. Let K be a convex curve, and let the plane be subjected 
to an area-preserving affine transformation T so as to minimize the length of 
TK. In [12] the present author showed that K minimizes the length if and 
only if the Fourier coefficients a, and 5, of p(@) vanish. The problem was then 
posed of finding the curve of area 1 which when transformed so as to minimize 
the length gives the greatest length. The partial result that the extreme figure 
must be a polygon of five or fewer sides was obtained by using variations of 
~(0) consistent with its sub-trigonometric nature and with the fact that a,=), 
13] The fact that it must be a triangle was subsequently proved by Gustin 

13}. 

It is interesting that new proofs of geometric theorems have been obtained 
by methods originally devised in connection with mathematical economics. 
Karlin and Shapley [14] have employed the following theorem, due to Bohnen- 
blust and themselves [15]: Let ¢2(x) be a family of continuous convex functions 
in a compact convex set A in n-space with infzcs sup «¢a(x) >0. Thcn there ex- 
ists a linear combination of +1 of the ¢. with non-negative coefficients which 
is positive in A. Using this fact, Karlin and Shapley gave a new proof of Helly’s 
theorem and also of the theorem that if a spherical surface in m-space is covered 
by a compact set of hemispherical sub-surfaces, then it is covered by a certain 
n+1 of these hemispheres. 

4. Function theory. E. F. Beckenbach [16; references to related problems 
may be found therein] has recently obtained an interesting result relating the 
mean values of the modulus of a function analytic within a unit circle along 
radii and along concentric circles. Let 


* It is curious that Santal6é [20] has shown that if the very same problem is considered on 
the surface of a sphere instead of in the plane, one is led to the conclusion that K is in all cases a 
circle, no matter what a is. 
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1 1 Qn 
f ale) = | | de. 


Then ifv(¢) <1, (orevenif theaverageofvis £1) theny(r) | (0) | +2r(i- | (0) | ) 
for r<} and u(r) | | +(1- | f(0)|)/2(4 —r) for This result fol- 
lowed from the fact that u(r) is a convex function and 3” u(r)dr <1. These two 
facts, without further reference to analytic function theory, imply the stated 
inequalities on uw. It seems quite remarkable that inequalities such as these can 
be achieved so directly by consideration of simple geometric properties of con- 
vex functions. 


5. Programming. The typical problem of linear programming is that of 


minimizing a linear function g(x) = of m variables x1, x2, , Xn, sub- 
ject to a certain set of m linear constraints, f;(x)=b;— 0. ajx:20, j= 
1,2, ---+, m (x means the vector with components x, - - - , X,) and also the n 


special linear constraints x;20. A problem readily cast in this form is tha: of 
minimizing the cost g of a diet composed of n foodstuffs consumed in quanti- 
ties x1, ---, x, and containing m basic dietary ingredients. The inequalities 
{;20 correspond to stating minimum requirements of the various ingredients. 
The interpretation of x;=0 is clear. 

The above problem is that of determining the maximum of a linear (and 
thus convex) function in a convex polyhedron. In fact, it amounts to finding 
the distance from the origin to the supporting plane to this polyhedron normal to 
the vector (1, C2, - - - , €n). It may be stated in the form of a minimax or saddle 
value problem by the introduction of Lagrangian multipliers as follows [17]. 
First define $(x, 2) = g(x) + >? A,f;(x). A vector x° yields a solution of the maxi- 
mum problem if and only if there exists a vector 4° with nonnegative compo- 
nents such that 


$(x, 2°) S 2°) S a) 


for \;20, x;20; that is, if (x°, 4°) is a saddle point of the “surface” z=(x, 2). 

Related to the above problem is the dual problem; namely that of deter- 
mining the minimum of subject to the constraints A4;=0, 
One sees immediately that this problem leads to identically the same saddle 
problem as before, and indeed the two extremum problems are equivalent; if one 
has a solution, the other does, and the one minimum equals the other maximum. 

Fenchel [18, 19], dropping linearity entirely, has made a very broad generali- 
zation of these notions, based upon a certain interesting conjugate relation be- 
tween convex functions. This conjugate relation depends on the er rela- 
tion with respect to the quadric 


in (n+1) space. 
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Let f be a convex function of x=(x, - + + , x,) in a convex set D. The polar 
plane of a point (x, 2) is given by 


where =((&, - - , &n), are coordinates in another space in which 
it is convenient to imagine the polar plane. Then the function $(&) conjugate 
to f(x) is determined as the upper envelope of these polar planes as (x, z) varies 
over the graph z =/(x); that is, 


o(€) = sup (x-& — f(x)). 


The domain A of definition of @ is those points where the sup is finite. Now 
Fenchel shows (under certain continuity conditions on f) that A is a convex set 
and @ a convex function; furthermore the conjugate of ¢ is f. 

Now let f be convex in C and g concave in D, where C and D are convex sets. 
Let ¢ in T and y in A be the conjugate functions of f and g respectively (—y is 
conjugate to —g). Then the following two extremal problems are equivalent: 
(a) maximize g(x) —f(x) in C(\D; (b) minimize $(€) —(é) in TMA. Either both 
or neither of the problems have solutions, and the one maximum equals the 
other minimum. In case g(x) = )> ci; in the positive orthant D: x;=0, and 
f(x) =0 in the convex region C determined by the inequalities }> ajx;=b,, this 
pair of problems reduces to the dual linear programming problems mentioned 
earlier. 

It would appear that this duality should be capable of further exploitation. 
Incidentally, Fenchel [19] shows how it can be used to prove many of the classi- 
cal inequalities—Hdlder’s, etc. 
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INVERSE AND COMPLEMENTARY SEQUENCES 
OF NATURAL NUMBERS 


J. LAMBEK, McGill University, and L. MOSER, University of Alberta* 


1. In 1926, S. Beatty** proposed the following problem: If x is a positive 
irrational number, the sequences m(1+x) and n(1+x7) jointly contain one 
and only one term between each pair of consecutive positive integers. (Here 
and elsewhere in this paper, m and m range over the set of positive integers.) This 
was proved by A. Ostrowski and A. C. Aitken; the latter obtained a similar re- 
sult for rational x. 

We shall call two sets of positive integers complementary if they have no com- 
mon elements and together exhaust all positive integers. As usual, let [x] de- 
note the greatest integer in x. Then another way of stating the above result is to 
say that the sets of positive integers [xm]+m and [x-'n]+-n are complemen- 
tary. 

This problem has re-appeared in the literature at various timesf, in particu- 
lar its connection with Wythoff’s game§ has been noted. It seems however 
that no really wide generalization has been given, in spite of the fact that other 
such pairs of complementary sequences are easily constructed. Thus [e"]-++m 


* This paper was written jointly at the Summer Research Institute of the Canadian Mathe- 
matical Congress. 

** This MONTHLY, problem 3173, vol. 33, 1926, p. 159. 

tJ. V. Uspensky and M. A. Heaslet, Elementary Number Theory, New York, 1939, page 98. 
See also this MONTHLY, problems 4270, vol. 54, 1947, p. 549, and 4399, vol. 57, 1950, p. 343. 

§ W. W. Rouse Ball and H. S. M. Coxeter, Mathematical Recreations and Essays, London, 
1944, pp. 38-39. See also D. B. Sawyer, A Property of the Golden Section, Eureka, March, 1948. 
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and [log »]+n are complementary, and so are pa-+m and (n—1)+n, where 
bm denotes the m-th prime, and +(m) the number of primes not exceeding n. 

2. In what follows it will be convenient to use the term number to denote 
a non-negative integer or infinity. If f= { f(1), f(2), - +> } is a non-decreasing 
sequence of numbers, it must be of one of three types: 

(A) f() is finite for a finite number of only, 

(B) f(m) is finite for all m and constant from a certain m on, 

(C) f(m) is finite for all m but tends to infinity with n. 
We then define 


(1) ft(m) = number of m such that f(m) < n. 


It is easily verified that f is of type A, B, or C if and only if f* is of type B, 
A, or C respectively. 

Two sequences f and g of numbers will be called inverses provided, for all 
pairs of positive integers m, n, 


(2) S(m) < n or g(n) < m, but not both. 
Note that the inequality f(m) <n is false when f(m) is infinite. We shall prove 


THEOREM 1. A sequence of numbers has an inverse if and only if it is non-de- 
creasing. The inverse of f is unique and coincides with f* 


Proof. Assume f has an inverse g. If g’ were another inverse, (2) would im- 
ply that g(m) <m if and only if g’(m) <™m, for all positive integers m, n. If g(n) 
is finite, take m=g(n)+1. Then g’(m) <g(m) +1 so that g’(m) Sg(m) is also finite. 
Similarly, if g’(m) is finite, then g(m) Sg’(m). Hence, in any case, g(m) =g’(m). 
Thus g is unique. 

Next we wish to show that f(m+1)<4f(m). Without loss in generality we 
may assume that f(m-+1) is finite, say f(m+1)=n—1<n. Then, by (2), 
g(n)<m-+1, and a fortiori g(n)<m. Again by (2), f(m) <n, so that f(m)sn—1 
=f(m+1). Thus f is non-decreasing. 

Conversely let f be non-decreasing, so that f' is defined by (1). Consider 
the inequality f'(n) <m. According to (1), this means that the inequality f(x) <n 
has fewer than m solutions for x. Since f is non-decreasing, x =m is a solution 
if and only if all positive integers x <m are solutions, that is, if and only if the 
number of solutions is m or more. Thus f*() <m if and only if f(m) <n. In view 
of (2), ft is inverse to f. 

This completes the proof of theorem 1. We may remark in passing that 
(f')t=f, as follows immediately from theorem 1 and the symmetry of definition 
(2). We shall consider some examples of pairs of inverse sequences. 

Example 1. f(m) = pm, f'(n) =2(n—1). In general, if f(m) =m-th largest posi- 
tive integer with property P, or f(m) =infinity if fewer than m positive integers 
have property P, then f'(m) =number of positive integers less than n with prop- 
erty P. 
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Example 2. f(m) = [e], f'(n) = [log ]; f(m) =m?, f'(n) = [Vn—1]. In gen- 
eral, let ¢(x) 20 be any real-valued strictly increasing function defined for all 
real x=0, f(m)=[¢(m)], then f'(m)=number of positive integers less than 
(If is undefined, then f'(m) =0.) 

Example 3. 1(m?) and [+/>,] are inverses. In general, if f and g are non-de- 
creasing sequences of positive integers, then f'(g(m)) and g'(f(m)+1) are in- 
verses. 

3. If F isa set of positive integers, let F(m) denote the m-th largest member 
of F, as long as m does not exceed the number of elements of F, otherwise let 
F(m) denote infinity. The equation F(m) =f(m) +m establishes a one to one cor- 
respondence between all non-decreasing sequences f of numbers and all sets F 
of positive integers, as is easily verified. (It is understood that infinity plus or 
minus a finite number is infinity.) This correspondence will be exploited in the 
following 


THEOREM 2. Two non-decreasing sequences f and g of numbers are inverses if 
and only if the corresponding sets F and G of positive integers, defined by F(m) 
=f(m)-+m and G(n) = g(n)-+n, are complementary. 


Proof. Let f and g be inverses. By (2), the inequalities f(m) <n and g(n) <m 
are contradictories, in the sense that one and only one of them is true. Hence 
also the inequalities F(m)<m-+n and G(n) <m-+n are contradictories, and we 
conclude that F(m) is never equal to G(m). In particular, the set of all F(m) with 
m Sa and the set of all G(m) with n Sf(a) have no elements in common. Let f(a) 
be finite, then none of the elements in either of these sets exceeds a+/(a); for 
we infer from (2) that g(f(a)) <a since f(a) £f(a). Thus the two sets just con- 
sidered contain together the first a+f(a) positive integers. If f is not of type A, 
this argument holds for any a, hence the sets F and G are complementary. If f 
is of type A, then g is of type B, and the same argument can be used with f and 
g interchanged. 

Conversely, let F and G be complementary sets. Then f is a non-decreasing 
sequence of numbers, hence possesses an inverse f', by theorem 1. By the above, 
f'(n)+n defines the complementary set of F, hence coincides with g(m) +n. 
Thus g=f" is the inverse of f. 

This completes the proof of theorem 2. 

Example 4. We wish to calculate the m-th non-square. Let F(m) =m? be the 
m-th square. Then f{(m) =m?—™m and f'(n) =number of m such that m*—m <n, 
which is the same as the largest integer m satisfying this inequality. Since both 
sides of the inequality are integers, we may just as well compute the largest 
integer m satisfying m?—m+1/4<n, that is m—1/2<~+/n. Thus f'(n) = 
[1/2+/n]={/n}, where % denotes the closest integer to x. The m-th non- 
square is therefore n+{/n 

Example 5. In the same way we find that the m-th non-triangular number is 


n+{V/2n}. 
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4. If F is a set of positive integers, let CF denote the complementary set. 
We have shown in theorem 2 that CF(n) =n+f"(n), where f(m) = F(m) —m. The 
set F may be identified with the sequence { F(1), F(2),---> }, which is a non- 
decreasing sequence of numbers, and therefore possesses an inverse F'. The ques- 
tion arises: Can CF be expressed in terms of F* rather than f'? It turns out to 
be more convenient, if instead of F' we introduce F* defined as follows: 


(2’) F*(n) = Ft(n + 1) = number of elements of F not exceeding n. 
This may be supplemented by defining F*(«) = «©. We shall need the following 
LEMMA 1. n+ F*(CF(n)) =CF(n). 


Proof. When CF(n) is infinite, this is clear. Otherwise F*(CF(n)) is the num- 
ber of elements of F not exceeding CF(n). Now the number of all positive inte- 
gers not exceeding CF(n) is CF(n), and the number of elements of CF not ex- 
ceeding CF(n) is m. The result follows. 

We shall define recursively 


(3) F,(n) =n; F,(n) = n + F*(Fy_1(n)) for k > 0. 
Then we obtain the following 
THEOREM 3. CF(n) =limy..F:(). 


Proof. We first show that Fi.4:(m) = Fi(m). Clearly this is true for k=0. Sup- 
pose k>O and F;,(m) = Fy-1(m); then 


Fiyi(n) =n + F*(Fi(n)) = + F*(Fea(n)) = Fi(n), 


by (3), induction hypothesis, and the fact that F* is non-decreasing. 
We next show that F,(n) S$ CF(n). Without loss in generality we may assume 
that CF(n) is finite. Suppose F,_1(m) S CF(m), then 


F,(n) = n + F*(Fi-s(n)) S n + F*(CF(n)) = CF(n), 


by (3), induction hypothesis, monotoneness of F*, and lemma 1. 
Thus we may write G(n) =lim,..F:(n) SCF(n). Without loss in generality 
we may assume that G(n) is finite. Then for sufficiently large k, 


G(n) = Fya(n) = Fi(n) = n + F*(Fy-1(n)) = n + F*(G(n)), 


so that n=G(n) — F*(G(n)). By definition, this is the number of elements of CF 
not exceeding G(m). Thus n is the number of m such that CF(m) SG(n), whence 
CF(n) SG(n). Since we have already shown the converse inequality, we have 
equality, as was to be proved. 

Example 5. The m-th non-prime is the limit of the sequence: n, »+7(n), 

5. For given nm, the sequence F,(m), (k=0, 1, - - +), is a sequence of inte- 
gers, hence attains its limit CF(m) in a finite number of steps, provided this 
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limit is finite. In fact, it follows from theorem 3 that we need not go beyond 
k=CF(n)—n. It may happen that the number of steps need not depend on n. 
We shall see that in many cases two steps are sufficient. 


THEOREM 4. If F(m+1)—F(m) 2m, then CF(n) = F.(n). 


Proof. In view of lemma 1 and (3), it suffices to prove that F*(CF(n)) 
= F*(F\(n)). We know from the proof of theorem 3 that Fi(m) SCF(mn), hence 
we need only show that F*(CF(n)) S F*(F,(m)), which means that any element 
of F not exceeding CF(n) does not exceed F,(m). Assuming F(m) <CF(n), we 
shall prove that F(m) S F,(n). 

Since F and CF have no elements in common, we even have F(m) <CF(n), 
so that (CF)*(F(m)) Sn—1, hence by lemma 1, 


(4) F(m) = m + (CF)*(F(m)) Sm+n— 1. 


If m=1, we deduce immediately that F(1) <n<F,(n), by (3). There remains 
the case m>1. 

If m>1, by hypothesis and (4), F(m—1)+m—1S F(m) sm+n-—1, so that 
F(m—1) Sn, whence m—1S F*(n). Hence, by (4) again, 


F(m) Sn+m—1Sn+F*(n) = F,(n), 


as was to be proved. 


Example 6. The n-th positive integer which is not a perfect k-th power 
(k=2) is n+[(n+[n'/*])1*]. For k=2 this gives a new formula for the n-th 
non-square. 

Example 7. The m-th positive integer not of the form [e™] with m21 is 


n+ [log (n+1+ [log (n+1)])]. 


Postscript. It has come to the authors’ attention that Viggo Brun | Rechen- 
regel zur Bildung der n-ten Primzahl, Norsk. Mat. Tidsskr., vol. 13, 1931, pp. 
73-79] proved the following: If m=n—m(n+m,+ for k>0, then 
n+n+ --++m—p, as ko. This result is an immediate consequence of 
Theorem 3. For, letting CF(n)=p,, one obtains F*(n)=n—2(n) and n+m+ 

++ D. H. Lehmer [An inversive algorithm, Bull. Amer. 
Math. Soc., vol. 38, 1932, pp. 693-694] generalized Brun’s result from p, to 


any increasing sequence of positive integers, thus obtaining a result essentially 
equivalent to Theorem 3. 
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CARNOT AND THE CONCEPT OF DEVIATION 
C. B. BOYER, Brooklyn College 


In this MONTHLY for October, 1952 (p. 535), one reads that “The concept 
of aberrancy ... was originally introduced by Transon (using the name de- 
viation)” [1]. It is the intention here to point out that this notion is older still, 
for it appeared in 1803 in the Géométrie de position of Lazare Nicholas Margue- 
rite Carnot, published thirty-eight years before Transon wrote. The reputation 
of this book has been well-known, yet references to it leave the impression that 
it is devoted exclusively to synthetic geometry. It is indeed true that in it 
the author contributed significantly to so-called “pure” geometry; but the last 
of the five main sections of the work is concerned primarily with the use of ana- 
lytic methods. In fact, the point of view toward coordinates here presented is 
the broadest of any mathematician since Newton, another who traditionally is 
placed among the synthesists. Carnot may well have been unaware of the sug- 
gestion of polar, bipolar, and other coordinate systems in Newton's Methodus 
fluxionum, for this has been overlooked by other mathematicians and historians 
[2] just as has a similar contribution in the Géométrie de position. Under the 
heading, “De la détermination d'un point dans l’espace, et du changement de ses 
coordonnées,” Carnot proposed a wide variety of systems: polar, bipolar, bi- 
angular, and various triangular types. In such cases the equations of a curve de- 
pend on the elements used as a frame of reference; and so Carnot looked for in- 
trinsic coordinates. He wrote, 


For any curve there exists at each point a certain line which does not depend on any 
particular hypothesis or any basis of comparison taken in absolute space; this is the radius 
of curvature ...I believe nevertheless that it is possible to find various variables which 
shall have the condition demanded. 


Carnot considered using arc length as a second coordinate, along with the 
radius of curvature; but he rejected this because it depends upon the choice 
of a fixed point from which the distance is measured. Instead he introduced 
what now is known as the angle of deviation—or, strictly speaking, its comple- 
ment: 


I propose, for example, the angle which is formed, at the point describing the curve, by 
the tangent and the line which bisects the infinitesimally small secants drawn through the 
curve parallel to this tangent. 


In illustrating this idea, he replaced the tangent by the normal, thus bring- 
ing his angle into conformity with the modern aberrancy. If M is a point on a 
curve and K is the center of curvature for this point, Carnot drew perpendicular 
to MK a variable secant mm’ with midpoint n. His angle of deviation was then 
the angle between MK and the limiting position of the line Mn as m and m’ 
approach M along the curve. As an example of the use of the angle as a coordi- 
nate, he cites the parabola y?= px, for which one finds r=4p csc* z, where r is 
the radius of curvature and z the angle of deviation. This latter equation, in 
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curvature and aberrancy as coordinates, is one form of intrinsic equation of the 
parabola. Carnot also derived a general formula for cot z, pointing out that for 
a circle this expression vanishes; i.e., in modern terminology, the deviation of a 
circle at every point is zero [3]. 

In 1841 Transon gave a formula for tan z, introducing the phrase “axis of 
deviation” and pointing out that the deviation is a geometric representation of 
the third derivative. He also called attention to some remarkable properties of 
the deviation of conics, including the relation tan z=r’/r, where r and r’ are 
the radii of curvature of the conic and its evolute; and he introduced the center 
of deviation [4]. Walker recently has added striking relations between the de- 
viation of a conic and its differential equation [5]. The concept which Carnot 
introduced so unobtrusively thus has led to developments comparable to those 
connected with the idea of curvature; yet the adumbration in the Géométrie de 
position has been so completely overlooked that there seems to have been no 
reference to it by mathematicians or historians. Even Transon himself, ancien 
éléve of the Ecole Polytechnique which Carnot helped to establish, made no 
mention of it. 

The year 1953 marks the bicentennial of the birth of Lazare Carnot (as 
well as the sesquicentennial of his Géométrie de position), and hence brief atten- 
tion may appropriately be called to the extraordinary career of this versatile 
mathematician. He was born at Nolay (Céte-d’Or), on May 13, 1753, one of 
eighteen children, several of whom came to occupy positions of prominence. 
Lazare in 1771 entered the Ecole de Méziéres, where he studied under Monge. 
His fellow military students regarded him as un original because he preferred 
the library to the cafe. His early published works were both literary and sci- 
entific, for he wrote considerable poetry and took a particular interest in bal- 
loons, fortifications, and politics. When the Revolution opened, Carnot, then a 
captain, accepted enthusiastically the views of the National Assembly; but 
he took no active role until in 1791 he was elected to the Legislative Assembly. 
Subsequently designated to the National Convention, he was among those who 
voted for the execution of Louis XVI on grounds of “justice and policy alike.” 
Becoming more and more deeply involved in politics, he was chosen a member 
of the Committee of Public Safety in 1793. Charged with the organization and 
operation of the revolutionary army, his genius achieved such striking success 
against the invading forces that he was honored spontaneously as “The Or- 
ganizer of Victory.” Carnot did not hesitate to seize a musket and lead his 
troops against the enemy. Robespierre was bitterly jealous of his military abil- 
ity and St. Just accused him of moderantism, but he survived the Terror 
through his indispensability to the defence of France. However, he was pro- 
scribed in 1797 as the result of false charges of implication in a royalist plot. 
His name was stricken from the rolls of the Institut, and his chair as geometer 
was filled, mirabile dictu, by General Bonaparte. Monge, the great geometer, 
was so mesmerized by Napoleon that even he joined in this unanimously ap- 
proved outrage against mathematical honesty. 
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When Bonaparte scored his coup d’ Etat of 1799, Carnot returned to France 
to become minister of war and was reelected to the Institut. He was elected a 
member of the Tribunate, and, being an inflexible republican, he exerted all his 
influence against Napoleon’s imperial designs. In the midst of extraordinary 
factional turmoil, Carnot remained aloof from political parties, opposing 
tyranny of every description. A man of courage and conviction, in 1804 he de- 
livered a brilliant speech against the conferring of the Imperial Dignity upon 
Napoleon [6], but his was the lone negative vote. Yet, when in 1814 Carnot saw 
France in jeopardy, he offered Napoleon the services of his “weak sexagenarian 
arm.” As commander at Antwerp he held out against the enemy with notorious 
stubbornness, withdrawing only after the abdication of the emperor. Carnot did 
not conspire to recall Napoleon; but during the Hundred Days he again served 
his country, this time as Minister of the Interior, and he struggled in vain 
against the second abdication of Napoleon. He became a member of the provi- 
sional government, but he was exiled on the return of Louis XVIII and died at 
Magdeburg on Aug. 2, 1823 [7]. Carnot’s influence, however, was felt long after- 
ward: through his memory, in his descendants, and through his published 
works. His eldest son was the famous physicist, Sadi Carnot, author of Ré- 
flexions sur la puissance motrice du feu (1824); and a grandson, Marie Frangois 
Sadi Carnot, became the fourth president of the Republic. 

Notwithstanding a life of intense military and political activity, Carnot 
contributed regularly to mathematics. The work which enjoyed greatest popu- 
larity was his Réflexions sur la métaphysique du calcul infinitésimal of 1797. In 
this he sought to reconcile the various approaches to the calculus—the method 
of exhaustion, Leibnizian differentials, Newtonian fluxions, the limits of D’Alem- 
bert, the series of Lagrange, and others—concluding, unfortunately, that the 
true unifying basis was to be found in the principles of the compensation of er- 
rors. The Réflexions appeared in half a dozen French editions, as well as in many 
other languages, including English [8]; but as a permanent contribution to 
mathematics it falls far below the Géométrie de position of 1803. 

Monge and Carnot may be looked upon as the two founders of modern syn- 
thetic geometry, the former through his influence as a great teacher at the 
Ecole Polytechnique, the latter through his published works, especially the 
Géométrie de position. This contains numerous theorems of pure geometry, a 
striking example of which may be cited here: If M, N, P, Q are the areas of the 
four faces of a tetrahedron, and if m, n, p are the angles between the faces N 
and P, M and P, and M and N respectively, then Q?= M?+N?+P?—2MN 
cos p—2NP cos m—2PM cos n [9]. Yet the best known of Carnot’s contribu- 
tions to pure geometry comes not from the Géométrie de position, but from his 
Essai sur le théorie des transversales of 1806. This contains the celebrated 
“theorem of Carnot”—Given any algebraic curve of order m which cuts a tri- 
angle ABC, let A, be the product of the m distances, real or imaginary, from A 
to the points of intersection of the curve with the side AB, and let B, and C, 
be defined similarly for the sides BC and CA; and let A>, C2, and Bz be the simi- 
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lar products corresponding to the sides AC, CB, and BA respectively. Then 
A,B,C; = A2B2C>. If the curve is a straight line, one has a well-known theorem 
of antiquity; if the curve is a cubic, it follows from Carnot’s theorem that the 
three points of inflection lie on a straight line, a familiar result of the eighteenth 
century. 

Carnot was aware that one of the advantages of analytic geometry was its 
generality; and hence he had sought, in De la corrélation des figures de géométrie 
(1801), to establish for synthetic geometry a comparable universality. Here he 
proposed a unifying idea somewhat akin to the controversial “principle of con- 
tinuity” of Poncelet. This idea was further developed in the Géométrie de posi- 
tion. 

Among Carnot’s earliest works was the Essai sur les machines en général 
(1783), which contained the Carnot theorem on the loss of energy due to abrupt 
changes in velocity. This work, revised, appeared again in 1803 as Principes 
fondamentaux de I’ équilibre. 

There have been greater scientists and geometers than Carnot; there may 
have been greater political and military figures (although none more upright 
and courageous). But as a brilliant example of a man of thought and of action, 
it is difficult indeed to find a mathematician who outshines Carnot the elder— 
“The Great Carnot.” 
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MATHEMATICAL NOTES 
EpitTep By F. A. FIcKEN, University of Tennessee 
Material for this department should be sent to F. A. Ficken, University of Tennessee, 
Knoxville 16, Tenn. 
INEQUALITIES CONNECTED WITH DEFINITE HERMITIAN FORMS, II 
A. OPPENHEIM, University of Malaya, Singapore 


1. Let A =(ax%), B=(bx%) be two non-negative definite Hermitian matrices 
of order n. Then by considering the corresponding non-negative definite quad- 
ratic forms we see plainly that the matrix 


(1) C = (cix), Cik = Qik + Dix 
is also non-negative definite Hermitian. 


THEOREM 1. Suppose that the latent roots of A, B, C are the (non-negative) 
numbers arranged in increasing order: 


(2) @ Sa, 63825 --- SB, 3% 

Then 

(3) >>> & ++ + (Bi 1, 2,---, 0). 
If we take i= this inequality becomes (if a=det A, etc.). 

(4) clin > gi/n 4 


an inequality due (for real forms) to Minkowski [2] from which as I have shown 
[3] inequalities of Hadamard and Fischer are deducible. (See also Hardy, Little- 
wood, and Pélya, Inequalities.) 

Suppose that \ and yw are arbitrary non-negative numbers whose sum is 
unity. Then from (4) we obtain 


(5) det + 2 + 
and so by the inequality of the arithmetic mean and geometric mean, 
(6) det (Aain + (det (aix))M(det (bix))*. 


This inequality was given recently by Ky Fan [1]. Plainly (5) is stronger 
than (6). 
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2. The proof of Theorem 1 depends on an interesting algebraic inequality 
which is worth stating as a theorem. 


THEOREM 2. Suppose that the numbers e;, f; (t=1, - - + , m) are non-negative. 
Suppose that ef’, ff are the same sets arranged in non-decreasing order. Then if 


(7) II* = (e*# + + feta) + 
and 
(8) I]. = + + +54), 


(9) = = 


i=1 


Equality holds throughout (9) if one of the sets consists of equal numbers. Equal- 
ity holds on the right of (9) if one of the permutations which carries (e;) into (ef) 
carries (f;) into (ff). Equality holds on the left if a permutation carrying (e;) into 
(ef) carries (f;) into the reverse set to (ff). 


The second inequality in (9) is a special case of an inequality of Ruder- 
man [4]. 

We may suppose without loss of generality that ef =e; for each 7. This is 
merely a renumbering. 

For n= 1 the theorem is trivial. For n = 2 it follows from the obvious identity 


(x + y’)(x’ + y) — (x + y)(x’ + y’) = — x)(y’ — 9) 


when we take 


Consider the product 
II = + + fa) (en + fn)s (n = 3). 


If in any pair (e;+fi) (¢ii1+fi41) we have unequal fj, fis1, then by the case 
n=2 we get the minimum if f;<f,4: and the maximum if f;>f(4:. Thus after a 
finite number of steps we reduce [| to []« by interchanges which do not increase 
[[ and likewise we can change | to []* by interchanges which do not decrease | J. 
The maximum number of steps needed will be n(n —1). It is clear that Theorem 
2 follows at once. 


3. Suppose that A is positive definite. There is a non-singular transforma- 
tion which carries A and B into diagonal form with elements a; (in some order) 
and #; (in some order). This transformation carries C into the diagonal form 
;: (in some order) so that 


where the \,; are a permutation of the a; and the yu; are a permutation of the A. 
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By Theorem 2 


2 TL 2 I] (a; +8), 


j=l j=l 


since \¥, - - - , Af are a selection (in non-decreasing order of magnitude) from 
the non-decreasing set - , a, so that 

2 a;, and likewise = B;. 

i Now apply Hélder’s inequality to obtain 


Theorem 1 follows therefore when A is positive definite, B non-negative defi- 
nite. But if A and B are both non-negative and neither is positive definite, there 
is nothing to prove. 

; If we apply (3) to the matrix \A+yB where X and yp are positive numbers 
with sum unity and then use the inequality of the arithmetic mean and the 


geometric mean we obtain 


THEOREM 3. LetyiSy2S Sn be the latent roots of the matrix C=\A+pB 
AZ0, w20, A+u=1, where A and B are non-negative Hermitian matrices with 
latent roots 


@ Say, Bi S Be S Ba; 
then, fori=1,2,---,n, 


Theorem 3 is due to Ky Fan (Problem 4430, this MONTHLY, proposed in 
vol. 58, 1951, p. 194; solution in vol. 60, 1953, p. 50). The inequality is weaker 
than that of Theorem 1. 


4. I am indebted to the referee for giving an alternative proof of Theorem 1 
which employs another result of Ky Fan (Problem 4429, ibid., proposed in vol. 
58, 1951, p. 194; solution in vol. 60, 1953, p. 48). 

“Let A, B be two non-negative Hermitian matrices of order n. Let the eigen- 
values of A, B, A+B be 


B: Bi S82 S Ba, 
A+B: °°: 3 


respectively. Let x1, - - - , Xx, be m orthonormal eigenvectors of A+B such that 
(A + B)x; = 5%; (l1Sj Sn). 
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According to the theorem in Problem 4429 we have 
as II (Ax;, Bi S II (Bx;, x5), 
j=l j=l 
and therefore 
i i l/é 
j=l j=l 
By Hélder’s inequality, the right side is not greater than 
i 
[ ((A + B)x;, | = i) 
j=1 


Hence 


But the proof of Theorem 1 given in Section 2 does not use the theorem of 
Problem 4429. 
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All material for this department should be sent to G. B. Thomas, Department of Mathe- 
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REMARK ON A GRAPHICAL PROCEDURE 
W. T. Guy, Jr., The University of Texas 


In a recent note, Russell [1] extended the graphical method of Liénard 
[2,3] (used in plotting phase trajectories in non-linear problems) to some differ- 
ential equations of the type 


dy _ + 


1 
dx 
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Russell’s construction is based on the curves I';: [x = —¢(y) ] and Tz: [y= —y(x)], 
where $(y) and ¥(x) denote continuous functions. 

The purpose of this remark is to point out that his construction, with a slight 
modification, is also valid for certain differential equations of the type 


dy 8% 
dx h(x, y) 


(2) 


We assume g(x, y) =0 defines x as a continuous (single-valued) function of 
y [say, x= —$(y) ] and h(x, y) =0 defines y as a continuous (single-valued) func- 
tion of x [say, y= —y(x) ] over the required domains. Hence, g(x, y) = [x +4(y) ]- 
c(x, y), h(x, y) = [y+ (x) ]k(x, y), and equation (2) becomes 


dy +x c(x, y) 
dx Y(x)+yk(x,y) (x+y 


Here R(x, y) might be called a rotation factor. The construction of Russell based 
on his curves I’; and I; is repeated except now the line PD must be replaced by 
a line PE whose slope is equal to that of PD when multiplied by the value of 
R(x, y) at point P. 

The justification and limitations of the method are clear. If R(x, y) =1, the 
Russell construction obtains. If R(x, y) =a constant, the indicated modification 
is easily carried out. A protractor graduated with numbers indicating slopes 
will be of help in measuring the slope of PD and in «!rawing the line PE, es- 
pecially for the case R(x, y)#¥ a constant. This precedure has some merit over 
the usual directional field scheme, even in this latter case, if T, and T, can be 
easily determined and if R(x, y) is simpler to use than g(x, y)/h(x, ¥). 


(3) R(x, 9). 
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COMMENTS ON THE POLAR PLANIMETER 
L. I. Lowe, Northrop Aircraft 


Figure 1 shows a popular design of planimeter used in engineering offices for 
measuring areas of various shapes. No doubt the question has occurred to many 
operators: “Is the total angle generated by the indicating roller, as the point 
traces a closed curve, in direct proportion to the area encompassed by the curve?” 

The following verification of the planimeter would seem worthwhile, first, 
because it is more straightforward and complete than the brief statement in 
Calculus books that “the polar planimeter is an instrument which makes use of 
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Guldin’s Formula for measuring areas.” Secondly, it provides a valuable ex- 
ample of a curvilinear coordinate system and the teacher of Advanced Calculus 
is frequently hard put for examples when covering that topic. 

For purposes of analysis the planimeter will be thought of as a coordinate 
system. The arms have arbitrary lengths, a and b, which are constant, but can 
be assumed long enough for coverage of any arbitrary area. The variables, a 
and B, may be related to the cartesian variables, x and y, by means of the rela- 
tions: 


(1) x = acosa+ bcosB 
(2) y = asina+ db sin B. 
The differential element of area is: 
(3) dA =|J| 
Ox Oy 
da da 
in which J= 
Ox dy 
ap ap 


Computing the partial derivatives from the basic equations, we find that 
the Jacobian becomes: 


(4) J = — absin (a — 8). 
The element of area is thus found to be: 
(5) dA = — absin (a — B)dadB. 


Now let the planimeter trace an element of area, AA, and an accounting 
will be made of the motion of the indicating wheel. We see from the drawing that 
while the point is moving from 1 to 2, a is being increased by Aa but there is no 
relative motion between the two arms. From 2 to 3, 6 is being increased by A8 
while a is held fixed. The remaining two sides of the curvilinear quadrilateral 
are traced by motions which are the exact reverse of those which moved the 
point from 1 to 2 and 2 to 3, respectively. The wheel movement resulting from 
the point’s moving from 2 to 3 is equal and in the opposite direction to that 
which occurs as the point moves from 4 to 1. Hence, these movements cancel 
each other and will be neglected. 

As the point moves from 1 to 2, the total wheel rotation will be some multi- 
ple of 


(6) a-Aa cos (a — 8); 
from 3 to 4 it will be the same multiple of 


(7) a-Aa cos (a — {8 + As}). 
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Subtracting (7) from (6) and rearranging will yield: 


(8) a-Aa[cos (a — 8)(1 — cos AB) — sin AB sin (a — 8) ]. 
Tracing point 
«-B b x 
Y | 
|3 
Indicating wheel 2 
a + 
ls 
§ 
; 


x 


Fic. 1. The Polar Planimeter as a Coordinate System. 


Tracing point 


Fic. 2. Element of Area. 


If we take Aa =da and A8=df and neglect all powers of A higher than the first 
power, we may replace (8) by 

— asin (a — 8)(dadg) 
which is a multiple of the right hand side of equation (5). 
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ELEMENTARY PROBLEMS AND SOLUTIONS 
Epitep By Howarp Eves, University of Maine 


Send all communications concerning Elementary Problems and Solutions to Howard 
Eves, Mathematics Department, University of Maine, Orono, Maine. This department 
welcomes problems believed to be new and demanding no tools beyond those ordinarily fur- 
nished in the first two years of college mathematics. To facilitate their consideration, solutions 
should be submitted on separate, signed sheets, within three months after publication of 
problems. 


4 


PROBLEMS FOR SOLUTION 
E 1126. Proposed by J. V. Pennington, Houston, Texas 


Smith: “Down in Todd County, which is a 23-mile square, I have a ranch— 
oblong in shape—measuring a whole number of miles each way.” 

Jones: “Hold on a minute. I happen to know the area of your ranch; let 
me see if I can figure its length.” He figures furiously. “I need more information. 
Is the width more than half the length?” 

Smith answers the question. 

Jones: “I now know the length.” 

Brown: “I too know the area and, although I did not hear your answer to 
Jones’ question, I too can tell you the length.” 

Green: “I did not know the area, but I can now tell you the length.” 

How many miles long is the ranch? 


E 1127. Proposed by Edmund DeWan, Harpur College 


Enumerate the total number of squares appearing on an m Xn checker 
board, m=n. Solve the analogous problem for rectangular parallelopipeds. 


E 1128. Proposed by C. S. Ogilvy, Hamilton College 


Show that the point P at eye level from which a picture high on a vertical 
wall subtends the maximum angle can be found as follows. With eye-level line 
as directrix and bottom of picture as focus, draw a parabola; from the intersec- 
tion of this parabola with a line through the center of the picture parallel to the 
directrix, a line perpendicular to the directrix meets it at P. 


E 1129. Proposed by M. S. Klamkin, Polytechnic Institute of Brooklyn 


Find an integral arithmetic progression with an arbitrarily large number of 
terms such that no term is a perfect rth power for r=2, 3, ---, or N. Is this 
still possible if N= 0? 


E 1130. Proposed by Victor Thébault, Tennie, Sarthe, France 


Let the perpendicular bisector of the median BB’ of triangle ABC and the 
tangent at B to the circumcircle of triangle ABC cut the line AC in points M 
and N respectively. Show that triangle A BC is isosceles with vertex at A if and 
only if AM/AN=3/4. 
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SOLUTIONS 


E 1051 [1953, 551]. Correction. The condition given in the solution is suf- 
ficient but not necessary. 


E 1088 |1954, 347]. Correction. In line 3 of the solution, the statement 
“P,=1” should read “P,=P.” 
Differential Equation of All Conics of Given Eccentricity 


E 1070 [1953, 332; 1954, 50]. Proposed by A. W. Walker, University of To- 
ronto 


Find the lowest order differential equation satisfied by all conics with a given 
eccentricity e. 


Solution by Alexander Oppenheim, University of Malaya. The required 
differential equation is 
(1) (814s) = k[(i+ (2ute - 4ymu 4u 


where 


(2) k=(e —1)"/(2-e), 


and suffixes denote differentiations with respect to x. 
Note that 


(3) 2k = tan ¢, e = sec 9/2, 


where ¢ is the angle between the asymptotes of the conic. 
To obtain (1) it is best to use the general equation of the conic 


(4) ax® + 2hxy + by? + 2gx + 2fy+c¢=0 

and the relation 

(5) (a + b) tan @ = 2(h? — ab)!/?. 

From (4) 

(6) by = — hx — f + X"?, X = — Cx? — 76x —- A, 


where, as usual, C=ab—h?, G=hf—bg, - - - , are the cofactors of c, g,+ ++, in 
the determinant A of (4). Differentiate (6). Then 


(7) by, = — h — (Cx +G)X-"?, 
(8) by: = = X, 
by (2). Differentiate (8) twice and obtain 

(9) A?/84, = — 2(Cx+Q), 


(10) 
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whence, since AC—G?=DA, 


(11) (2uus — = 4b. 
We now employ (5) in the form 

since C=ab—h?. From (7), (8), (9) 


1/2 


(13) h = — — + (1/2)0 
so that, by (10), 


—2/3 


(14) 160 "(C+ = — 91 — — 


after a little reduction and by (11). From (10), (12) it follows that 


(15) = — + yx) — — J, 


since we may ignore }=0. Equation (15), of order four, is the differential equa- 
tion of all conics of given eccentricity e where & and e are related by the second 
equation in (2). 

Asa check take e= 1 org =0. Then (15) gives u2= 0, the well known differential 
equation of all parabolas. If we put ¢=7/2, we see that the differential equation 
of all rectangular hyperbolas is 


(16) (2uts — + yi) — — 4u = 02 


The derivative of the left hand side of (16) is 2uus(1+~7), so that it is easy to 
verify that elimination of k from (15) leads to u;=0, the differential equation 
for all conics. It is worth adding that (10) gives the differential equation of all 
conics with given area or with given product for semi-axes, since, if a, 8 are the 
semi-axes, = A?, 


The Counterfeiters of Lower Slobbovia 
E 1096 [1954, 46]. Proposed by L. R. Ford, Jr., University of Illinois 


As is well known, Lower Slobbovia is too poor a country to afford its own 
mint. There are N coiners engaged in making Rasbuckniks, the local currency, 
to government specifications. However it is suspected that some of them may be 
counterfeiting by introducing some base metal into the alloy. Any pair of 
counterfeits will weigh the same, although slightly different from the weight 
of a good coin. Each coiner produces either all good coins or all counterfeits. 
With one guaranteed good coin, a set of infinitely refinable weights, a beam 
balance, and as many coins from each coiner as may be needed, determine in 
three weighings whetuer any of the coiners is dishonest, and which ones. 
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Solution by Julian Braun, China Lake, Calif. First, determine the weight of 
the good coin, W,. Second, put one coin from each coiner on one side of the bal- 
ance and determine their total weight, T. The discrepancy is D=T—NW,,. If 
D=0, all the coiners are honest. Thirdly, if D#0, put 2*" coins from the ith 
coiner, 1=1, 2, --+, N, on one side of the balance and determine the total 
weight of these coins, 7’. The discrepancy in this case is D’ = T’(2" —1) W,. 

Find the integer S such that D’/D=5S/8(S), where 8(S) is the number of 
ones in the binary representation of S. Let S= >-™,B,2*, where B;=0 or 1. 
Then B(S) is the number of dishonest coiners and the ith coiner is honest or 
dishonest according as B;=0 or 1. 

Also solved by G. B. Charlesworth, Octave Levenspiel, P. B. Maggs and 
D. B. Mumford (jointly), R. H. Marquis, D. C. B. Marsh, L. V. Mead, C. S. 
Ogilvy, C. R. Perisho, H. W. Resnikoff, L. A. Ringenberg, Thomas Saaty, A. J. 
Tingley, and the proposer. These solutions showed interesting variations from 
one another. 


A 3:4:5 Right Triangle 
E 1097 [1954, 46]. Proposed by Leon Bankoff, Los Angeles, Calif. 
Arcs AB and CD are quadrants of circles tangent externally at their mid- 
points, E, and such that AC and BD when extended meet perpendicularly in F. 
A circle is inscribed in the mixtilinear triangle EDB, touching ED in M. G is 


the projection of M upon EF. Show that triangle MGF is a 3:4:5 right tri- 
angle. 


Solution by W. J. Cherry, North Riverside, Illinois. Let K be the center of 
arc AB. Taking BF and BK as x and y axes, denote the incenter of the mixti- 
linear triangle by P: (x, y). Let N and H be the projections of Pon BF and BK. 
Then, by applying the Pythagorean Theorem to triangles KHP and NPF, we 
find 


x?=4rny and (x — 7)? = r2(2y + 12), 


where r; and rz denote the radii of arcs AB and CD respectively. Taking r:=1, 
arbitrarily, and rz=+/2—1, and solving the equations of the two parabolas 
simultaneously, we get 


x= 2A2/2—1)/7, (9 — 4/2)/49, 
whence 
< PFN = arc tan y/(1 — x) = arc tan 1/7, 


and 


X< MFG = 7/4 — arc tan 1/7 = arc tan 3/4, 
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and the conclusion is established. 

Also solved by Hiiseyin Demir, E. I. Gale, A. R. Hyde, Beckham Martin, 
L. A. Ringenberg, David Rothman, C. W. Trigg, Roscoe Woods, and the 
proposer. 


A Set of Finite Geometries 


E 1098 [1954, 47]. Proposed by L. A. Ringenberg, Eastern Illinois State Col- 
lege 
Construct all geometries satisfying the following axioms: 
Axiom I. Space S is a set of m points, m a positive integer. 
Axiom II. A line is a non-null subset of S. 
Axiom III. Any two distinct lines have exactly one common point. 
Axiom IV. Every point lies in exactly two distinct lines. 


Solution by the Proposer. Since there is at least one point, there are at least 
two distinct lines (I, IV). Two such distinct lines, say a and 6, have a point A 
in common; to be distinct, at least one of them, say a, contains a point other 
than A, say B. Then B must lie on some third line c (III, IV) distinct from a 
and 6. It follows that the number m of lines in the geometry is an integer ex- 
ceeding 2. 

Two distinct pairs of distinct lines cannot have the same common point 
(IV). Hence »=C(m, 2)=m(m—1)/2. Thus a geometry with m lines has 
n=m(m—1)/2 points, each line consisting of m—1 points and each point lying 
on two lines. 

Following is a diagram showing how to construct the geometry with m=5, 
n= 10, the extension to an arbitrary m 23 is obvious. 


linea: A BC D 


b: A EF G 
CRF 
ec: DGI J 


Also solved by W. E. Briggs, Myles McConnon, D. C. B. Marsh, G. H. 
Meisters, J. D. Miller, D. B. Mumford, and D. R. Sudborough. 


Editorial Note. The postulates form an interesting but simple set for class- 
room purposes. The student should examine the set for consistence, independ- 
ence, categoricalness, and fertility. 

A Gambling Table 
E 1099 [1954, 47]. Proposed by W. R. Van Voorhis, Fenn College 


On a certain gambling table there are NV squares marked “2 to 1,” “3 to 1,” 
..., “N+1 to 1.” A sum of money is placed on each square and one of the 
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squares is selected as a winner, paying the player at the odds marked on that 
square. The player loses the amounts placed on the other squares. What is the 
maximum WN for which it is possible for a player to place money so that he can 
never suffer a net loss? 


Solution by C. F. Pinzka, Princeton, N. J. Let S, be the sum placed on the 
kth square and T the total amount bet. In order that there be no net loss, it is 
necessary and sufficient that (k+2)S,2T7 for k=1, 2,---, N. Writing this 
as 1/(k+2)<S,/T, and summing both sides from k=1 to k=WN, gives 

N 1/(k+2) <1. It is easily verified that 4 is the maximum N. 

Also solved by W. E. Briggs, S. H. Eisman, A. L. Epstein, J. M. Howell, 
P. G. Kirmser, M. S. Klamkin, Bart Park, F. D. Parker, L. A. Ringenberg, 
D. C. Russell, N. Shklov, Hans Ury, and R. Z. Vause. 

Ringenberg pointed out that not only can the player insure a net gain for 
N=1, 2, 3, 4, but that he can so bet that the net gain is independent of the 
square selected. If N=2, S,=4, S,:=3, the gain is 5; if N=3, S:=20, S.=15, 
S;=12, the gain is 13; if N=4, S;=20, S,=15, S;=12, S,=10, the gain is 3. 
Russell generalized the problem so that the stated odds are p; to1, pptoi,---, 
to 1. 


An Inequality 

E 1100 [1954, 47]. Proposed by L. E. Ward, Sr., Naval Ordnance Test Sta- 
tion, China Lake, Calif., and L. E. Ward, Jr., University of Nevada 

For x20 and #22, prove that (1+<x')!/*—(1+x')"/'<x. 

Solution by J. V. Whittaker, U.C.L.A. Setting 1+x‘'=y', raising each side 
to the tth power, and dividing through by y‘', we obtain the equivalent in- 
equality 

(1 — 1/y*)' 1 — 1/y', 
for y21, #22. But now 
(1 — 1/y*)* S (1 — 1/y')* S 1 — 


Equality is attained only when y=1. Hence in the given inequality, equality 
occurs only when x=0. 

Also solved by Tom Black, W. E. Briggs, A. E. Danese, H. M. Feldman, 
Norman Greenspan, A. R. Hyde, M. S. Klamkin, Viktors Linis, A. E. Living- 
ston, Lee Lorch, L. L .Pennisi, B. E. Rhoades, L. A. Ringenberg, R. Z. Vause, 
Chih-yi Wang, L. E. Ward, Jr., Albert Wilansky, R. E. Wild, and J. R. Winkel- 
man. 

Wilansky established the improved inequality 


where L is the original left hand side. 
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ADVANCED PROBLEMS AND SOLUTIONS 
EpitTep By E. P. Starke, Rutgers University 


Send all communications concerning Advanced Problems and Solutions to E. P. Starke, 
Rutgers University, New Brunswick, New Jersey. All manuscripts should be typewritten 
with double spacing and margins at least one inch wide. Problems containing results believed to 
be new or extensions of old results are especially sought. Proposers of problems should also en- 
close any solutions or information that will assist the editor. In general, problems in well 
known textbooks or results in readily accessible sources should not be proposed for this de- 
partment. 


PROBLEMS FOR SOLUTION 
4598. Proposed by I. J. Schoenberg, University of Pennsylvania 
Let the Laurent series 


fay as, <r, 


converge in an open ring containing the circle | z| =1. Show that a sequence 
{b,} exists satisfying the conditions 


(i) |b, | < K, for some K, for all = 0, +1, +2,---, 
(ii) DY anbm—n = 0 for all m, 
(iii) b, # 0 for some n, 


if and only if f(z) vanishes somewhere on the circle | 2| =1, 
4599. Proposed by Leonard Carlitz, Duke University 


Let p be a prime >2, a40 (mod p). Show that the number of solutions of 
the congruence 


(x + y + = 2axyz (mod 


is p?+1. 
4600. Proposed by Leon Bankoff, Los Angeles, California 
Vertices A—-C and B—D of square ABCD are joined by quadrants of circles 3 


(B) and (C). A semi-circle (0,) is described internally on the diameter BC and 
a circle (O2) is drawn tangent to the three arcs. Another circle (03) is drawn 
tangent to circle (O02) and to arcs AC and BC, and a right triangle is formed by 
joining Os to O, and dropping a perpendicular from O; upon BC. Successively 
tangent circles are drawn in the same manner (with (O,) tangent to (On_1) 
and to arcs AC and BC) and right triangles are formed (with 0,0, for hy- 
pothenuse). 
Show that the infinitude of triangles so constructed are Pythagorean. 
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4601. Proposed by Harry Goheen, Iowa State College, Ames, Iowa 


What is the necessary and sufficient condition that as io the limit of n, 
exists, if 2; is defined by the recurrence relation 


= bn; + 
no being given? 


4602. Proposed by C. M. Ablow and D. L. Johnson, Seattle, Washington 
Show that 


> A; cos (Bit + C3) 
tml 


changes sign as ¢t varies from zero to positive infinity. The A;, B;, C; are real 
constants; 


SOLUTIONS 


Squarefull Integers 


4459 [1951, 636]. Proposed by D. J. Newman, Republic Aviation Corporation, 
Farmingdale, N. Y. 


Find an asymptotic expression for the number of integers, not exceeding x, 
each of which has the property that each of its prime divisors divides it to the 
second power at least. 

Editorial Note. As has been pointed out by P. T. Bateman and others, the 
evaluation of the integral (1) of an earlier solution [1953, 719] is more than a 
direct application of Cauchy’s theorem on residues and its complete justification 
requires tools far from simple. 

Actually a rather good result can be obtained in this problem by a relatively 
elementary argument, as follows: 

2. Solution by P. T. Bateman, University of Illinois. If c, is 1 or 0 according 
as n does or does not have the property in question and if A(x) is the number of 
integers not exceeding x having this property, then 
(1) A(x) = » Cn 


naz 
and 


= 


as in the former solution. Let 5, be the number of ways in which the positive 
integer m can be expressed in the form m*n!, that is, as a product of a square and 
a cube. Thus 


i 
- 
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i=1 m=1 


n=1 


for R(s)>4. We consider the summatory function 
B(x) = 


mn?S 


If m and n are any positive integers such that m?’n*<x, then either m<x'/ 
or n Sx", possibly both. Thus 


B(x) = > 1+ > 


+ O(1)} + + O(1)} + O(1)}? 


mS 
= 4. 1/2 — 42/5 4 O(x1/5), 
ms ng 


Now we use the asymptotic formula 


n-* = — 1) + + O(y™), 


where r is some fixed positive number other than 1, and y21. (This formula is an 
easy consequence of the Euler-Maclaurin sum-formula of the first order; cf. 
Titchmarsh, The Theory of the Riemann Zeta-Function, Oxford, 1951, p. 14.) 
Using this in (2) we obtain 


B(x) = 4 ¢(2/3) + O(a-2/5)} 
(3) + { ¢(3/2) — 2x71/10 O(x-8/10) } — x2/5 4+ O(x!/8) 
= + ¢(2/3)x1/8 + O(x1/5), 
Now if R(s)>1, we have 
1/¢(s) = u(m)m-*, 
m=1 
where yu denotes the Mébius function. Accordingly 


(for R(s)>4) and thus 


Cn = 
jon 
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Hence by (3) 
A(x) = p(m)b; = B(x/m*) 
= nlm) + £(2/3) + } 
ms 
4) = { u(m)m-* + o( Zn m- 


+ (2/3)21 + 0( 


m> zi/6 


m3 


Now if r is a fixed positive number greater than 1, and y is positive, >. »>,n" 
=O(y'-"). Thus (4) gives 


A(x) = x1/2¢(3/2)/f(3) + 4+ + + 
+ + O(a"), 


This seems to be the best result readily obtainable by elementary methods. By 
more delicate methods, however, it is possible to sharpen the error term in 
(5) somewhat, e.g., O(x!/* log? x). 


(5) 


Inequality Involving Coefficients of Star-Like Functions 
4535 [1953, 268]. Proposed by M. S. Robertson, Rutgers University 
Let the function f(z)= a:*0, be regular for || and let it 


map |z| $1 onto a simply-connected domain D, star-like when viewed from 
the origin. Show that 


2 
2( SD (m+ 1)| + + +++ + |?, 


n=1 n=1 
with equality only for f(z) =aiz. 


Solution by W. C. Royster, Alabama Polytechnic Institute. Let f(z) = Re® and 
denote by C, the map of | z| =1 effected by w=f(z), and by C, the map of 
|z| =1 effected by w= [f(z)]*. Let A; be the area enclosed by C, and Az the 
area enclosed by C2, the area on separate sheets counted separately. Then 


1 1 
A,=—]| Rdg=—] R?—@, 
2 


* @ 
A -f Rig = | Rt 
2 do 4 


~ 
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where z=re, Since f(z) is starlike with respect to the origin, 0¢/092=0 and 
Jc" (0/00)d6 = 27. By the Schwartz integral inequality one has 


(j a0) f te ap, 
0 06 0 06 o 00 0 00 
Thus A,?S(2/2)As2, equality occurring only when R is constant. From this the 
required inequality follows since 


A,= (m+ 1)| aan +--+ + 
n=l 
with equality holding only for f(z) =a,z. 

It is interesting to note that if f(z) is p-valently starlike with respect to the 
origin then the right hand member of the given inequality is multiplied by p. 
This happens since />"(0¢/00)d0 = 2p for functions of this class. In this case 
equality holds only for f(z) =a,2?. 

Also solved by the Proposer. 


Uniform Convergence of Trigonometric Sums 


4536 [1953, 268]. Proposed by O. E. Stanaitis, St. Olaf College, Northfield, 
Minn. 


Prove that 


sin 76 sin n6 cos #76 sin 
n=1 n n=1 n 


are uniformly convergent in any interval, and that 


sin 276 cos 


n=l n 
is divergent. 


Solution by E. M. Wright, University of Aberdeen, Scotland. The first series is 
uniformly convergent. We have 


sin sin nO 
n=1 n n=1 n 


cos N(N + 1)@ cos n(n + 1)6 
N n(n + 1) 


cos n(n — 1)@ — cos n(n + 1)0 


which, as VN, becomes 


cos n(n + 1)6 
n=l n(n + 1) 
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uniformly in @, and the last series is uniformly convergent. Similarly for the 
second given series, since 

2 cos m6 sin n@ = sin n(n — 1)@ — sin n(n + 1)8. 


; The third series converges for some values of @ and diverges for others. Let 
us put = p7/gq, where , g are relatively prime positive integers and 


N 2 
S(N) = >> 
q q 


Since sin (m*pmr/g) cos (npw/q) has period g, we have S(rg) =rS(q) and 
N N N 4 
sin) = +s w - o[—]} = + sen. 
q q q 


Thus f(N), as just defined, satisfies |f(N)| <2g. Hence, if 
N sin 70 cos n@ N. S(n) — S(n — 1) 


n=l 


n 


As N—» this last series is absolutely convergent by the Dirichlet test. There- 

fore the given series with 6=pm/g converges or diverges according as S(g) =0 

or S(q) 

By combining the terms in S(qg) which correspond to n=r and n=q-—r, it is 
easy to see that S(g) =0 for all even g (p is then odd). For g=3, 5; S(g) #0, evi- 
dently. It is somewhat more difficult to see that S(g) #0 whenever q is odd; so 
that, in every interval, the given series converges for infinitely many values of 
6 and diverges for infinitely many others. 

The problem is given in Bromwich, Infinite Series, 1926, p. 68, ex. 16, except 
for a misprint: the third series is printed whereas the second is intended. 

Also solved by Joshua Barlaz, S. Parameswaran, and the Proposer. 


Inversion of Order of Summation 
4537 [1953, 268]. Proposed by Albert Wilansky, Lehigh University 


Given that converges (as an iterated sum) whenever isa 


sequence such that <0. Show that for 
each such {x,}. 


Solution by R. P. Agnew, Cornell University. The required result and, more- 
over, necessary and sufficient conditions that a matrix a, have the property 
: in question, follow immediately from theories of linear transformations orig- 
inated by Silverman and Toeplitz in 1913 and greatly extended by Hans Hahn 
(Uber Folgen linearer Operationen, Monatshefte fiir Math. und Physik, vol. 32 

(1922), pp. 3-88) and others. 
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If we set 
p @ P 

(1) Sp= Ape = Do 

n=1 kel n=1 
then 
(2) Sp A pk 


and the hypothesis of the problem becomes the hypothesis that S;, S2,---, 
and lim S, exist whenever >>| x;| < ©. Necessary and sufficient conditions for 
this are existence of constants A;, As, - - , such that 


(3) lim A px = Ax 


and existence of a constant M such that 
(4) |Ap.| $M k= 
Moreover, when (3) and (4) hold and >| xe| <0, 


(5) lim Sp = >> Axce. 
| k=l 


This is the required result. In terms of the given matrix a,,, the necessary and 
sufficient conditions (3) and (4) take the forms 


Pp 
Onk 


n=1 


(6) Gnk = Ax, M. 
n=1 


The latter condition implies, but is not implied by, the condition 
(7) | on. | 
An example, in which a,, is the amusing matrix 
1 1 1 1 1 1 1 1 


and >| xz converges so slowly that Dal x:| diverges, shows that the series 
involved in the problem may fail to be absolutely convergent. 

Also solved by V. F. Cowling, A. E. Livingston, and the Proposer, all of 
whom make explicit use of the properties of Banach spaces and refer to S. Ba- 


nach, Theorie des Operations Lineaires, Warsaw (1932), p. 67. 
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EpiTeEp By E. P. Vance, Oberlin College 
All books for review should be sent directly to E. P. Vance, Oberlin College, Oberlin, Ohio. 


Calculus, A Modern Approach, Second Enlarged Edition. By Karl Menger, The 
Bookstore, Illinois Institute of Technology, 1953. 


Of books on calculus there is no end. Yet seldom, if ever, does one of them 
make a new and significant contribution to the eternal problem of presenting 
the fundamental ideas correctly and systematically and “at a level of clarity 
and simplicity comprehensible to college students of mathematics, science and 
engineering.” This book is exceptional; it is really new. And to some experienced 
teachers its newness may be disturbing, especially to those who until now have 
labored under the impression that existing textbooks are all right as they are. 

Professor Menger explains in the preface the circumstances which prompted 
him to write the book: “Most mathematicians consider calculus as a domain 
about which they know everything. Twenty years ago the author, all in all, 
concurred in this view. Later, teaching the subject repeatedly and striving for 
greater clarity of the conceptual framework, he began to realize that he had 
underestimated the difficulties. In the late 1940’s he ventured on a project to 
rescue calculus from one of the greatest dangers that may befall a discipline— 
the danger of petrification. A situation which rules out any essential progress 
can easily arise when specialists abstain from a critical study of bool:s on a sub- 
ject and at best skim the pages with an eye toward their possible use as texts 
for sophomore classes. In 1952, even after publishing the first edition of this 
book, the author felt that he should study calculus further.” 

The book is an outgrowth of a sequence of works by Professor Menger which 
have appeared since 1944, beginning with his Algebra of Analysis, and which 
have had as their aim a systematization and clarification of the foundations of 
analysis. These earlier works were not primarily intended for beginners. The 
present work is a textbook for beginners and for their instructors. In fact it is 
through its effect on this latter class that mathematical instruction will benefit 
most from the ideas which are here presented for the first time in a relatively 
simple manner and in textbook form. 

It is not surprising that a reform of this kind should present striking changes. 
Among these the most striking at first are those of notation, and the author 
states in justification that “ ... the principal objective . . . is rather a clarifica- 
tion of the basic ideas and of their applications. In pursuing this aim the author 
found (to his regret) that changes in the notation were inevitable. In a certain 
sense these changes are only minor.” It becomes clear after reading the book 
that the invention of the new notation was an essential step toward the clarifi- 
cation of the basic ideas and their applications and is thus amply justified. How- 
ever, it is likely that these changes will be regarded by many readers as by no 
means minor ones. 
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There is so much in the book that is new that it is difficult to give a fair and 
comprehensive review of its content in a limited space. At the risk of injustice 
it may be said, briefly, that the main objectives of the book are the following: 

(1) To clarify, in a manner which accords with modern points of view, the 
ideas underlying the basic terms: number, class, function, variable quantity (herein 
designated by V.Q.), domain of a function, range of a function, domain of a V.Q., 
range of a V.Q., the derivative (function) of a function, anti-derivative (function) 
of a function, integral (function) of a function, inverse (function) of a function. 

(2) To study the operations performed on functions; the algebraic opera- 
tions and that of substitution, and the basic operations of calculus—differentia- 
tion, integration—and their reciprocity laws. 

(3) To present a new and adequate notation, free from ambiguities and in- 
consistencies, with which to represent the various concepts enumerated in (1) 
and (2). 

(4) To provide an introduction to the application of these concepts to the 
physical sciences. 

It goes without saying that all of the elementary functions are adequately 
presented. 

Never before, in the reviewer’s opinion, have these objectives been pursued 
in an elementary book with such painstaking care or with as much enthusiasm 
and vigor. As the author implies, this is no book to be skimmed over. Nearly 
every page must be read with scrutiny. The careful reader will detect several 
misprints not included in the list of errata, but these are almost all of the harm- 
less variety or of the kind that will puzzle the reader to a point where he will 
be rewarded by a better understanding. 

In Chapter I the author attacks the problems of differentiation and inte- 
gration from the geometrical point of view. He defines a simple curve (graph) C 
as one which is intersected in at most one point by any line perpendicular to the 
X-axis. There is no Y-axis. The two problems are those of finding the slope of 
the curve C and the area under C. New notation is introduced at the outset: 


b 
f C (read: the integral of C from a to 6) = the area under C from a to b. 


DCx (read: the derivative of C at x) = the slope of C at x. 
The horizontal line at the altitude c is denoted by c. Thus 


Dex = 0. 


The letter J denotes the line which, for any point x on X, has the altitude x 
above the point x. Thus 


Ix = &, Dix = 1, 
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If L is any non-vertical line, there are two numbers c and c’ such that the 
altitude of L above any point x on X is xc+c’. If this line is denoted by cI +c’ 
then 


+ 


for every x on X. 
If c is the slope line of L, then L is an area line of c, 


L is an ff if and omy if ¢ = DL. 


A line whose slope line or derivative is the horizontal line c is called an anti- 
derivative of c or Dc. 


L isaD~"c if and only if DL = ¢. 


Thus every fc is a D~c, and vice versa. 

The word function is not discussed in Chapter I, but the notation used for a 
simple curve (graph) is the same as that used later to represent the correspond- 
ing function. 

In Chapter II the study of the slope and area problems is continued by first 
introducing step lines and polygonal lines, the former being slope curves of the 
latter, and the latter being area curves of the former. Then the approximate area 
under a curve is found by drawing a step line, which approximates the given 
curve, and constructing the polygonal line which is the area curve of the step 
line. Similarly, the slope curve of the given curve is found approximately by 
means of an inscribed polygon whose slope curve gives approximately the slope 
of the given curve. 

In Chapter III the student is shown how to perform the calculations corre- 
sponding to the approximate graphical solutions found in Chapter II. The 
formula giving the approximate solution of the area problem is expressed in the 
form 


f Q = — + — + — 
where Q is a step line with steps at heights c,, cz - - - . This leads to the funda- 


mental inequalities for estimating the accuracy of the approximate solution, 
namely 


foas-fesfa 


where Q,, Qn, are step lines whose steps are at heights equal to those of the low- 
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oF est and highest points of C, respectively, in the m intervals. These ideas are il- 
lustrated by copious exercises and examples. 
x Chapter IV is entitled The Idea and the Use of Functions. This chapter could 


be read with profit by almost every teacher of elementary calculus. In it the au- 

thor, with admirable skill, leads the reader toward a thorough understanding of 

the function concept. He limits himself to functions whose domains are classes of | 
numbers and whose ranges are also classes of numbers. Thus he states, among 
many examples, “If, with each number, 16 times the square of the number is 
paired (e.g. with 0 the number 0, with 1 the number 16, with 3/2 the number 36), 
then mathematicians say that a function has been defined. More precisely, the 
class of the pairs of numbers so defined is called a function. Traditionally, these 
: definitions are expressed by saying that with each number x the number 16x? is 
paired. The function is described as the class of the pairs (x, 16x?) for all num- 
bers x.” 

. “The pairs (0, 0), (1, 16), (1.5, 36), (—1.5, 36), (4/2, 32) belong to the func- 

e tion, the pairs (16, 1), (3, 1/3), - - - do not, since 16-1671, 16-3341/3,---. 
ss Clearly the same concept can be defined by saying that with each number ¢ the 
ss number 16#? is paired, or by considering the class of pairs (a, 16a?) for all num- 
: bers a - - -. In each pair belonging to the function the first element is the num- 
* ber with which a number is paired. The class of all these first numbers is called 

; the domain of the function. The second element in each pair is the number which 
ie is paired with the first. Each such number is called a value of the function... . 
ie The class of all values of a function is the range of the function... .” 

“Clearly each function is a class of pairs of numbers. Two functions are 
equal if each pair of numbers belonging to the first function also belongs to the 
second, and each pair belonging to the second function also belongs to the 
first. Thus the function of which the pairs are (x, 16x?) for every x is equal to 
the function whose pairs are (¢, (44—1)(4#+1)+1) foreveryt....” 

‘ “But not each class of pairs of numbers is a function. In order that a class 
| of pairs of numbers be a function it is necessary and sufficient that in no two 
different pairs of the class are the first elements equal. ... ” 

“Can a function be equal to a number? Would anyone confuse a party of 
married couples with a spinster?” 

In this long chapter the basic properties of the elementary functions are 
studied with unusual care: the constant functions, the functions J, J*, polyno- 
mial functions, the exponentialand logarithmic functions, trigonometric functions 
and their inverses. A great deal of space is devoted to the difficult general con- 
cept of the inverse, in which the author makes a clear distinction between pairs 
of functions f, f* which are mutually inverse and those in which f* is inverse to 

f but f is not inverse to f*. This involves careful study, mainly geometrical, of the 
2 domains and ranges of functions, and is worth while. But the reader will be 
* puzzled at first by the proposal, at the very end of the chapter, to define the 
a class of elementary functions by means of the categories 
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I. (a) The constant functions c; 
(b) expo, logio; 
(c) arc r=2 arc tan, r=tan (J/2); 
(d) (m=1, 2,3,--+). 

II. Functions which can be obtained from those in (1) by the fundamental 

operations: addition, multiplication, division, and substitution. 

This is elegant but somewhat sophisticated for the beginner. He will have 
to read on in order to discover that J* can be proved to be an elementary func- 
tion in the domain [>0], for every c, and will have to use this result to prove 
for himself that the function arc sin is elementary: 


I 
arc sin = arc “( ). 


In fact the author gets the wrong answer. 

Chapter V, On Limits, gives an easy introduction to the theory of limits. 
The author illustrates, in a common sense manner, all of the essential ideas 
without actually using deltas and epsilons. He explains precisely what is meant 
by writing for example 


if y~ 2 and s~ 6, then y+2 ~ 8; 

suf suf arb 
which is read: “if y is sufficiently close to 2 and z is sufficiently close to 6, then 
y+z is arbitrarily close to 8.” Using this notation he makes it easy for the stu- 
dent to understand what is meant by 


lim (x? — 1) = 3, lim (2? — 1) = 3, 
22 
or, in his notation, 
(I? — 1) = 3, 
and, in general, 
lim f = b. 


The application of these ideas to the study of continuity is well done and 
the author takes the opportunity to identify “simple curves” with the graphs of 
continuous functions. 

The author’s intention here, as in other chapters, is to make the concepts 
clear and to state the facts correctly without supplying always the details of a 
complete proof. For example the number e is investigated by careful computa- 
tions of the values of the expression 4=(1+x)/* for small values of x. From 
which, with the use of tables of logarithms, it is then made clear that logie lies 
between .434 and .435. 
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Chapter VI, The Basic Concepts of Calculus, deals with derivatives, anti- 
derivatives, and integrals, by utilizing the analysis developed in Chapter V. 
Here the fundamental theorems of integral calculus are explained and made 
plausible. Applications involving elementary functions are plentiful and the ex- 
ercises are good. The chapter ends with a fairly thorough discussion of the Reci- 
procity Laws: 


L@ f f=f (f continuous), 


Law | is proved. Law II is not proved at this point because the Theorem of 
the Mean, needed in the proof, does not appear until Chapter X. The reviewer 
believes that the postponement of this basic theorem is unwise as it is essential 
to the study of the increase and decrease of functions and of problems in max- 
ima and minima, a subject which has important applications, and which, on ac- 
count of its appeal to the student, should come as early as possible in the course. 
It is noticeable that the author introduces Law II with a misleading statement 
to the effect that it will be proved in a later chapter for any function f having a 
derivative Df. Of course the proof, in Chapter X, is based on the continuity of 
Pf. 

In Chapter VII, Applications to Physical Science, the author treats the diffi- 
cult concepts of mensuration, variable quantity, and function, as applied to 
physical sciences. It is not easy to give an adequate idea of the effectiveness 
with which these concepts are developed. By easy stages, starting with a class = 
of elements ¢, o’, - -- of undetermined nature, which can be compared with 
each other by means of the usual order and equivalence relationships ¢ <0’, ¢~o’, 
there is established a pairing of the elements of 2 with numbers of a numerical 
scale so that with each element o of 2 a number no is paired—just as in the 
definition of function—but such that if o<o’ then noSno’, and no=no’ if and 
only if s~o’. By postulating the additive property for the elements of 2 and the 
existence of a unit element a; (an element with which the number 1 is paired), 
there is established a mensuration scale. Thus, if e~o,+01+0;, then with o is 
paired the number 3 and one writes ¢ = 30,; and if ¢+o0+o0=0,+0;, then with 
o the number # is paired and one writes ¢ = 30. And by an obvious extension the 
statement ¢ = +/20;, is defined. The ideas are copiously illustrated of course. The 
notion of the unit element is illustrated, for example, by such common state- 
ments of mensuration as: 

“The rod R is 7.5 ft.,” “the weight P is } lb.,” 
which yield the corresponding equivalent formulae: 


(a) R= 7.5 ft., P =} bb. 


Each such formula is an equality of an observable and a so-called denominate 
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number (a pure number followed by an observable comparable with the first). 
Menger then defines a preferred notation, which consists in replacing o’ = ko by 
(o’ in a) =k. So that, in the example, (a) above is replaced by: 


(b) (R in ft.) = 7.5, (P in Ibs.) = 3. 
By easy stages we then arrive at the Fundamental Law of Mensuration: 
If ine) =a 
and if (oc ino’) = 5, 
then ina) = ab. 


The author is now prepared to answer the question “What is a variable as 
studied in Science?” The answer is: 

“If a number is paired with each element of a class >—the same number with 
any two equivalent elements—then we call the class of pairs (element, number 
paired with the element) a variable, more specifically, a variable with the do- 
main 2. For example, if 2 is additive and a unit ¢; has been selected in 2, the class 
of pairs 


in for all elements o of 


is a variable with the domain 2. More generally, every numerical scale in a class 
2 of comparable elements gives rise to a variable, namely the class of pairs 
(c, no), for allo, where nga is the label of ¢ in the numerical scale. The numbers 
paired with the elements are called the values of the variable, and the class of all 
values of the variable its range.” 

Thus, for example, if in the class of rods having length, a foot is selected as 
the unit, for each rod p we have the pairing, (p, p in ft.), and the class of all 
such pairings, for all rods p, is an observable variable, called the length of rods. 
The variable is denoted by the letter /. Its values are /p = (p in ft.). 

Functions are the only variables whose domains are numbers. From this it 
follows that into a function a variable can be substituted provided that the range 
of the variable is a part of the domain of the function. 

The chapter culminates in a penetrating exposition of the central question: 
When is a variable m a function of another variable n? 

The fundamental idea comes from the definition of functionally related 
classes: Let 2 and T be two classes. A class of pairs (r, 7) such that r belongs to 
T and a belongs to Z is called a pairing of 2 with T. 2 is functionally related to T 
with regard to the pairing II if any two elements of 2 which are paired with two 
equivalent elements of T are equivalent; in other words, if whenever (7, 7) and 
(r’, 0’) belong to II and r~7’ then o~o’. 

This gives the central theorem: If the class 2 is functionaliy related to the 
class T with regard to the pairing II, and numerical scales are established, pair- 
ing labels mo and nr with all elements of = and T respectively, then the varia- 
ble m is a function of the variable n with regard to the pairing II; that is to say, 
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m=fn for some function f. [There is an unfortunate misprint in the text, not 
noted in the Errata; the result is given as n =fm. | 

The rest of this long chapter is devoted to the study of derivatives and inte- 
grals of functions connecting variable quantities. To illustrate: Let s be the 
distance travelled by a particle (in units of length), ¢ the time elapsed since the 
beginning of the observation of the motion (in units of time), s, the function 
connecting s with ¢. If 7, 7’, are the instants of time at which ¢ has the values x, 
x’, respectively, a, a’, the positions corresponding to those instants of time, so 
that 


So = Sir, so’ = 
then the difference quotient 


so’ — so sit’ — — 


— tr tr’ — tr x’ — 


is the average speed of the particle between x and x’ (that is, between the in- 
stants 7 and 7’). If this quotient has a limit as x’—<x then this limit, which is 
Ds.x, is called the instantaneous speed of the particle at x (that is, at the in- 
stant 7). In symbols, 
so’ — se — 
Ds.x = lim ———— = lim ————— 
tr’ — tr vor — 


is the instantaneous speed at x. Using v to represent the variable quantity in- 
stantaneous speed of the particle, then 


= Ds; x, for every number x, 
or 
= 
Another way to describe this relation is to say that: 
If then » For example; 
If s,;= 167? then » = 


These are functions. The traditional way of expressing these ideas would be: 
If f is the function connecting s with ¢ then 


s = ft, v = Dft. 
These quantities are variables. In particular: 
if s = 16/ then v = 322. 


To illustrate the use of integrals the author considers, among other examples, 
the work done by a force p acting on a particle (the rectilinear case). If p, is the 
function connecting with the distance travelled s, we have 
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This is the work function from a on. The variable ,w is the work from a on. 


For example, if p,= then fa(—J) = —}/*?+4a? and 
aw = — $5? + 3a’. 


More generally, if « and v are two variables, v a function of u, then gwu= Ju 
is called the cumulation function of v with u from a on. The variable aw, the cumu- 
lation of v with u from a on, is connected with u by the cumulation function as 


follows: 
( f re) u. 


The rest of the book need not be discussed in detail. In fact the most impor- 
tant parts are those, already discussed, in which the fresh point of view of the 
author is set forth. In the remaining chapters the ideas presented in the first 
seven chapters are applied to the techniques of calculus which, in turn, are 
applied to problems in analytic geometry and mechanics. 

Chapter VIII, The Calculus of Derivatives, gives a thorough exposition of the 
techniques used in calculating the derivatives of functions expressed in terms 
of differentiable functions by means of the fundamental operations (addition, 
multiplication, taking the reciprocal, and substitution). 

The operation of differentiation is represented by the symbol D when ap- 
plied to a pure function. Df+g means (Df)+g, Df-g means (Df) -g, Dfg means 
(Df)g; the derivatives of f+g, f-g, fg are denoted by D(f+g), D(f-g), D(fg). 

One of the main features of this chapter is the illuminating discussion of the 
meaning of the various notations used in representing the derivative of a func- 
tion and the derivative of one variable with respect to another (with which the 
first is functionally related). 

Chapter IX, The Calculus of Antiderivatives, is a lucid exposition of methods 
of finding elementary antiderivatives. It affords thorough training in the basic 
techniques and in the use of the Menger notation. 

Chapter X, The Mean Value Theorem and Its Consequences, begins with a 
common sense treatment of the approximate value of a function near a point a 
by means of its value at a and the value of the derivative at a. This is followed 
by a Mean Value Theorem, which is proved by means of Rolle’s Theorem which, 
in turn, is rendered plausible by a clear geometrical argument. The chapter cul- 
minates in Taylor’s Theorem with the Lagrange form of the remainder. The rest 
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of the chapter deals with the usual applications to computation with due atten- 
tion to estimations of the remainder term. 

Chapter XI, Two-Place Functions or Functions of Two Variables. This final 
chapter introduces new notation for the study of functions of several variables. 
Capital letters are used to denote two-place functions and all one-place functions 
are denoted by lower case letters. Here J is the two-place function such that 
I(x, y) =x for every (x, y), J is the two-place function such that J(x, y) =y for 
every (x, y). The one-place function denoted in the preceding chapters by /, 
such that Jt=¢ for every ¢ is now denoted by +. Graphs (surfaces) of simple two- 
place functions are studied. The idea of substituting one-place and two-place 
functions in a two-place function is given careful treatment. This idea is applied 
to the differentiation of one-place functions defined implicitly. Partial differ- 
entiation of functions is expertly presented but only the simpler techniques are 
treated in the text proper. The differentiation of a one-place function obtained 
by substituting two one-place functions in a one-place function is the only case 
of the Substitution Rule considered. There is an interesting but brief Appendix 
devoted to the study of partial differentiation of m-place functions and a state- 
ment of the general Substitution Rule. 

There is little doubt that this book deserves the most careful attention of 
teachers of calculus. It is written in the lively, vigorous style which is character- 
istic of the author and on this account it should appeal to any earnest student. 
Opinions among teachers may vary widely as to the seriousness of the defects of 
the traditional point of view of which the author complains. But there can be 
little difference of opinion as to the effectiveness of the remedies which he 
supplies. The fact is that even an experienced teacher who reads the book with 
care will probably find that he has acquired a deeper insight into the significance 
of the fundamental ideas and a better understanding of the art of presenting 
them to beginners. It must not be dismissed with the pointless remark that it is 
a book on function-theory, not calculus. The two are inseparable. But there is 
a debatable question: How deeply should one go into the foundations of analysis 
in writing a book for beginners? Some teachers prefer a book that skims over 
the surface of the subject because it affords them free opportunity to present 
the deeper ideas in the class room according to their own tastes. This is not that 
kind of book. But it still leaves opportunities for the teacher to introduce sup- 
plementary material, since much that is traditionally included in a calculus 
book is omitted; for example, drill exercises and problems of a so-called practical 
nature drawn from geometry, mechanics, etc. 

Unfortunately the typography of the book is unattractive, but this fault 
will undoubtedly be remedied by an enterprising publisher who will print the 
book in the attractive format in which so many lesser books are now being of- 
fered. 


H. E. Bray 
The Rice Institute 


“¢ 
4 
q 


1954] RECENT PUBLICATIONS 


NEW BOOKS RECEIVED 


History of the Theories of Aether and Electricity. By Sir Edmund Whittaker. 
New York, Philosophical Library, 1954. xi+319 pages. $8.75. 

College Algebra. New Third Edition. By P. K. Rees and F. W. Sparks. New 
York, McGraw-Hill Book Company, Inc., 1954. xiii+460 pages. $4.25. 

Analytic Geometry. Second Edition. By E. S. Smith, Meyer Salkover and 
H. K. Justice. New York, John Wiley and Sons, Inc., 1954. xiii+306 pages. 
$4.00. 

Understanding College Algebra. By E. R. Smith, Samuel Selby and Murray 
Kleiman. New York, The Dryden Press, 1954. xvi+571 pages. $3.50. 

Elements of Statistics. By H. C. Fryer. New York, John Wiley and Sons, 
Inc., 1954. vi+262 pages. $4.75. 

An Introduction to the Calculus of Finite Differences. By C. H. Richardson. 
New York, D. Van Nostrand Company, Inc., 1954. v+142 pages. 

Theory of Equations. By C. C. MacDuffee. New York, John Wiley and Sons, 
Inc., 1954. vii+120 pages. $3.75. 

Tables of Integral Transforms. By A. Erdélyi. New York, McGraw-Hill Book 
Company, Inc., 1954. xx+391 pages. $7.50. 

Matter Energy Mechanics. By Jakob Mandelker. New York, Philosophical 
Library, 1954. ix+73 pages. $3.75. 

First Course in Calculus. By H. R. Cooley. New York, John Wiley and Sons. 
Inc., 1954. xii +643 pages. $6.00. 

Introductory College Mathematics. By Adele Leonhardy. New York, John 
Wiley and Sons, Inc., 1954. ix+459 pages. $4.90. 

Introductory College Mathematics. By C. G. Jaeger and M. M. Bacon. New 
York, Harper and Bros., 1954. xii+382 pages. $4.75. 

The Mechanism of Economic Systems. By Arnold Tustin. Cambridge, Har- 
vard University Press, 1954. xi+161 pages. $5.00. 

Number. Fourth Edition. By Tobias Dantzig. New York, The Macmillan 
Company, 1954. ix+340 pages. $5.00. 

Grundzuge der Ausgleichungsrechnung. By Walter Grossmann. Berlin, Ger- 
many, Springer-Verlag, 1953. 261 pages. $4.56. 

The Japanese Abacus, its Use and Theory. By Takashi Kojima. Rutland, 
Vermont, Charles E. Tuttle Company, 1954. 102 pages. $1.25. 
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NEWS AND NOTICES 
EpITeD By EpitH R. SCHNECKENBURGER, University of Buffalo 


Readers are invited to contribute to the general interest of this department by sending 
news items to Edith R. Schneckenburger, University of Buffalo, Buffalo 14, New York. Items 
must be submitted at least two months before publication can take place. 


SALE OF MATHEMATICAL LIBRARY 


The library of the late Professor Otto Sz4sz is available for sale. It contains 
over 100 mathematical books of importance, complete sequences of the Bulletin 
of the American Mathematical Society, the Journal of the London Mathematical 
Society, and the Proceedings of the London Mathematical Society from 1933-1952, 
and other useful periodicals. To be mentioned also are the Proceedings of vari- 
ous mathematical congresses from 1912 on. In addition, about 4500 reprints are 
available. Interested persons may contact Mrs. Brigitta Dobratz, 2911 Knoll 
Drive, Concord, California. 


PERSONAL ITEMS 


The following have received Guggenheim Fellowships for 1954-55 for the 
studies indicated: 

Professor C. B. Boyer, Brooklyn College, studies in the history of the theory 
of the rainbow; 

Associate Professor D. G. Chapman, University of Washington, studies of 
the size and dynamic structure of wild animal populations; 

Professor S. S. Chern, University of Chicago, studies in the field of differen- 
tial geometry; 

Professor M. R. Hestenes, University of California, Los Angeles, study of 
the applications of the theory of linear spaces to the calculus of variations; 

Associate Professor E. R. Kolchin, Columbia University, study of certain 
questions in the Galois theory of differential fields; 

Professor C. C. Lin, Massachusetts Institute of Technology, studies in the 
field of theoretical fluid dynamics; 

Assistant Professor F. I. Mautner, Johns Hopkins University, research on 
unitary representations of locally compact topological groups; 

Professor Hans Rademacher, University of Pennsylvania, studies of gen- 
eralized modular functions belonging to real quadratic fields; 

Assistant Professor M. A. Rosenlicht, Northwestern University, study of 
non-complete group varieties; 

Research Professor A. H. Taub, University of Illinois, studies of the gravi- 
tational field in relativity and of curved shocks in hydrodynamics; 

Professor Alexander Weinstein, University of Maryland, study of initial 
and boundary value problems for singular hyperbolic equations and for partial 
differential equations of mixed type. 

National Science Foundation Postdoctoral Fellowships have been awarded 
to the following: G. E. Baxter, University of Minnesota; R. C. Blanchfield, 
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Princeton University; D. A. Buchsbaum, Columbia University; C. W. Curtis, 
Yale University; Jacob Feldman, Columbia University; W. M. Huebsch, Uni- 
versity of Notre Dame; A. P. Mattuck, Princeton University; J. C. Moore, 
Princeton University; Edgar Reich, University of California; M. F. Ruchte, 
University of Wisconsin; V. L. Shapiro, Institute for Advanced Study; Edwin 
Weiss, Massachusetts Institute of Technology. 

Professors B. W. Jones of the University of Colorado and Einar Hille of Yale 
University represented the Mathematical Association of America at the Inter- 
national Congress of Mathematicians which was held in Amsterdam, The Neth- 
erlands, September 2~9, 1954. 

Assistant Professor A. N. Aheart of West Virginia State College represented 
the Association at the inauguration of President W. J. L. Wallace of the College 
on April 9, 1954. 

Professor H. J. Ettlinger of the University of Texas was the representative 
of the Association at the inauguration of President D. H. Morgan of Texas 
Agricultural and Mechanical College on May 20, 1954. 

Professor Fritz John of the Institute for Mathematics and Mechanics, New 
York University, was appointed to represent the Association at the Second Uni- 
versity Convocation of the Bicentennial Year of Columbia University which 
was held on June 1, 1954. 

Brown University announces the following; Dr. John Wermer has been ap- 
pointed to an assistant professorship; Assistant Professor Joanne Elliott of Mt. 
Holyoke College has been appointed Visiting Assistant Professor for the year 
1954-55; Professor R. E. Gilman has retired with the title of Professor Emeritus. 
The University reports also that a program of instruction leading to the Bache- 
lor of Science in Applied Mathematics will be offered beginning Fall, 1954; Pro- 
fessor E. H. Lee, chairman of the Division of Applied Mathematics, is in charge 
of the program. 

At Montana State University: Dr. A. S. Merrill, dean of the faculty and 
professor of mathematics, has been appointed Vice-President of the University 
while retaining his other titles; Professor Harold Chatland, chairman of the 
Department of Mathematics, has been appointed Dean of the College of Arts 
and Sciences. 

Randolph-Macon Woman's College announces: Professor Evelyn P. Wiggin 
is on leave during the academic year 1954-55 and is studying at the University 
of Chicago; Miss Sara L. Ripy has been appointed to an instructorship. 

Rutgers University reports the following: Assistant Professor R. M. Cohn 
has been promoted to an associate professorship; Dr. B. H. McCandless, Dr. 
Solomon Leader, Dr. V. L. Shapiro, and Dr. K. G. Wolfson have been promoted 
to assistant professorships; Mr. N. R. Stanley is transferring from the New 
Brunswick Branch to the Newark Colleges of the University; Mr. A. G. Wootton 
has resigned. 

University of Connecticut announces the following appointments: Dr. Allen 
Devinatz of the Institute for Advanced Study to an assistant professorship; Dr. 
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R. L. Ingraham of the Institute for Advanced Study and Dr. H. C. Griffith of 
the University of Tennessee to instructorships. 

University of Illinois announces: Professor D. G. Bourgin has a Fulbright 
Grant to Rome for the academic year 1954-55; Professor S. S. Cairns has a Ful- 
bright Grant to the University of Strasbourg in France and a sabbatical leave 
for the year 1954-55; Assistant Professor Irving Reiner is on sabbatical leave 
during 1954-55; Assistant Professor Lowell Schoenfeld is on sabbatical leave 
during 1954-55 and is dividing his time between the University of Pennsylvania 
and the Institute for Advanced Study; Dr. Alex Heller of Harvard University 
has been appointed to an assistant professorship. 

At the University of Wisconsin: Associate Professor R. C. Buck has been 
promoted to a professorship; Assistant Professor W. F. Eberlein has been pro- 
moted to an associate professorship; Dr. R. C. T. Smith of New England Uni- 
versity College, New South Wales, has been appointed Visiting Lecturer for the 
academic year 1954-55; Assistant Professor Elizabeth Sokolnikoff Hirschfelder 
has resigned. 

University of Wichita reports the following: Assistant Professors Harold 
Huneke and Ferna Wrestler have been promoted to associate professorships; 
Miss Jeneva J. Brewer has been promoted to an assistant professorship; Miss 
Ann Klein has been appointed to an instructorship. 

Wayne University announces the following appointments to instructorships: 
Mr. W. S. Bicknell, Mr. Charles Briggs, Dr. Shu-Teh C. Moy, and Dr. E. S. 
Northam. 

Associate Professor B. H. Bissinger of Lebanon Valley College has been 
elected Chairman of the Department of Mathematics. 

Mrs. Elizabeth B. Bockelman, formerly with the Acoustics Research Labo- 
ratory of Harvard University, has been appointed to an instructorship at 
Wellesley College. 

Professor H. K. Brown of Northeastern University has been promoted to the 
position of Dean of the Graduate Division. 

Associate Professor D. E. Christie has returned to Bowdoin College after a 
year’s leave of absence as a visiting Fellow at Princeton University. 

Mr. Albert Derin of Kansas State College is now in Military Service. 

Assistant Professor D. G. Duncan of the University of Arizona has accepted 
an appointment at San Jose College. 

Mr. J. J. Gehrig of the Ballistics Research Laboratories, Aberdeen Proving 
Ground, has accepted a position as research engineer with Northrop Aircraft, 
Incorporated, Hawthorne, California. 

Mr. M. N. Haller, Jr., previously an engineer with Southern Bell Telephone 
and Telegraph Company, is now engaged as an engineer by Capehart-Farns- 
worth Company, Fort Wayne, Indiana. 

Mr. P. W. M. John, who has been teaching at the Casady School, Oklahoma 
City, has been appointed to an instructorship at the University of Oklahoma. 

Mr. H. T. LaBorde of the University of North Carolina has been appointed 
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to an assistant professorship at the University of the South. 

Dr. H. N. Laden has accepted a position as senior methods research officer 
with the Chesapeake and Ohio Railway Company, Cleveland, Ohio. 

Mr. P. M. Landry, formerly with the Vitro Corporation, Fort Walton, Flor- 
ida, is now in Military Service. 

Dr. H. M. MacNeille, who has resigned as Executive Director of the Ameri- 
can Mathematical Society effective September 1, 1954, has been appointed 
Professor and Chairman of the Department of Mathematics of Washington 
University. 

Associate Professor R. R. Otter has returned to the University of Notre 
Dame after a year as Visiting Professor at Princeton University. 

Assistant Professor L. J. Paige, who spent the year 1953-54 at the Institute 
for Advanced Study, has returned to the University of California. 

Associate Professor T. S. Peterson, Portland State Extension Center, has 
been promoted to a professorship. 

Mr: S. L. Pollack of the National Bureau of Standards has accepted a posi- 
tion as mathematician with the Raytheon Manufacturing Company, Waltham, 
Massachusetts. 

Dr. Edgar Reich of the Rand Corporation is spending the academic year 
1954-55 at the Institute for Advanced Study as a National Science Foundation 
fellow. 

Mr. Yomei Sawanobori, formerly at the Cornell Aeronautical Laboratory, 
Buffalo, New York, has a position as a mathematician with the I. B. M. Corpo- 
ration, New York City. 

Dr. Alice T. Schafer has been appointed to an assistant professorship at Con- 
necticut College. 

Assistant Professor M. R. Spiegel of Rensselaer Polytechnic Institute has 
been promoted to an associate professorship. 

Dean G. W. Starcher of the College of Arts and Sciences of Ohio University 
has been appointed President of the University of North Dakota. 

Assistant Professor B. B. Townsend of Louisiana State University has been 
promoted to an associate professorship. 

Professor J. F. Wampler, previously professor of mathematics and registrar 
at York College, has been appointed to an associate professorship at Nebraska 
Wesleyan University. 


Dr. R. C. Briant, project director at Oak Ridge National Laboratory, died 
on April 25, 1954. He had been a member of the Association for twenty-seven 
years. 

Professor Emeritus H. S. Brown of Hamilton College died on March 27, 1954. 
He was a charter member of the Association. 

Professor J. A. Daum of Texas Agricultural and Mechanical College died on 
March 23, 1954. 

Professor Emeritus Arnold Dresden of Swarthmore College died on April 10, 
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1954. He had served as President of the Association during 1933-34, as Vice- 
President during 1931, and as a member of the Board of Governors for the pe- 
riods 1935-40 and 1943-45. 

Professor Emeritus William Findlay of McMaster University died on Mav 7, 
1953. He was a charter member of the Association. 

Professor D. L. Holl, head of the Department of Mathematics of Iowa State 
College, died on May 22, 1954. He was a member of the Iowa Section for nine 
years. He served as Section Chairman in 1950-51 and had been elected Governor 
of the Iowa Section for the period 1953-56. 

Professor E. L. Post of City College of the City of New York died on April 
21, 1954. He had been a member of the Association for thirty-one years. 

Mr. J. E. Sanders, retired meteorologist of Knoxville, Tennessee, died on 
April 8, 1953. He was a charter member of the Association. 

Professor J. M. Stetson, William and Mary College, who was a charter mem- 
ber of the Association, died on May 5, 1954. 


THE MATHEMATICAL ASSOCIATION OF AMERICA 


Official Reports and Communications 
NEW SECTIONAL GOVERNORS OF THE ASSOCIATION 


The following have been elected Governors of the Association for a three- 


year term beginning July 1, 1954 by a mail vote of the membership of the Asso- 
ciation in the Sections indicated: 


Allegheny Mountain Morris Ostrofsky, Westinghouse Electric Corp. 
Indiana Ralph Hull, Purdue University 

Kentucky H. H. Downing, University of Kentucky 
Metropolitan New York R. M. Foster, Polytechnic Institute of Brooklyn 
Nebraska M. A. Basoco, University of Nebraska 

Northern California W. H. Myers, San Jose State College 

Oklahoma L. W. Johnson, Oklahoma A. and M. College 
Rocky Mountain C. A. Hutchinson, University of Colorado 
Wisconsin R. C. Huffer, Beloit College 


Over 50% of the members eligible to vote actually cast their votes in the 
Kentucky Section and the Wisconsin Section. In general, a high percentage of 
votes is cast in the elections of Sectional Governors, thus indicating the desire 


of members of the Association to have their Sections represented by outstanding 
candidates. 


H. M. Geuman, Secretary-Treasurer 
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THE 1954 COMBINED MEMBERSHIP LIST 


Each year the Mathematical Association of America and the American 
Mathematical Society issue a Combined Membership List. The 1954 edition of 
this List will be sent to members of the Society. 

The office of the Association should be notified promptly of all changes in 
rank, position, and address which have not previously been reported. The final 
date for receipt of changes is October 15. Any errors in the 1953 Membership 
List should also be reported before October 15. 


H. M. Geuman, Secretary-Treasurer 


NEW MEMBERS 


Professor H. M. Gehman, Secretary-Treasurer, announces that the following 
89 persons have been elected to membership by the Board of Governors on ap- 


plications duly certified. 


S. S. ABHYANKAR, A.M.(Harvard) Teaching 
Fellow, Harvard University. 

A. G. AsHoox, A.B.(Boston C.) Teacher, 
Oliver W. Holmes Jr. High School, Dor- 
chester, Massachusetts. 

SHELDON BALK, Student, Arizona State College, 
Tempe. 

G. E. Bartow, Jr., B.S.(V.P.I.) Assistant 
County Agent, Gate City, Virginia. 

G. C. Barton, B.A.(U. of Washington) 
Teacher, Concrete Senior High School, 
Washington. 

R. S. BEarD, M.S.(Kansas) Col. U.S.A. re- 
tired; Part-time Prof., Yeshiva Institute of 
Mathematics. 

C. P. Benner, M.S.(C.I.T.; Houston) Asst. 
Prof., University of Houston. 

H. F. Bennett, B.S. (Illinois Wesleyan) Pa- 
tent Agent, Eastman Kodak, Rochester, 

RENE Btocu, Student, University of Basel, 
Switzerland. 

J. D. Bucknottz, Student, Little Rock Junior 
College. 

REv. Fr. FLORENT CARTUYVELS, Lic. in Math. 
(Louvain) Teacher, St. Xavier's College, 
Ranchi, India. rg 

C. P. M.Ed. (St. Thomas) _Instr., 
Hermantown High School, Duluth, Minn. 

J. R. Cox, Student, Lebanon Valley College. 

J. B. Crawrorp, B.S. (U.S. Military Academy) 
Asst. Prof., Alabama Polytechnic Institute. 

C. L. Davis, M.A.(Michigan) Instr., General 
Motors Institute, Flint, Michigan. 


D. E. Dean, B.S.(Akron) Technical Super- 
visor, Firestone Tire & Rubber, Akron, 
Ohio. 

T. P. Denneny, M.S.(Notre Dame) _Instr., 
John Carroll University. 

R. S. Dick, Student, Queens College. 

A. Doucury, B.S. (Northeast La. S.C.) 
Mathematician, David Taylor Model 
Basin. 

B. M. Drucker, Ph.D.(North Carolina) Asst. 
Prof., Georgia Institute of Technology. 

G. E. Duncan, Student, Georgia Institute of 
Technology. 

G. A. Easterry, B.S.(Allen) Teacher, Mor- 
ristown College High School, Tennessee. 

E. J. FARRELL, M.A.(Stanford) Asst. Prof., 
University of San Francisco. 

M. A. FELDsTEIN, Student, Arizona State Col- 
lege, Tempe. 

C. T. Fixe, Student, University of the South. 

L. R. Forp, Jr., Ph.D. (Illinois) Research In- 
structor, Duke University. 

L. R. Foy, M.S.(St. John’s) Instr., Marian 
College, Poughkeepsie, N. Y. 

J. D. Gates, A.B.(Hiram) Teacher, Maple 
Heights High School, Ohio. 

J. B. GorBEL, Student, University of Oregon. 

S. W. Gotoms, M.A.(Harvard) Graduate 


Student, Harvard University; Instr., Bos- 
ton University. 

T. E. HaGensee, M.S.(DePaul) Teacher, 
Crane Technical High School, Chicago, 
Ill. 

D. I. Hammer, M.A.(Columbia) Instr., Mont- 
clair State Teachers College. 
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W. G. Hanks, Student, University of Mani- 
toba. 

M. E. Harris, Student, Yale University. 

H. L. Henpricxs, B.A.(Pomona) Teacher, 
Pasadena City College. 

Prec. R. A. Huttguist, M.S.(Purdue) Re- 
search Assistant, Fort Sill, Oklahoma. 
Rev. L. E. Isenecker, S.J., Ph.B. (Loyola) 

Graduate Student, Catholic University. 

DiaANE M. JouNsoNn, Student, University of 
Manitoba. 

Vireinta L. JoHNsoNn, Student, Harpur Col- 
lege. 

Mrs. CATHERINE S. Kay, M.A.(Michigan) 
Instr., North Central College, Naperville, 
Illinois. 

C. S. KrmsBati, Tech. Representative, Philco 
Corporation, Philadelphia, Pa. 

J. E. Kist, M.S.(Purdue) Research Assistant, 
Purdue University. 

S. M. Kiernan, Brooklyn, New York. 

Diane L. Kocuer, Student, Harpur College. 

A. R. LamMontaGne, M.A.(New Hampshire) 
Graduate Assistant, University of New 
Hampshire. 

ADELE LEONHARDY, Ed.D.(Missouri) Chair- 
man of Dept., Stephens College. 

K. R. Lucas, B.S.(Washburn) Graduate 
Student, University of Kansas. 

DALE Mangss, M.A.(Vanderbilt) Asst. Prof., 
Baylor University. 

M. D. Marcus, A.B.(California) Research 
Assistant, University of California at 
Berkeley. 

H. J. Jr., M.S. (Illinois) Training 
Instr., Meteorology, Dept. of Weather 
School, Chanute AFB, III. 

McKay, B.S. (Drexel) Instr., Drex- 
el Institute of Technology. 

D. E. McOwen, M.A. (Michigan) Instr., Gen- 
eral Motors Institute, Flint, Michigan. 

B. T. MENDELSOHN, Student, City College of 
New York. 

A. A. MuLLin, Student, Syracuse University. 

E. R. Muttuins, Jr., Ph.D.(Illinois) Asst. 
Prof., Swarthmore College. 

W. E. NANcE, Student, University of the South. 

J. H. Ortman, B.A.(Wisconsin) Mathemati- 
cian, Douglas Aircraft Co., Santa Monica, 
California. 

G. K. OverHOLTzER, Ph.D.(Indiana) Asst. 
Prof., Georgia Institute of Technology. 
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J. L. Pattone, A.B.(Fordham) Asst. Actu- 
ary, Woodward & Fondiller, New York 
City. 

N. C. Perry, Ph.D.(U.S.C.) Asst. Prof., 
Alabama Polytechnic Institute. 

W. L. Puitures, Jr., Student, Yale University. 

T. B. Pitts, B.A.(Texas) Instr., University of 
Texas. 

S. D. RasusHKa, Student, Washington Univer- 
sity. 

P. A. REICHLE, M.S.(Pennsylvania) M.A. 
(Columbia) Col., U.S.A., retired Asst. 
Prof., Lenoir Rhyne College. 

Mrs. Susan L. RED, Student, Queens College. 

Mrs. Joy B. Russex, M.S.(New York Univ.) 
Instr., University of Buffalo. 

N. B. Smita, M.A.(Oregon S. C.) Graduate 
Student, Iowa State College. 

LawrENcE Soxotorr, B.A.(Brooklyn) Jr. 
Mathematician, Republic Aviation Corp. 

J. M. Srarx, S.M.(M.I.T.) Graduate Stu- 
dent, Massachusetts Institute of Technol- 


ogy. 

R. F. Sremuart, M.A.(Columbia) Instr., 
Montclair State Teachers College. 

W. A. Stock, M.A.(Northwestern) Instr., 
Purdue University. 

C. R. Srrarn, M.S.(Purdue) Teaching Fel- 
low, University of Michigan. 

C. M. Sruart, M.A.(Duke) Asst. Prof., 
Clemson College. 

A. C. SuGar, Ph.D. (California) Consultant, 
Ryan Aeronautical Co., San Diego, Calif. 

F. B. Taytor, A.M.(Columbia) Asst. Prof., 
Manhattan College. 

L. C. TENG, M.A.(Michigan) D.I.C. Staff, 
Massachusetts Institute of Technology. 

E. A. TraBant, Ph.D.(C.1.T.) Assoc. Prof., 
Purdue University. 

C. K. Tsao, Ph.D.(Oregon) Asst. Prof., 
Wayne University. 

J. J. URBANcEK, M.A.(Northwestern) M.S. 
(DePaul) Chairman, Chicago Teachers 
College. 

W. A. VEzEAu, Ph.D. (St. Louis) Assoc. Prof., 
St. Louis University. 

G. E. Voyzes, B.A.(North Texas S. C.) 
Teacher, Decatur Baptist College. 

E. A. WALKER, M.A.(Sam Houston) Gradu- 
ate Student, University of Kansas. 

T. C. WaLtace, Student, Arizona State Col- 
lege, Tempe. 
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IRENE S. WELNA, Student, St. Joseph College. OswaLp WyLer, Sc.D.(Swiss Federal Inst. of 


Roy WEsTWICK, Student, University of British Tech.) Asst. Prof., University of New 
Columbia. Mexico. 

E. F. WaittLesey, A.B.(Princeton) Instr.. Monica J. WyzaLexk, Student, Harpur College. 
Bates College. H. J. ZassenHAus, Ph.D.(Hamburg) Prof., 


McGill University. 


THE MARCH MEETING OF THE METROPOLITAN NEW YORK SECTION 


The thirteenth annual meeting of the Metropolitan New York Section of the 
Mathematical Association of America was held at St. John’s University, Brook- 
lyn, New York, on March 27, 1954. Professor H. F. Fehr, Collegiate Vice-Chair- 
man of the Section, presided at the morning session and Professor W. H. Fager- 
strom, Chairman of the Section, presided at the afternoon session. 

One hundred two persons attended the meeting, including the following 
seventy-five members of the Association: 


R. G. Archibald, C. Y. Bartholomew, Jonas Beraru, Samuel Borofsky, C. B. Boyer, A. B. 
Brown, W. F. Cassidy, G. B. Charlesworth, K. P. Clancy, Charles Clos, P. J. Cocuzza, T, F. Cope, 
Demetrios Counes, W. H. H. Cowles, I. A. Dodes, J. N. Eastham, W. H. Fagerstrom, H, F. Fehr, 
J. M. Feld, A. B. Finkelstein, A. D. Fleshler, William Forman, R. M. Foster, G. C. Francis, 
D. H. Frank, E. T. Frankel, E. J. Germino, William Gonzalez, Bernard Greenspan, Harriet 
Griffin, W. T. Hamilton, C. M. Hebbert, G. C. Helme, E. Marie Hove, Aida Kalish, G. W. Kays, 
L. S. Kennison, H. S. Kieval, A. E. Kinney, M. S. Klamkin, Charles Koren, Edna Kramer-Lassar, 
C. H. Lehmann, D. R. Lintvedt, D. M. MacEwen, May H. Maria, Emanuel Mehr, F. H. Miller, 
Morris Morduchow, A. J. Mortola, Eugene Odin, Maria T. Pan, J. J. Quinn, J. K. Reckzeh, 
R. M. Reed, H. D. Ruderman, J. P. Russell, John Salerno, Charles Salkind, A. H. Sarno, Seymour 
Schuster, Abraham Schwartz, Aaron Shapiro, Edward Shapiro, James Singer, Geraldine D. Smith, 
Morris Smith, E. R. Stabler, Mildred M. Sullivan, R. L. Swain, P. M. Treuenfels, R. M. Walter, 
Alan Wayne, M. E. White, R. C. Yates. 


The following officers were elected for the year 1954-55: Chairman, Pro- 
fessor H. F. Fehr, Teachers College, Columbia University; Collegiate Vice- 
Chairman, Professor A. B. Brown, Queens College; High School Vice-Chair- 
man, Mr. D. H. Frank, Forest Hills High School; Secretary, Professor E. Marie 
Hove, Hofstra College; Treasurer, Mr. Aaron Shapiro, Midwood High School, 
Brooklyn. 

At the business meeting, the following report on the activities of the Com- 
mittee on Contests and Awards for the year 1952-53 was given by its chairman, 
Professor W. H. Fagerstrom: 

The fourth annual contest was held on May 14, 1953. There were 476 schools 
registered for the contest. Of these, 291 were in the Metropolitan New York 
area, and the other 185 were distributed throughout 30 states and provinces. 
Among the schools in this latter group were those in the units operated by Brit- 
ish Columbia, Oregon and Western N. Y. each conducting its own state 
wide contest using the contest questions of the Metropolitan New York Section. 

The maximum score available for the school was 450 points, and the maxi- 
mum for the individual was 150 points. The three highest ranking schools were: 
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James Madison High School, Brooklyn, N. Y., score 349; Phillips Acad- 
emy, Andover, Mass., score 336; Talmudical Academy, New York, N. Y., score 
325. The three highest ranking students were: E. Alan Phillips, Weston High 
School, Middlesex, Mass., score 130; Jerrold Rubin, James Madison High 
School, Brooklyn, N. Y., score 126; Isaac J. Sharon, Talmudical Academy, 
New York, N. Y., score 126. There were 28 recipients of Certificates of Merit. 
These certificates are awarded to the highest ranking student in each area, pro- 
vided his score was 85 or more. Honor keys, given in recognition of the students 
having won the award for two consecutive years, were awarded to 13 students. 

The Committee’s receipts for the year were $1687.33. Its expenses were 
$1680.84, leaving a balance of $6.49. 


Very Reverend John A. Flynn, C. M., S. T. D., President of St. John’s Uni- 


versity, welcomed the people at the meeting, and then the following papers were 
presented: 


1. Mathematics in communication, by Dr. Brockway McMillan, Bell Tele- 
phone Laboratories, Murray Hill, New Jersey. (By invitation.) 


The technology of communication is divided broadly into the fields of transmission (the actual 
conveying of messages) and switching (the setting up of transmission paths, when needed, out of 
equipment common to many possible paths). Most transmission is accomplished by devices gov- 
erned by linear time-invariant differential equations. The relevant applicable mathematics is then 
centered around the Laplace and Fourier transforms. Switching devices are discrete and quantized; 
their relevant mathematics is that of logic, general algebra, and combinatorics. In both fields one 
distinguishes between the study of apparatus per se and the study of its performance in the en- 
vironment of use. The latter study always invokes probability theory in some form. Almost all 
communication problems invoke linear graphs in their statement or study. 


2. Fundamental preparation in mathematics for college study in engineering 
science—What the colleges need and expect, by Professor F. H. Miller, Cooper 
Union. 


Topics, concepts and techniques needed by the engineering student for college work in mathe- 
matics were discussed. Some of the items are essential prerequisites; others are highly desirable if 
the requirements of physics and engineering departments are to be adequately met. Specific 
examples of various topics and applications of techniques were cited, reference being made to 
A Report on Mathematics Preparation for Engineering Colleges, by F. H. Miller and S. G. Roth, 


published in the Journal of Engineering Education, April, 1947, and to other reports bearing on the 
subject. 


3. Fundamental preparation in mathematics for college study in engineering 
science—What the high school can do, by Mr. D. H. Frank, Forest Hills High 
School. 


Difficulties in meeting the needs of the engineering school are: 1) Too little mathematics 
taken too long before entering college, 2) too many non-essentials taught, 3) lowering of standards, 
4) requirements the same for all students, 5) small numbers entering engineering resulting in diffi- 
culty in separating them from the others, 6) interest in science and engineering not necessarily ac- 
companied by interest in mathematics. 

The ways to meet these problems are: 1) separation in classes or segregation within classes, 
2) integration where possible, 3) rewriting of syllabus so that each student is exposed to all the 
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branches of mathematics for every term of the subject, 4) discontinuation of solid geometry as 
separate subject and distribution over all grades, 5) weeding out all non-essentials, 6) teaching 
elements of differential and integral calculus in last year with stress on all necessary concepts and 
techniques, 7) at all times teaching for the necessary disciplines so well exemplified in mathe- 
matics. 


4. On informal symbolic logic and its place in mathematical education, by Pro- 
fessor E. R. Stabler, Hofstra College. 


The speaker outlined some key topics for an informal treatment of symbolic logic. These 
included propositions, propositional functions, quantified propositions, and classes. Then he pre- 
sented examples of the potential usefulness of these topics in promoting the objectives of high 
school and college mathematics courses. Finally, considering the possible role of informal logic in 
the curriculum, he proposed: 1) an experimental twelfth grade course in logic and logical aspects of 
mathematics; 2) incorporation of some informal symbolic logic in college mathematics courses for 
general students; 3) a basic college course including informal symbolic logic, and related founda- 
tional topics, for students majoring in mathematics. 


E. Marte Hove, Secretary 


THE MARCH MEETING OF THE MICHIGAN SECTION 


The annual meeting of the Michigan Section of the Mathematical Associa- 
tion of America was held on March 27, 1954, at the University of Michigan, 
Ann Arbor, in conjunction with the meetings of the Michigan Academy of Sci- 
ence, Arts and Letters. Professor J. S. Frame of Michigan State University 
presided at both morning and afternoon sessions and at the luncheon and busi- 
ness meeting. 

A total of one hundred thirty-two persons attended the meetings, including 
the following sixty-six members of the Association: 

A. C. Andersen, K. J. Arnold, J. L. Bagg, J. W. Baldwin, R. C. F. Bartels, F. A. Beeler, 
J. H. Bell, W. S. Bicknell, J. W. Bradshaw, Fred Brafman, C. H. Butler, A. T. Butson, C. D. 
Calhoon, Y. W. Chen, R. V. Churchill, D. F. Coffey, S. D. Conte, J. W. Coy, C. C. Craig, D. A. 
Darling, P. S. Dwyer, C. M. Erikson, N. C. Fisk, K. W. Folley, J. S. Frame, J. W. Gaddum, 
Casper Goffman, G. W. Grotts, V. G. Grove, Frank Harary, Gerald Harrison, G. E. Hay, H. M. 
Heater, Fritz Herzog, T. H. Hildebrandt, E. F. Ingalls, L. G. Johnson, L. S. Johnston, P. S. Jones, 
Leo Katz, A. F. Lampen, J. F. Lanahan, Harry Langman, H. D. Larsen, K. B. Leisenring, J. F. 
Manogue, J. E. McLaughlin, L. E. Mehlenbacher, G. E. Meike, D. M. Mesner, E. E. Moise, 
W. K. Moore, H. W. Nace, A. L. Nelson, P. A. Nurnberger, Mary H. Payne, Emily C. Pixley, 


H. H. Pixley, J. H. Powell, P. H. Raker, H. F. Stelson, Brother Andrew Stephen, B. M. Stewart, 
J. G. Sowul, Leonard Tornheim, R. L. Wilder. 


At the business meeting the nominating committee consisting of Professors 
H. D. Larsen, B. M. Stewart, and L. E. Mehlenbacher proposed Professor R. V. 
Churchill, University of Michigan, for chairman and Professor S. D. Conte, 
Wayne University, for Secretary-Treasurer for the year 1954-55. The slate was 
elected unanimously. 

A motion was made by Professor R. V. Churchill and seconded by Professor 
V. G. Grove that the By-Laws of the Section be changed, if necessary, to allow 
the Sectional Governor to be a member of the Executive Committee of the Sec- 
tion. The motion was passed. 
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It was voted to empower the Executive Committee of the Section to make 
such resolutions as are necessary to effect the deposit and withdrawal of funds 
from banks. 

The following papers were presented at the morning and afternoon sessions: 

1. The enumeration of connected linear graphs: an analogy to rooted trees, by 
Professor Frank Harary, University of Michigan. 

Using a combination of Pélya’s counting theorem and a functional identity involving ex- 
ponentials and the cycle indexes of the symmetric groups, it is easy to show that the method of 


Cayley for enumerating rooted trees is identical in form with a derivation of the number of con- 
nected linear graphs in terms of the known number of all (connected or disconnected) linear graphs. 


2. Distribution of binomial coefficients modulo p, by Mr. D. M. Mesner, Mich- 
igan State College. 

The set of (n?+-n)/2 binomial coefficients for exponents 0, 1, - - - , #—1 is considered. For p 
a prime, the distribution of this set into residue classes modulo p is obtained. If m has the form pt 


the distribution is quite simple; in particular, the combined frequency of non-zero residues is 
[(p?+p)/2]*. If fa(é) denotes the frequency of numbers =i (mod p), then 


(1) lima = 0; 


and 


It is shown that (1) holds even when is replaced by an arbitrary modulus m. 


3. An elementary proof of an old number theorem, by Professor Harry Lang- 
man, Detroit Institute of Technology. 

The theorem states that it is impossible to have four different integral squares in arithmetic 
progression. The proof is entirely elementary and is based on the argument of “infinite descent.” 
It is shown that if a set of four integers exists whose squares in order have a common difference, 


another set exists conforming to the same conditions and such that the largest of these is less than 
the largest of the assumed set. 


4. Tschebyscheff polynomial approximations in high speed computing, Dr. Leo 
Razgunas, University of Michigan, Willow Run Research Center. 


Presented by title. 


5. The angle sums in a simplex, by Dr. J. W. Gaddum, Michigan State Col- 
lege. 

In vol. 59, no. 6 of this MonTHLY, bounds were found by the speaker for the angle sums (di- 
hedral and trihedral) in a tetrahedron. In this paper, similar bounds are discussed for the dihedral 


angle sum in a simplex. The bounds on a trihedral angle in terms of the surrounding dihedral angles 
are also discussed. 


6. Teaching of mathematics at Russian universities, by Professor G. G. 
Lorentz, Wayne University, introduced by the Secretary. 


Programs, teaching methods and the organization of Russian Universities were discussed. 
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7. Some mathematical problems in industry, by Dr. Morris Ostrofsky, West- 
inghouse Electric Corporation, Pittsburgh. (By invitation.) 


Several industrial problems in mathematics were discussed, two of them in some detail. The 
first of these involved the numerical solution of a heat conduction equation du/dt = DV*u+heat 
source where the heat source was non-linear in both time and space coordinates. It was pointed 
out that numerical solutions can be used as an experimental tool to crystallize physical concepts. 

The second problem involved the solution of the integral equation known as the transport 
equation. In course of that solution numerical integrations in a form convenient for digital com- 
puters were necessary. A useful and simple formula for calculating the coefficient for numerical inte- 
gration was presented. (In conclusion, the opportunities for mathematicians in industry were 
discussed.) 


8. Solving problems on Wayne’s UDEC (Unitized digital electronic computer), 
by Dr. E. P. Little, Wayne University, introduced by the Secretary. 


The Wayne University UDEC (Unitized Digital Electronic Calculator built by the Burroughs 
Corporation) has magnetic drum storage and five hole teletype input and output. Pulse duration 
is 0.1 sec. and pulse frequency is 120KC/sec. The drum stores 5300 words of 9 decimal digits plus 
sign and check number. Information within the machine is expressed in excess-three binary coded 
decimal notation. Numerical data and instructions are both stored on the drum. Random access is 
possible by a decision code but instructions are normally stored in sequence as far as possible. 
The tabulation of a simple linear function is described. 


9. On a probability distribution arising in group organization theory, by Mr. 
J. H. Powell, Michigan State College. 


In the very simplest model of a group organization we consider an organization of individuals 
in which connections in either direction exist or do not exist between pairs. Such an organization 
may be represented by an incidence matrix with 0’s on the principal diagonal. The immediate 
problem confronting us is that of the chance distribution of the number of mutual connections, 
i.e. pairs of symmetrically placed 1’s in the incidence matrix. In particular, consider the distribu- 
tion relative to a fixed set of row and column totals of the incidence matrix. Using a result by 
Katz and Powell in a paper on directed graphs, it is shown how to compute as many factorial 
moments of the exact distribution as desired and, implicitly, the exact probability of K mutual 
connections. 


10. Use of bordering in matrix inversion, by Professor P. S. Dwyer, Univer- 
sity of Michigan. 


This paper shows how bordering, which is frequently useful in evaluating determinants, may 
be applied to the problem of matrix inversion. Bordering is especially useful when relations involv- 
ing elements of the matrix to be inverted are present. The proof is based on the formula for the 
inverse of a partitioned matrix and the fact that each term of the inverse may be interpreted as the 
ratio of two determinants. Application is made to the problem of the size of the errors of the in- 
verse resulting from errors in the elements and to the determination of the inverse of a correlation 
matrix when the factors of the reduced correlation matrix are known. 


S. D. Secretary 


THE MARCH MEETING OF THE SOUTHEASTERN SECTION 


The annual meeting of the Southeastern Section of the Mathematical Asso- 
ciation of America was held March 19-20, 1954 at the University of South Caro- 
lina, Columbia, South Carolina. Professor W. V. Parker, Chairman of the Sec- 
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tion, and Professor W. L. Williams, Vice-Chairman, presided over the general 
sessions; Professors H. K. Fulmer, J. C. Currie, J. D. Novak, D. F. Barrow and 
L. A. Dye presided over subsections. 

There were about two hundred in attendance including the following one 
hundred sixteen members of the Association: 


P. L. Armstrong, Emil Artin, R. W. Bagley, D. F. Barrow, Helen Barton, R. C. Blackwell, 
R. G. Blake, A. T. Brauer, J. P. Brewster, G. M. Brown, J. W. Brown, B. F. Bryant, C. C. Buck, 
E. A. Cameron, N. A. Childress, R. S. Christian, D. H. Clanton, J. B. Coleman, J. A. Cooley, 
C. L. Cope, F. E. Cothran, R. W. Cowan, R. R. Croxton, J. C. Currie, N. E. Dodson, E. G. 
Douglas, B. M. Drucker, G. E. Duncan, L. A. Dye, E. D. Eaves, J. C. Eaves, D. O. Ellis, M. H. M. 
Esser, M. K. Fort, Jr., Tomlinson Fort, H. K. Fulmer, L. L. Garner, Leslie J. Gaylord, D. R. Good- 
ner, M. O. Gonzalez, S. T. Gormsen, B. F. Hadnot, S. W. Hahn, J. C. Harden, Jr., E. A. Hedberg, 
Marguerite Z. Hedberg, T. F. Hickerson, T. R. Norton, L. P. Hutchison, C. W. Huff, G. B. Huff, 
J. B. Jackson, F. W. Kokomoor, Stephen Kulik, G. B. Lang, J. W. Lasley, Jr., C. G. Latimer, T. H. 
Lee, Anne L. Lewis, F. A. Lewis,G. H. Lundberg, R. A. Lytle, E. L. Mackie, C. E. Martin, S. T. 
Martin, W. N. Mebane, Jr., J. T. Moore, R. H. Moorman, J. C. Morelock, J. A. Nohel, J. D. 
Novak, G. K. Overholtzer, W. V. Parker, P. B. Patterson, A. H. Payne, Lillian G. Perkins, N. C. 
Perry, C. G. Phipps, R. L. Plunkett, Ellen F. Rasor, P. A. Reichle, G. E. Reves, J. M. Robertson, 
H. A. Robinson, F. Virginia Rohde, H. R. Rouse, W. C. Royster, ‘’. A. Rutledge, F. W. Saunders, 
P. C. Scott, E. B. Shanks, Olivia H. Shanks, D. C. Sheldon, C. Eucebia Shuler, J. L. Sloan, 
C. B. Smith, F. W. Stallard, E. L. Stanley, L. W. Stark, C. M. Stuart, C. S. Sutton, H. S. Thurston, 
R. Z. Vause, Jr., J. H. Wahab, D. T. Walker, E. F. Ward, W. W. Weber, W. M. Whyburn, 
M. C. Wicht, W. L. Williams, N. K. Williamson, R. A. Willoughby, R. L. Wilson, F. J. Witt, 
J. W. Young, B. K. Youse. 


The following officers were elected for the coming year: Chairman, Professor 
W. L. Williams, University of South Carolina; Vice-Chairman, Professor R. H. 
Moorman, Tennessee Polytechnic Institute; Secretary-Treasurer, Professor 
H. A. Robinson, Agnes Scott College. 

The following program was presented: 

1. A gap in the treatment of dependent functions, by Professor C. G. Phipps, 
University of Florida. 


In treating dependent and independent functions a situation unaccounted for occurs where 
the Jacobian of the functions is only conditionally zero. This condition usually implies that one 
function has a maximum or a minimum subject to the restraint of fixed values for the others. 


2. Marginal notes, by Professor J. W. Lasley, Jr., University of North Caro- 
lina. 


Synthetic division of one polynomial by another, evaluations of a determinant, complex 
numbers as vectors, rotation of axes, the equations of conics, Newton’s correction formula, deriva- 
tives of vectors, systems of linear equations, rank of a matrix are topics for which special tech- 
niques are developed. 


3. Lissajous figures in the teaching of undergraduate mathematics, by Professor 
J. H. Wahab, Georgia Institute of Technology. 


A method of drawing Lissajous figures is adapted to the study of parametric equations, com- 
posite functions, and the solutions of differential equations with variables separated. 
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4. Mathematics in England, by Professor Tomlinson Fort, University of 
Georgia. 
Professor Fort spent the academic year 1952-1953 in Great Britain. He lectured at various 


British universities and had an unusual opportunity to see British mathematicians and to observe 
mathematical organizations and schools. His paper is concerned with his personal observations. 


5. A geometric device to facilitate the use of hyperbolic substitutions, by Pro- 
fessor J. W. Young, University of Florida. 


As a geometric device to facilitate the use of hyperbolic substitutions in integrals involving 
Vut+a?, /uv2—a’, /a?—x? and a rational function of sh ¢ and ch ¢, Dr. Young suggested one 
write \/a*—(iu)?, 1\/12?—a?, /a*+(iu)*, respectively and place results on a right triangle in the 
complex plane modified by taking the argument pure imaginary. The appropriate substitution 
follows by the Theorem of Pythagoras and the relations 


tan ; 
ont 


sin , cos , 
sin’ ch cos” th 


sh 
6. Summer conferences in collegiate mathematics, by Professor E. A. Cameron, 
University of North Carolina. 


Professor Cameron told of the two 1954 summer conferences and outlined the work of the 
one to be held at his University. 


7. Some notes on matric computation, by Professor W. V. Parker, Alabama 
Polytechnic Institute. 
This paper is devoted to a discussion of various techniques in computation with matrices with- 


out resorting to the use of determinants. Several well-known results are obtained by use of these 
techniques. 


8. Exchange of information. 


Representatives from 28 mathematics departments participated in a discussion on certain 
common administrative problems which had been suggested prior to the meeting. Professor H. K. 
Fulmer rendered a very helpful report to the section of the findings of this discussion group. 


9. Pliicker numbers and the field of values, by Professor J. C. Currie, Georgia 
Institute of Technology. 

Let A be a square matrix with complex elements. Then the complex numbers xAZ’ with the 
vector x subject to xz’ =1 fill up a convex region of the complex plane, called the field of values of 
A. The boundary of this region is an algebraic curve, or, more often, a part of an algebraic curve, 
some of whose properties are examined with the aid of the Pliicker relations. 


10. Criteria for a logarithmic solution of a certain type of linear differential 
equation of the second order with a regular singular point, by Professor R. W. 
Cowan, University of Florida. 


Using the method of Frobenius, a compact expression is obtained for a general coefficient of 
an infinite series solution of the differential equation. The criteria desired arise from a considera- 
tion of various values of the difference of the roots of the indicial equation. 
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11. On differential equations with variables separated, by Professors J. A. 
Nohel and G. K. Overholtzer, Georgia Institute of Technology. 


In most popular differential equations texts the problem of finding a solution of a differential 
equation with variables separated is treated as a mere term-by-term antidifferentiation process. 
In this teaching note it is pointed out that an existence theorem can be proved by elementary 
means. 


12. Vibrating nets, by Professor Tomlinson Fort, University of Georgia. 


Professor Fort shows how the study of a vibrating loaded net leads to the study of a boundary- 
value problem for a partial difference equation of the second order. The solution of these problems 
in their simpler forms depends upon a study of the characteristic roots of a matrix. In the more 
complicated cases progress has been made primarily by kinematical methods. In some cases 
methods similar to those used in the study of analogous problems in the theory of differential equa- 
tions have proved effective. 


13. An almost regular form, by Professor P. B. Patterson, University of Flor- 
ida. 


An almost regular form is defined; the form (2, 2, 5, —1, —1, 0) is given as an example; and 
it is shown that this form does not represent primitively any m* where m=1 (mod 4). The method 
of approach is essentially that used on indefinite forms by B. W. Jones and E. H. Hadlock in a 
joint paper presented to the Society by E. H. Hadlock on November 29, 1952. 


14. The Frenet formulas for a ruled surface, by Professor J. D. Novak, Uni- 
versity of South Carolina. 


Extending Aoust's notion of courbure inclinée V.G. Grove introduced the curvature matrix and 
the Frenet Formulas for a ruled surface with respect to a curve. In this paper the “first and second 
curvatures” are used to study ruled surfaces, particularly the line of striction, and to characterize 
certain classes of ruled surfaces. 


15. On quasi-idempotent matrices, by Professor G. B. Huff, University of 
Georgia. 


Professor Huff defined a matrix A to be a quasi-idempotent if there exists a matrix polynomial 
F(x) such that A’ = F(r) for every positive integer r. It was then shown that a given matrix A isa 
quasi-idempotent if and only if there is an integer k such that A(E—A)*=0 and that the associ- 
ated F(x) is readily computed for a quasi-idempotent A. The explicit formulas lead to a miniature 
theory of exponentials and logarithms of certain special matrices. 


16. The numerical quartic, as a synthesis within elementary mathematics, by 
Professor F. E. Cothran, Greenwood, South Carolina. 


The method of Ferrari indicates parabolic and “hemi-quartic” loci which when “normalized” 
produce a graphical solution of the resolvent cubic. The devices used were elementary but lead the 
student to an interest in numerical calculus, the nature of real numbers, combinations of roots, 
difference and differential calculus, and hyperbolic functions, thus synthesizing several compart- 
ments of mathematics. 


17. A criterion for primary crossed products of degree p* over fields of formal 
power series in four variables, by Professor J. T. Moore, University of Florida. 


The paper provides a test for the primary nature of a crossed product algebra of the type 
indicated, which probably can be extended to apply to a much more general situation. 
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18. A note on unstable homeomorphisms, by Professor B. F. Bryant, Vander- 
bilt University. 

Utz (Unstable Homeomorphisms, Proc. Amer. Math. Soc., vol. 1, 1950, pp. 769-774) calls a 
homeomorphism f of a compact, dense-in-itself, metric space X onto X unstable, with instability 
constant, 5, provided for each pair of points x, yCX there exists an integer n(x, y) such that 
d(fr(x), f"(y)) >&. It was established that if w(x) [a(x) ] does not contain periodic points then there 
exists u, v@ w(x) [a(x)] whose orbits are negatively [positively] asymptotic, and with the aid of 
this it was shown that there exist p, g, r, x©X such that the orbits of p and g are positively asymp- 
totic and the orbits of r and s are negatively asymptotic. 


19. The theory of braids, by Professor Emil Artin, Princeton University. 


The theory of braids gives a good example that illustrates problems of topology and can also 
be used as an example to the concept of abstract groups. The composition of braids is done in 
the obvious way by tying one to the other. It can be shown that this leads to a well defined group 
generated by especially simple braids with several relations between them. The problem of classi- 
fication amounts to a solution of the word problem in this group. The solution of the word problem 
can be given in two ways—an algebraic and a geometric way. 


20. An integral transform, by Professor B. F. Hadnot, Florida State Univer- 
sity. 

The integral = [o*+!/ T(k+1) for ¢>0 is proved to occur as the limit of 
certain (R, ¢, x) sums of }-a,e~*». The integral J(e) is viewed as an integral transform of Ad*(w) 
for the limit e—0. One use of this idea has been to simplify certain known results in the theory of 
Riesz and Abel summation. 


21. Elliptic integrals in terms of Legendre polynomials, by Professor M. O. 
Gonzalez, University of Alabama. 


For the elliptic integral of the first kind u= fei-k sin? ¢)~"%d¢ the author obtains the 
expansion 


= (—1)* 


where P,(A) is the Legendre polynomial of order m and \=k’?—k?. 
It is also shown that 


= (—1)* 
and that 
(- 1)**1 


nm — 1)(2m + 1)(2n + 3) Pe 


E= + Prs(d)] = (a). 
@ 


Expansions for other elliptic integrals are also given. 


22. Necessary and sufficient conditions for convergence of infinite series, by 
Professor E. B. Shanks, Vanderbilt University. 


In this paper, the writer proves that most known tests for convergence or divergence of infinite 
series of positive terms can be deduced simply from the following two theorems. 

A necessary and sufficient condition that a series }-a, of positive terms converge (diverge) is 
that there exist real numbers C, (an unbounded sequence of real numbers D,) such that 
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for some fixed and n=1, 2,-- +. For example, Kummer’s test for convergence follows by 


taking C,=a,D,/p, p>0, v=1 in the first theorem. Numerous other standard tests are derived 
similarly in the paper. 

The theorems are also of theoretic importance as well as of practical importance, unlike most 
theorems similar in nature. 


23. Summation of series, interval of differencing one half, by Professor F. Vir- 
ginia Rohde, University of Florida. 
Given a function ¢(x) whose finite integral is f(x), interval of differencing h=}. If all differences 


of $(x) are zero after the kth one, an expression for , andl ¢(x) can be found in terms of f(x), o(x), 
and the k differences of ¢(x) taken between suitable limits. 


24. A mathematics program for high speed computation, by Professor B. M. 
Drucker, Georgia Institute of Technology. 
Dr. Drucker announced that his institution had received a grant for the purchase of a new 


electronic digital computer. He outlined the types of courses which will be offered in the fall for the 
preparation of programs, coding and the operation of an automatic computer. 


25. On definitions of trigonometric functions, by Professor W. A. Rutledge, 
Alabama Polytechnic Institute. 


The standard text-book definitions of the trigonometric functions as functions of angles are 
considered and some of the difficulties that arise, both technical and pedagogical, are discussed. 
Methods of avoiding these difficulties are presented. 


26. The proper approach to trigonometry, by Professor Tomlinson Fort, Uni- 
versity of Georgia. 


Professor Fort has developed the subject of trigonometry with the basic idea that the infinite 
series is a concept with which the student is already informally familiar through his work with 
infinite decimals, infinite geometric series, infinite binomial series and numerous other limits in his 
elementary geometry. A review of these topics is proposed, still strictly in an informal manner, 
to be followed by a discussion of functions and the definition of sin x and cos x by infinite power 
series. Applications to the solution of triangles appear as one of many uses. 


27. On the least primitive root (mod p), by Professor Alfred Brauer, Univer- 
sity of North Carolina. 

In a paper presented at the meeting of the American Mathematical Society at East Lansing, 
Michigan, on September 4, 1952, the following theorem was proved. 

Let p be a prime of form 4n+1 and k the number of different prime divisions of p—1. We set 
2*=r. If g is the smallest positive primitive root (mod p), then 
(1) 


In this paper the following theorem is proved which often gives a better estimate than (1). 
Let s=s(p) be the maximal number of consecutive integers which are all not relatively prime 
to p—1. Then 
(2) g< (st1), 


The bounds (1) and (2) are greater than the well known bound of Vinogradoff, but are obtained 
by purely elementary methods. 
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28. Finitely generated rings, by Professor Emil Artin, Princeton University. 


A ring R is called Noetherian if its ideals satisfy the ascending chain condition. Let R be a 
Noetherian ring with unit element, S an extension ring with the same unit element such that 


S = Ruy + Rur+ +++ + Ron 


for certain elements w, of S. Such a ring S will be called a modul finite extension of R. 

If on the other hand one can find a finite number of elements £1, &, - - - , & of S such that 
every element of S can be expressed as polynomial of £1, &, - - - , & with coefficients in R then we 
call S a ring finite extension of R. The following theorem is proved: 

Let RC TCS such that S is a ring finite extension of R but a modul finite extension of 7, then 
T is a ring finite extension of R. 

As an application one can show: A ring finite extension of a field which is a field, is a finite alge- 
braic extension. 


29. Some inequalities concerning Mills’ ratio, by Professors D. F. Barrow and 
A. C. Cohen, Jr., University of Georgia, presented by Professor Barrow. 

Mills’ ratio is obtained by dividing the area under the normal curve to the right of any point 
by the ordinate of the normal curve at that point. The authors developed a method by which 
certain inequalities involving Mills’ ratio can be more easily proved than by existing methods; 
and one new inequality was proved which had been only surmised before. This method depends 
upon the following obvious theorem: 

If F(x) is everywhere positive, and if the product F(x)f(x) has a non-negative limit as x 
grows infinite and a derivative which is everywhere negative, then f(x) is everywhere positive. 


30. A normalizing transformation of the fiducial distribution of point-biserial 
correlation, by Professor N. C. Perry, Alabama Polytechnic Institute. 

Lev has shown that ry»; (point biserial correlation) has a distribution function related to non- 
central Student’s distribution. Johnson and Welch have pointed out a normal curve approxima- 


tion to the non-central ‘t’ frequency curve. The present paper combines these two results, obtain- 
ing a useful normal approximation to a transformation of roi. 


31. Moisture stresses in a certain type of plywood, by Professor C. B. Smith, 
University of Florida. 

Certain types of plywood plates having the directions of the grain of adjacent plies not 
perpendicular are considered. Moisture stresses produced by a uniform change in moisture content 


are discussed mathematically, and it is shown that by a suitable choice of the angles formed by the 
directions of the grain of adjacent plies warping can be avoided. 


32. Abridged series for numerical evaluation, by Professor B. K. Youse, Uni- 
versity of Georgia. 


This paper has appeared in this MONTHLY, March, 1954, page 184. 


33. Function space approximation in laminated orthotropic material under 
plane strain, by Professor R. G. Blake, University of Florida. 


The method of approximation in elasticity based on the concept of function space is used to 
obtain a first approximation to the solution of certain problems of plane strain in a rectangle con- 
sisting of three layers of orthotropic material cemented together. 

This material was prepared under a contract with the Office of Naval Research. 
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34. A Cremona invoiutorial transformation without fundamental curves of the 
Sirst species, by Professor L. A. Dye, The Citadel. 


The transformation involves a homoloidal web of cones with collinear vertices on a funda- 
mental line of contact of the second species. The principal surfaces are images of isolated funda- 
mental points on this line of contact. 


35. On the non-equivalence of two indefinite ternary forms, by Professor T. R. 
Horton, University of Florida. 


L. E. Dickson's Studies in the Theory of Numbers furnished a table of reduced, indefinite, non- 
zero, ternary quadratic forms of determinant d$83 except for d=68, d=81. Dickson stated 
(p. 147) that of four possible reduced forms of d =68, the equivalence of two forms A and B was 
not decided. This paper gives a necessary and sufficient condition that a form be equivalent to a 
form of a certain type and applies this condition to show A and B not equivalent. 


36. A theorem concerning mapping of a topological group into the circle, by Pro- 
fessor R. L. Plunkett, Vanderbilt University. 


A factor group of the group of all mappings of a compact, connected, commutative topological 
group, G, satisfying the second axiom of countability, into the circle is shown to be isomorphic 
with the character group of G. This is shown to imply that every mapping of such a group into 
the circle is homotopic to an interior mapping. 


37. The mapping of a second system of algebraic curves invariant under a cyclic 
involution of period eleven, by Professor S. T. Gormsen, University of Florida. 


In this paper the second of a series of systems of algebraic curves, invariant under an Ji, is 
mapped onto a surface F in S;. It is found that these curves at one of the branch points on F 
have five tangent lines at this point. Four of these tangent lines belong to a fifth order tangent 
cone to F, and the fifth tangent line belongs to the remaining tangent plane to F at this point. 


38. Distinction between some geometric properties arising in non-analytic map- 
pings, by Professor W. C. Royster, Alabama Polytechnic Institute. 


Interesting distinctions between convexity and starlikeness arise when mapping by non- 
analytic functions are considered. In this note the mapping function is applied to certain families 
of curves. 


39. On infinite series of sets, by Professor D. O. Ellis, University of Florida. 


Continuing a series of studies (Autometrized Boolean algebras, Boolean functions, Boolean 
metrized spaces, etc.) of similarities and contrasts between Boolean rings and the real field, one 
considers the convergence (in the Kantorovitch sequential topology) of }>;_, Ai in a set algebra, 
where addition is taken as symmetric difference. The main result, quite clear after a little thought, 
is: Do¢., As converges if and only if lim:.. A;=0. Certain simple consequences of this result are 
obtained. 


40. An alternative formalism for the basic notions of algebraic topology, by 
Professor C. C. Buck, University of Alabama. 


The peculiarities of this notation are that the boundary of an n-cell is a set of (—1)-cells, 
rather than a chain, and that an n-chain of a complex K over a commutative group G is written as 
a mapping of the n-cells of K into G. This notation exhibits clearly the relationship between the 
boundary chain and the coboundary of a cell. 
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41. On orbital topologies, by Professor R. W. Bagley, University of Florida. 
Two problems (Hausdorff property and continuity of the topology-inducing mapping) sug- 
gested by the paper, David Ellis, Orbital topologies, Quart. Jour. of Math. (Oxford), vol. 4, 1953, 


pp. 117-119, are solved. The topological structure is determined in the special case for which the 
topology-inducing mapping is a permutation. 


42. Open topological disks in the plane, by Professor M. K. Fort, Jr., Univer- 
sity of Georgia. 

A subset X of the plane has the “almost fixed point property” if corresponding to each con- 
tinuous function f of X into X and each positive number e there is a point p in X such that the 
distance from p to f(p) is less than e. It is shown that if an open topological disk in the plane 
is bounded and has a locally connected boundary, then the open disk has the almost fixed point 
property. 


43. Projections of surfaces obtained from an involution of period 13, by Pro- 
fessor N. A. Childress, University of Florida. 


The image of a planar cyclic involution of period thirteen can be represented as a surface of 
order thirteen in Ss. This paper concerns five successive projections of this image surface from a 
certain point on the surface being projected each time. Two results of a peculiar nature are ex- 
plained. 


H. A. Rosinson, Secretary 


THE APRIL MEETING OF THE IOWA SECTION 


The forty-first annual meeting of the Iowa Section was held jointly with the 
sixty-sixth annual meeting of the Iowa Academy of Science and the twenty-first 
convention of the Junior Academy of Science at Iowa State College, Ames, Iowa, 
on April 30 and May 1, 1954. The Chairman, Professor J. O. Chellevold, Pro- 
fessor H. T. Muhly, acting for the Vice-President, and the Secretary, Professor 
Fred Robertson, presided in turn. 

An amendment to the by-laws was approved by a unanimous vote. The 
amendment reads as follows: 

That Section 2 of Article II of the By-Laws of the Iowa Section of the Mathe- 
matical Association of America be replaced by: 

Section 2. The Executive Committee shall consist of the officers of the Sec- 
tion and the Sectional Governor. 

The registered attendance was sixty-eight including the following forty-four 
members of the Association: 


E. W. Anderson, Fred Brandner, R. H. Bruck, I. H. Brune, E. L. Canfield, J. O. Chellevold, 
Marian Daniells, A. W. Davis, W. M. Davis, C. B. Germain, B. E. Gillam, H. E. Goheen, Cornelius 
Gouwens, F. S. Harper, J. J. L. Hinrichsen, D. L. Holl, R. S. Jacobsen, N. W. Johnson, Jr., G. E. 
Kaldenberg, A. A. Karwath, Dora Kearney, Don Kirkham, O. C. Kreider, R. J. Lambert, C. E. 
Langenhop, B. F. Laposky, C. H. Lindahl, W. D. Lindstrom, F. W. Lott, C. G. Maple, J. V. Mc- 
Kelvey, H. T. Muhly, E. N. Oberg, Fred Robertson, Hazel M. Rothlisberger, J. A. Schumaker, 
Augusta L. Schurrer, E. R. Smith, R. D. Stalley, H. P. Thielman, Henry Van Engen, G. P. Weeg, 
D. V. V. Wend, Roscoe Woods. 
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Officers elected for the year are: Chairman, Professor H. T. Muhly, State 
University of Iowa; Vice-Chairman, Professor F. A. Brandner, Iowa State Col- 
lege; Secretary-Treasurer, Professor Fred Robertson, Iowa State College. 

The following papers were presented: 


1. Recent advances in the foundations of euclidean plane geometry, by Professor 
R. H. Bruck, University of Wisconsin. 


Unfortunately for the wide audience whose interest in geometry was awakened in high school 
or university, recent significant advances in the foundations of the subject are hidden behind a 
screen of abstract algebra. The present paper offers a pictorial account of some of the geometric 
axioms, a simple explanation of the algebraic problems which these pose and a brief account (with 
references rather than proofs) of the answers. The main topic is coordinatization of euclidean 
planes. Touched on are planar ternary rings, Veblen-Wedderburn systems, division rings with the 
right inverse property, right alternative and alternative division rings. 


2. Curve families and finite asymptotic values, by Professor D. V. V. Wend, 
Iowa State College. 


The author gives some relationships between the topological structure of a branched regular 
curve family F and the finite asymptotic values of the analytic functions contoured by F. Let a 
b-curve be an element of F separated from some other element of F by no element or chain of F. 
If F contains a b-curve, then any analytic function contoured by F has a finite asymptotic value. 
If the trees of F are isolated and F contains no b-curve, then there exists an analytic function 
contoured by F which has no finite asymptotic value. A precise characterization of the level curve 
families of harmonic polynomials is given. 


3. Peculiar derivative functions, by Mr. N. B. Smith, Iowa State College. 


Let f(x) be a function defined on an interval x in E;, and let p be a point property of f(x). The 
function f(x) is said to be peculiar with respect to the property » if there exists a partition of X 
into two subsets X; and X2 each everywhere dense in X and such that f(x) has the property p at 
no point of X. Peculiar derivative functions are discussed in the cases where the property ?p is 
continuity, neighborliness of W. W. Bledsoe, Proc. Amer. Math. Soc., vol. 3, 1952, p. 114, and 
cliquishness of H. P. Thielman, this MonTHLy, vol. 60, 1953, p. 156. In particular it is shown 
that there exist no derivative functions peculiar on X with respect to neighborliness. 


4. Mathematical theory for the utilization of tagged atoms in determining plant 
nutrient transformation rates in soils, by Professors Don Kirkham and W. V. 
Bartholomew, Iowa State College. 


Differential equations are derived and integrated, for determining mineralization rates and im- 
mobilization (demineralization) rates, in a tagged atom system of mineral and non-mineral plant 
nutrients undergoing simultaneous interchange in soil. The results are 


st b 


t 
s= + 
( a b ab 
where m= mineralization rate; i=immobilization rate; t=time; a=mass of tagged atoms, mineral 
and non-mineral; b=mass of all atoms, mineral, non-mineral, tagged and untagged; x = mass (at 


time ¢) of tagged and untagged mineral atoms; y=mass of tagged mineral atoms; xo and yo, 
values of x and y at ¢=0. The theory is in good agreement with experiments. 
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5. The gradient method for the Minkowski metric in the solution of simultane- 
ous linear algebraic equations, by Professors H. E. Goheen and J. R. Winkelman, 
Iowa State College. 


The system oi simultaneous equations 
Azx—c=0, 


in which A is a non-singular »Xm matrix and c is an mX1 matrix, can be solved by minimizing 
the positive definite function 


fix) 


in which ¢; is the ith component of the 2X1 matrix Ax—c. The authors have devised a method for 
obtaining this minimum. 


6. A modification of the Hartley-Politz scheme for dealing with not-at-home per- 
sons, by Professor Alan Ross, Iowa State College, introduced by the secretary. 


A method for obtaining unbiased estimates of population totals and estimates of their reli- 
ability is derived which does not require that callbacks be made on the not-at-homes in sample 
surveys. The theory of sampling with unequal probabilities of selection is expanded to cover the 
situation in which one utilizes a priori probabilities of obtaining schedules from all persons in the 
universe when these probabilities are actually determined only for those persons interviewed on 
first call. Horvitz and Thompson, Journ. Am. Stat. Assoc. vol. 47, 1952, p. 260, have considered 
sampling with arbitrary probabilities from the standpoint of a fixed sample size. The extension 
made here is to consider the case in which the sample size is fixed, but from the “sample” a vari- 
able number of schedules is obtained. Tables are presented which give comparisons of efficiency of 
the no-callback scheme with the usual survey technique involving repeated callbacks on the not-at- 
homes. 


7. The randomization analysis of a generalized randomized block design, by 
Professor M. B. Wilk, lowa State College, introduced by the secretary. 


The problems of estimation of effects and of tests of significance are considered for an experi- 
mental design in which ¢ treatments are allotted at random to r blocks, each containing pt experi- 
mental units, under the restriction that each treatment appears p times in every block. The sta- 
tistical inference is referred to the actual set of experimental units employed; a linear model is 
used, based on means and deviations over a finite population of (conceptual) yields. No distribu- 
tional assumptions are made. 

The means, variances, and covariances under randomization of the analysis of variance mean 
squares are given. Unbiased estimates of certain linear combinations of the parameters of the model 
are given, together with the variances of the estimates. Attention is also given to the estimation of 
these variances. Tests of significance, based on a randomization test, are presented. The possible 
approximation to the randomization test of procedures which derive from normal theory is con- 
sidered. 


8. An extension of some limit theorems in probability theory, by Professor John 
Gurland, Iowa State College. 


Some well-known limit theorems are proved and extended by using properties of sequences of 
characteristic functions. Special attention is given to stochastic convergence. 


9. Motivation in mathematics, by Miss Ethel Cain and Students of Roosevelt 
High School, Des Moines, introduced by the secretary. 


The speaker discussed motivation as a two-way path. The teacher may motivate the student 
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but the student may also motivate the teacher. Three students in the speaker's class, Mr. Jud 
Harper, Mr. Bob Helmick and Mr. Dick Rothrock, gave illustrations such as short hand multipli- 
cation, a method of obtaining certain combinations, a solution of a cubic equation by a system 
of bodies immersed in a liquid and controlled by a lever, and a set of geometric figures. 


10. A survey course for seniors, by Professor J. O. Chellevold, Wartburg 
College. 


In this paper the author described a senior-level course designed to show the relationship 
among various fields of mathematics and to provide a brief introduction to selected topics in modern 
mathematics. 


11. A study of student achievement and mortality in college mathematics, by 
Professor Fred Robertson, Iowa State College. 


The author discussed the achievements and mortality, as of the end of the winter quarter 1954, 
of the 1004 students who entered the college algebra course in September 1952. Success was defined 
as normal progress or a 2.0 or better grade point in each subject. An achievement and mortality 
table was shown for the group. 


FRED ROBERTSON, Secretary 


THE APRIL MEETING OF THE OHIO SECTION 


The thirty-eighth annual meeting of the Ohio Section of the Mathematical 
Association of America was held at the Ohio State University, Columbus, Ohio, 
on Saturday, April 17, 1954. Professor P. R. Rider, Chairman of the Section, 
presided at the morning and afternoon sessions. 

Eighty-nine persons registered in attendance, including the following sev- 
enty-five members of the Association: 


J. E. Adney, Jr., C. E. Amos, Grace M. Bareis, I. A. Barnett, H. M. Beatty, Foster Brooks, 
O. E. Brown, Emalou Brumfield, L. E. Bush, C. D. Calhoon, V. B. Caris, M. N. Chase, W. G. 
Clark, G. M. Clough, C. C. Crell, H. K. Crowder, R. C. Davis, B. B. Dressler, Wade Ellis, R. W. 
Emmert, P. L. Evans, H. E. Fettis, D. T. Finkbeiner, M. P. Fobes, H. W. Godderz, Marshall 
Hall, Jr., H. G. Harp, Frances Harshbarger, C. H. Heinke, R. G. Helsel, R. Y. Iwanchuk, L. A. 
Jehn, E. D. Jenkins, M. L. Johnson, Margaret E. Jones, John Kaiser, Chosaburo Kato, L. C. 
Knight, Jr., F. A. Kros, H. W. Kuhn, Nathan Lazar, R. F. Liskovec, L. L. Lowenstein, A. D. 
Martin, Margaret E. Mauch, S. W. McCuskey, E. J. Mickle, L. H. Miller, Knox Millsaps, Emma 
J. Olson, E. A. Peters, H. S. Pollard, Tibor Radé, P. V. Reichelderfer, P. R. Rider, R. F. Rinehart, 
D. L. Robb, G. deB. Robinson, Louis Ross, H. J. Ryser, Charles Saltzer, Samuel Selby, L. L. 
Shetler, R. L. Shively, Ruth B. Smyth, E. T. Stapleford, R. E. Thomas, Harold Tinnappel, 
H. S. Toney, W. R. Transue, D. R. Whitney, C. O. Williamson, J. A. Wilson, Alberta Wolfe, 
A. D. Ziebur. 


The following officers were elected for the coming year: Chairman, Professor 
W. R. Transue, Kenyon College; Secretary-Treasurer, Professor Foster Brooks, 
Kent State University; third member of the Executive Committee, Professor 
E. J. Mickle, The Ohio State University; Program Committee: Chairman, Pro- 
fessor Marion Wetzel, Denison University; Professor Wade Ellis, Oberlin Col- 
lege; Professor H. D. Lipsich, University of Cincinnati. 
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The following papers were presented: 
1. Statistical distributions, by Professor P. R. Rider, Wright-Patterson Air 
Force Base (Chairman’s address). 


Some of the statistical distributions in common use were described, methods of deriving them 
were discussed, and various applications were mentioned. 


2. Calculation of integrals of the form Jj sin? cost ® d®, by Mr. H. E. Fettis, 
Wright-Patterson Air Force Base. 


Attention is directed toward the evaluation of the integral 
f sin? cost ¢d@ 


where p> —1 if 0<x/2 and p>—1, g>-—1 if @=/2. Exact expressions are known to be possible 
when @=2/2 for arbitrary values of p and gq, and for integral values of p and gq, with @ arbitrary. 
The present discussion shows that reduction of elliptic integrals can be effected for p=+}, 
p=+4, g=0 and @S7/2. Also, a rapidly converging series is given for computation with values 
of p and q for which closed expressions are not possible. Several practical applications are noted. 


3. On the directional derivative, by Professor Wade Ellis, Oberlin College. 


The existence of the directional derivative df/ds of the function f =f(x, y) at the point (x, ) 
implies the existence of a plane x tangent to the surface S: z=f(x, y) at the point (x, y, z). The 
relation between df/ds and the dihedral angle formed by x and the z-plane can be exploited to en- 
hance the student’s intuitive grasp of this and related material in courses in the calculus. 


4. n-th roots of integral binary matrices, by Mr. Harvey Weitkamp, Univer- 
sity of Cincinnati, introduced by Professor I. A. Barnett. 


If X is a binary matrix, then X* may be written linearly in terms of X and I by the Cayley- 
Hamilton equation, and the coefficients may be given explicitly in terms of the trace and de- 
terminant of X. Necessary and sufficient conditions are given in order that a binary matrix with 
rational or integral elements shall have an mth root with rational or integral elements. The mth 
roots may be unique or may be infinitely many. 


5. Linear programming, by Dr. E. Leonard Arnoff, Operations Research 
Group, Case Institute of Technology. (By title.) 


Linear programming, through the transportation and simplex methods, presents recently de- 
veloped techniques for solving the problem of optimizing a linear functional, subject to restraining 
conditions which are in the form of linear inequalities. This paper first presents a brief description 
of the transportation and simplex methods. It then compares the types of optimization problems 
to which these methods are applicable with those optimization problems which can be handled by 
the methods usually taught in the college classroom, namely, those methods which use analytic 
geometry, differentiation, Lagrange multipliers, and the like. 


6. Functional mathematics for the secondary school teacher, by Professor Ema- 
lou Brumfield, Kent State University. 


The proposal is made that present mathematics-education students would be better pre- 
pared to teach secondary school mathematics if they could study, concurrently with their content 
mathematics courses, the importances, the applications, and the significance of these courses with 
respect to secondary school mathematics. The proposed study could be organized and effected by 
means of a mathematics-education course without affecting the content offerings of the department 
and without requiring the sectioning of students according to their professional objectives. 
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7. The place of algebra and geometry in the undergraduate curriculum, by Pro- 
fessor Gilbert deB. Robinson, University of Toronto. (By invitation.) 


The place of modern algebra in the undergraduate curriculum was discussed by Professor 
MacLane in Baltimore last Christmas; the author here discusses integrating the algebra with the 
appropriate geometry. This could be done by treating linear algebra and geometry in a one-term 
course, say in the junior year; the ideas of linear transformation, lattice of subspaces of a projec- 
tive space, determinants, matrices, the full linear group, efc., would all be coordinated. Quadratic 
ideas involving the reduction of a quadratic form and the application of Taylor’s series to a study 
of polarity and collineation, the inversion transformation and stereographic projection would fall 
naturally in the senior year. A one-term course in the first postgraduate year dealing with axio- 
matics, the significance of Pappus’ and Desargues’ theorems, the non-euclidean metric, the or- 
thogonal group, and the groups of the regular solids in three dimensions would tie the ideas to- 
gether. Such a three-term sequence could be extended and elaborated, but it would seem to provide 
a necessary minimum of geometrical background for a student proposing to go on to graduate 
work, and could replace existing courses in modern algebra and the foundations of mathematics. 


8. The convergence of sequences determined by a simple non-linear difference 
equation, by Professor Samuel Goldberg, Oberlin College. 

Let a real number xo be prescribed and determine x, x2, - - « by the difference equation %n41 
=x,?—k. Assume k= —1/4. The values of xo for which the sequence (x,) converges are found. In 
the special case k =2, these values form a denumerably infinite set and are determined from an 


explicit solution of the difference equation. The mode of convergence of (x,) to its limiting value is 
also discussed. 


9. Laplace transform and differential operator relations, by Dr. L. V. Robin- 
son, Wright-Patterson Air Force Base. 


Among other relationships, it is shown that 
1 a 
Lis] = =p (>, = =} 
and 
1 
where {'[f(¢)] is the Laplace transform of f(é). 


10. Some applications of the finite-difference analogue of Green's third identity, 
by Professor Charles Saltzer, Case Institute of Technology. 

The finite-difference analogue of Green’s third identity is deduced from the finite-difference 
analogue of Green's second identity (due to Courant, Friedricks and Lewy) by the use of the 
fundamental solution of the Laplace difference equation developed by McCrea and Whipple, Stohr 
and Duffin. This identity is applied to the Laplace, Poisson, and biharmonic difference equations. 


11. How mutation geometry says “No” to the trisection problem, by Dr. Beck- 
ham Martin, Owens-Illinois Glass Company. 


Two mutation views were developed: one starting from two prototypes resulting in the muta- 
tion trisection expression 


2b-rd-r = (pb — sd)-r 
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where 6 and d are the orthonormalized bisectors of the given angle and its supplement. The equa- 
tion was transmuted to a primordial prototype 


A-r=0, 
whose solution was 
r= 0/0. 


A second view started from a single prototype which led to a factorable 4th degree equation 
with a free choice parameter with which to control its irreducibility. A mutation curve was drawn 
for this which was shown to trisect the angle. 


Foster Brooks, Secretary 


THE APRIL MEETING OF THE ROCKY MOUNTAIN SECTION 


The thirty-seventh annual meeting of the Rocky Mountain Section of the 
Mathematical Association of America was held at Colorado Agricultural and 
Mechanical College, Fort Collins, Colorado, on Friday afternoon and evening 
and Saturday forenoon, April 30 and May 1, 1954. Professor M. L. Madison, 
Chairman of the Section, presided at all three sessions. 

Seventy-five registered for the meeting, including the following forty-nine 
members of the Association: 


C. F. Barr, D. L. Barrick, J. R. Britton, R. G. Buschman, R. K. Butz, F. M. Carpenter, 
A. G. Clark, Sarvadaman Chowla, G. S. Cook, Rev. F. T. Daly, W. E. Dorgan, H. T. Guard, R. R. 
Gutzman, C. L. Harbison, Leota C. Hayward, I. L. Hebel, LeRoy Holubar, J. E. Householder, 
P. F. Hultquist, C. A. Hutchinson, B. W. Jones, A. J. Kempner, Claribel Kendall, R. B. Kriegh, 
J. S. Leech, M. L. Madison, W. E. Mientka, M. W. Milligan, W. K. Nelson, Greta Neubauer, 
D. O. Patterson, O. M. Rasmussen, O. H. Rechard, A. W. Recht, L. W. Rutland, Jr., Nathan 
Schwid, S. R. Smith, W. N. Smith, L. C. Snively, K. H. Stahl, J. McD. Staley, P. O. Steen, 
E. P. Tovani, E. L. Vanderburgh, W. W. Varner, J. F. Wagner, F. J. Wall, C. R. Wylie, Jr., 
A. Zirakzadeh. 


Officers elected at the meeting for 1954-1955 were: Chairman, Professor 
Nathan Schwid, University of Wyoming; Vice-Chairman, Professor C. R. 
Wylie, Jr., University of Utah; Secretary-Treasurer, Professor F. M. Carpenter, 
Colorado School of Mines. 

The following papers were presented: 

1. On expressing the matrix A‘ as a polynomial in t, by Professor R. K. Butz, 
Colorado Agricultural and Mechanical College. 


This paper discusses the notion of E,F(X) matrices as introduced by G. B. Huff (Matrices 
such that A‘ is a polynomial in t and principal idempotent elements, Bull. Amer. Math. Soc., vol. 59, 
1953, p. 54). Emphasis is placed on the fact that the proofs of the main theorems require only the 
more elementary concepts of matrix theory and on methods of finding F(X) given a matrix A with 
elements in the field of complex numbers. The clarity with which some classical results follow by 
the use of this notion is pointed out. 


2. An optimum solution of N equations in M unknowns with N greater than M, 
by Mr. Leon Rutland, University of Colorado. 


A problem in engineering design led to a consideration of the system of equations 
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= 1,2,---,m), (n> m), 


j=l 


where the desired optimum solution of the system is that set of x’s for which the largest absolute 
value of any of the deltas in the set of equations 


be +O = 
j=l 


is as small as possible. Several theorems giving solutions to the problem under various conditions 
were either proved or stated with the proofs being omitted. A numerical example, carried through 
on a digital computer, was cited to indicate that even though there are as many as forty equations 
the method is feasible and the answer can be readily attained. 


3. The method of Frobenius, by Professor R. H. Cook, South Dakota School 
of Mines and Technology, introduced by the secretary. 

The usual textbook presentation of the method of Frobenius effectively camouflages two 
important points: (1) that the method involves a Taylor’s expansion, and quite often leads to just a 
Taylor’s expansion; (2) the conditions under which the method is applicable. This paper sug- 
gests a modified approach which has neither of the above disadvantages, is easily taught, and is 
sufficiently flexible to be applicable to many non-linear problems. 


4. Approximate solutions to a certain functional equation, by Professor C. A. 
Rogers, Colorado A and M. College. 

The following is investigated: Given a non-negative g(x), defined for all x >0, and which is 
strictly monotonic increasing and everywhere differentiable, to find a closed-form f(x), reasonably 


computable, such that f[f(x)] is at least approximately equal to g(x). It was indicated how this 
approximation could be accomplished for certain g-functions, with examples. 


5. Some infinite series, by Dr. W. E. Briggs, Professor S. Chowla, Professor 
(Emeritus) A. J. Kempner, and Research Assistant W. E. Mientka, University 
of Colorado, presented by Mr. Mientka. 


It is proved that 
= 23), 
1 
where 


6. Effect of rotation on the normal mode frequencies of transverse vibration of a 
cantilever beam, by Professor R. H. Cook and Mr. L. J. Eatherton, South Da- 
kota School of Mines and Technology, presented by Mr. Eatherton. 


The differential equation, 
ay ay ay, ay 
EI — —s= 


describes the vibrations, in a vertical plane, of a cantilever beam which rotates about a vertical 
axis through its clamped end. This equation and appropriate boundary conditions are considered 
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by the use of Taylor’s expansion. The normal mode frequencies are calculated in terms of Q, the 
rotational speed. Results show the dependence upon both 9? and % for the first and second modes 
and upon ? for the third mode. They are in excellent agreement with existing experimental data. 


7. Block designs, by Professor Burton W. Jones, University of Colorado. 


The definition and significance of balanced incomplete block designs are briefly described and 
methods of exclusion sketched. 


8. The mapping of the circles of S2 into the points of S3, by Professor C. R. 
Wylie, Jr., University of Utah. (By invitation). 

If the coefficients a, b, c, d in the general equation of a circle 

d(x? + y*) — 2ax — 2bv+c¢=0 
are interpreted as homogeneous point coordinates in S;, point circles are mapped into points on the 
paraboloid 
V cd =0, 

proper real circles are mapped into finite points outside V, improper real circles (lines) are mapped 
into real points at infinity, and imaginary circles are mapped into finite points within V. Pencils 
and bundles of circles are represented in S; by lines and planes, respectively, and may be classified 
according to the intersection of their images with V. Two circles which are orthogonal are repre- 
sented by points each of which lies in the polar of the other with respect to V. Conjugate pencils of 
circles are represented by lines conjugate with respect to V. Various theorems from college geometry 


were interpreted in S;, and the classical constructions for circles satisfying three conditions were 
considered as problems in descriptive geometry in 53. 


9. The University of Colorado Engineering Experiment Station analog com- 
puter—The UCEESAC, by Mr. Walter W. Varner, University of Colorado. 

A description of the Boeing analog computer recently installed at the University of Colorado 
was given. Types of problems that can be solved as well as restrictions on their solution were given. 


Finally a simple pair of simultaneous differential equations were considered and the simplicity of 
forming the wiring diagram shown. 


10. Fitting empirical equations to fluid meter data, by Professor S. R. Smith, 
University of Wyoming. 
Empirical equations of the form 
R 
Gp 
ad a+bR 


were fitted to both flow nozzle and orifice meter data and residuals determined. R, C curves were 
fitted to data for the fluids stream, oil and water for 0<R 33,040,000. C is the coefficient of dis- 
charge of the meter and R its corresponding Reynolds number, both dimensionless. 


11. On f(x) = F(x), F(x) given (real), f(x) unknown, by Professor (Emeritus) 
A. J. Kempner, Dr. W. E. Briggs, Professor S. Chowla and Mr. W. E. Mientka, 
University of Colorado, presented by Professor Kempner. 


For the real case the geometrical interpretation employed in studying f(x) —x =0 can be ex- 
tended so as to lead immediately to results such as: there exist totally discontinuous functions f(x) 


* Equation first used by I. D. Ruggles, former graduate student in mathematics, University 
of Wyoming. 
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for which ff(x) is single-valued and continuous; or, the function y=g(x) given by x—y+7x/2 
= cos (x+y), | | 31/2, is a single-valued inverse iterate of f(x) =x-++-x. It also makes plain the 
plausibility of introducing iterations of f(x) of any rational (or even any real) index n. 


12. Figure it out for yourself, by Professor A. W. Recht, University of Den- 
ver. 


Trained mathematicians are at a premium; everyone realizes the great part mathematics plays 
in a technical civilization. Yet election of mathematics in high schools is waning, despite efforts 
of government and mathematical groups stressing urgent need of training of the talented. The prob- 
lem is to reach the 85% who will use mathematics in normal pursuits and also train the talented 
15%. The paper suggested revisions of policy, including training of really professional mathe- 
matics teachers, emphasis on students doing own work, insistence on 100% accuracy, promotion of 
“first the problem, then the mathematics,” and strengthening of fundamental concepts to avoid 
mathematical accidents in home and industry. 


13. Textbooks on elementary mathematics, by Professor J. S. Leech, Colorado 
College. 


Attention is called to the fact that in a large majority of textbooks many terms are defined 
erroneously or ambiguously. Many theorems are stated and “proved” without any or with in- 
complete hypotheses. In all of these cases, the author believed that the correct definitions and 
statements of theorems not only do not increase the difficulty of the subjects, but serve to clarify 
concepts and contribute greatly to understanding. 


14. Freshman mathematics separation, by Professor I. L. Hebel, Colorado 
School of Mines. 


This paper dealt with the mathematics department’s experience over the last seven years in 
assigning and classifying freshman mathematics students. 


F. M. CARPENTER, Secretary 


THE APRIL MEETING OF THE SOUTHWESTERN SECTION 


The fourteenth annual meeting of the Southwestern Section of the Mathe- 
matical Association of America was held at Arizona State College, Tempe, Ari- 
zona, on April 16, 1954. Professor M. S. Hendrickson, Chairman of the Section, 
presided. 

Forty persons attended the meeting including the following twenty-seven 
members of the Association: 


O. B. Ader, C. E. Aull, C. E. Buell, J. H. Butchart, D. G. Duncan, J. F. Foster, Jr., R. S. 
Fouch, F. C. Gentry, R. F. Graesser, R. E. Graves, W. P. Heinzman, M. S. Hendrickson, Carol 
Karp, Max Kramer, Lincoln LaPaz, R. B. Lyon, W. W. Mitchell, Jr., E. D. Nering, J. L. Olpin, 
E. J. Purcell, L. C. Snively, A. H. Steinbrenner, Deonisie Trifan, Earl Walden, D. L. Webb, Charles 
Wexler, Oswald Wyler. 


The following were elected officers for the year 1954: Chairman, Professor 
D. L. Webb, University of Arizona; Vice-Chairman, Professor R. L. Westhafer, 
New Mexico College of Agriculture and Mechanic Arts; Secretary-Treasurer, 
Professor W. W. Mitchell, Jr., Phoenix College (four year term); Lecturers, 
Professor Max Kramer, New Mexico College of Agriculture and Mechanic Arts 
(one year), Professor D. G. Duncan, University of Arizona, (two years). 
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An invited address “The Story of Palomar” was given following the banquet 
by Dr. S. B. Nicholson of the Mount Wilson and Palomar Observatories. Dr. 
Nicholson also read a paper on “Results from Palomar” during the morning 
session. 

The following papers were presented: 

1. Some problems in electrostatics treated synthetically, by Professor J. H. 
Butchart, Arizona State College. 


The author showed that synthetic methods, notably inversion, are appropriate in problems 
concerning electrostatic forces in the plane. He discussed in particular the case of concyclic charges. 


2. Nonlinear differential equations that can be reduced to linear, by Professor 
M. S. Hendrickson, University of New Mexico. (By title.) 


The most general nonlinear equation of the first order which is reducible by a substitution of 
the form u=¢(x, y) to the linear equation u’=up(x)+q(x) is of the form »’ =—(¢,/¢,) 
+(1/¢,)[p@+q]. In particular, any equation of the form »’ =h(y)[F(x)+G(x) exp —/1/hdy] is 
reducible by the substitution u =exp [fi /hdy—Jf. to the equation u’ =G(x) exp —/Fdx which 
can be solved by integration alone. 

The most general second order equation reducible by a substitution of the form F(y’) =G(y) to 
a first order equation linear in G with yas independent variable is of the form y” +p(y) F(y’)y’/F’(y’) 
+q(y)y’/F’(y’) =0. In particular, any equation of the form y’’+g(y)(y’)*+h(y)(y’)?* =0 can be 
reduced by the substitution G(y) =(y’)*. 


3. Set representation theorems in implicative models, by Mrs. Carol R. Karp, 
New Mexico College of Agriculture and Mechanic Arts. 


Implicative Models were introduced by L. Henkin in his paper Am algebraic characterization 
of quantifiers, Fundamenta Mathematicae, 1950. They are the algebraic models for the System of 
Basic Implication, a calculus of propositions with only implication as a propositional connective, 
yielding as formal theorems only wffs valid in both the classical and intuitionistic calculi of proposi- 
tions. 

The result announced was that every Implicative Model ¥ = < X,0, +> can be extended to a 
Brouwerian Algebra with preservation of infs and exactly those tops which satisfy the following 
conditions: (i) s=top Y in the sense of Henkin, 7.e., (1) z is an upper bound of Y; (2) for any 2’, a 
in X, y+aSz' for all y in Y implies s+a $2’. (ii) for any 2’, a, a2, ° ++, @, in X, m any finite in- 
teger, + for all y in Y implies z+a,;+ - - +a, 

By a result of Tarski and McKinsey [Closed elements in closure algebras, Annals of Mathemat- 
ics, 1946] it follows that every Implicative Model is isomorphic to a subspace of closed sets of a 
topological space. 


4. Axioms of congruence for absolute geometry in a bounded domain, by Pro- 
fessor Oswald Wyler, University of New Mexico. 


A system of axioms of congruence is given which is valid in any open convex domain (whether 
bounded or not) of a euclidean, hyperbolic or elliptic geometry. Together with axioms of incidence 
and of order given elsewhere by the author (Duke Math. Jour. vol. 20, 1953; Composito Math., vol. 
11, 1953), and with an axiom of continuity, these axioms permit the construction of the full ab- 
solute geometry from the geometry of any open convex domain. No axiom of parallelism is used, 
but the type of a geometry (euclidean, hyperbolic, or elliptic) is given by properties of its absolute 
polarity. 
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5. A method of approximating moments in thin elastic plates, by Professor 
Deonisie Trifan, University of Arizona. 


A method of approximating the bending and twisting moments in thin elastic was presented. 
Based on a theory developed by Prager and Synge who regarded a state of stress in a deformed 
body as a point in function space, the actual solution which satisfies the equation of equilibrium, 
compatibility conditions, and boundary conditions is approximated by a series of so-called artificial 
states which satisfy only some of the above conditions. 


6. The case for intuitionism, by Professor E. J. Purcell, University of Ari- 
zona. 


7. A three dimensional model of a lattice, by Professor D. L. Webb, University 
of Arizona. 


The display of a symmetric model in three dimensions of a free modular lattice of three gen- 
erators. 


8. On normal subgroups of even order, by Professor K. A. Fowler, University 
of Arizona, introduced by the secretary. 


Let G be a finite group of even order m which contains exactly m elements of order 2. The 
following results are obtained: (1) h—12(m*+-m)/n where h is the number of conjugate classes 
each of which is equal to its inverse class; (2) there exists a real character, distinct from the 1- 
character, of degree less than m/n; (3) G properly contains a normal subgroup of bounded index, 
the bound depending only on the ratio m/n. 


9. Families of lines, by Dr. Andrew Sobczyk, Los Alamos Scientific Labo- 
ratory. 


The family of all lines in the plane is in one-to-one correspondence with the Mdébius cylinder, 
since, as may be seen from the normal form, each line is characterized by two coordinates p, a 
with —« <p<o, 0Sa<z and the boundary a=z of the p, @ strip must be identified into the 
boundary a=0 with reversal of the sign of ». The family of all lines in space is in one-to-one cor- 
respondence with the fiber-bundle of the tangent planes to a sphere. Many new and interesting 
results on line families are obtained by study of their graphs in these fiber bundles. 


10. Coalition formation in n-person games, by Professor E. D. Nering, Uni- 
versity of Minnesota and Goodyear Aircraft Company. 


The von Neumann and Morganstern theory of m-person games is essentially an extrapolation 
of an intuitive analysis of the zero-sum 3-person game. Except for this case the theory gives little 
information about the process of formation of coalitions nor, on the other hand, is the effect of such 
a process on the outcome of the game analyzed. The author presents an intuitive analysis of the 
zero-sum 4-person game which is an extension of the analysis of the zero-sum 3-person game and 
which leads to results different from those given by the von Neumann and Morganstern theory. 


11. Mathematics grades of engineers, by Professor F. C. Gentry, University of 
New Mexico. 


A survey of mathematics grades of a selected group of engineering graduates compared to 
their grades in other subjects was discussed. It was found that on the whole grades in mathematics 
are slightly lower than overall grade averages, that very few of these students failed a first year 
course, but about 25% failed one course in calculus. 


R. L. WESTHAFER, Secretary 
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EMPLOYMENT OPPORTUNITIES 


Mathematics Research, Inc., State College, Pa. Faculty consultant on industrial 
machinability research project, mathematical analysis of drilling, etc. 


The MonTSLy is devoting this space to paid announcements of employment 
opportunities for mathematicians. The text of such announcements should be in 
want-ad form and must be in the hands of the editor (C. B. Allendoerfer, Mathe- 
matics Department, University of Washington, Seattle 5, Wash.) before the first 
day of the month preceding the issue in which the notice is to appear. Announce- 
ments should indicate the academic rank or similar description of the opening, but 
should not mention a specific salary. Blind ads are permissible which direct replies to 
a specific box number in care of the Mathematical Association of America, Buffalo 
14, N. Y. In order to conserve space and achieve uniformity, the privilege is reserved 
to rearrange advertisements. Advertisers will be billed at the rate of $1.50 per line. 
Rates for display advertising may be obtained from the Advertising Manager. 


CALENDAR OF FUTURE MEETINGS 
Thirty-eighth Annual Meeting, University of Pittsburgh, Pittsburgh, Penn- 


sylvania, December 30, 1954. 


Thirty-sixth Summer Meeting, University of Michigan, Ann Arbor, Michi- 


gan, August 29-30, 1955. 


The following is a list of the Sections of the Association with dates of future 
meetings so far as they have been reported to the Associate Secretary. 


ALLEGHENY Mountarn, Pittsburgh, Pennsyl- 
vania, May, 1955. 

ILLtNoIs, Monmouth College, Monmouth, May 
13-14, 1955. 


InpDIANA, Butler University, Indianapolis, May, 


1955. 

Iowa, St. Ambrose College, Davenport, April 
15-16, 1955. 

KANSAS 

KENTUCKY 

LovuIsIANA- MISSISSIPPI 

MARYLAND-DIsTRICT OF COLUMBIA-VIRGINIA, 
Georgetown University, Washington, 
D. C., December 4, 1954. 

METROPOLITAN NEw York, Queens College, 
Flushing, New York, April 30, 1955. 

MIcHIGAN, Michigan State College, East Lan- 
sing, Spring, 1955. 

Minnesota, University of Manitoba, Winni- 
peg, October 16, 1954. 

Missour!, University of Kansas City, Spring, 
1955. 

NEBRASKA 


NORTHERN CALIFORNIA, University of Cali- 
fornia, Berkeley, January 15, 1955. 

OHIO 

OKLAHoMA, Oklahoma City University, Octo- 
ber 29, 1954. 

Paciric NORTHWEST 

PHILADELPHIA, Princeton University, Prince- 
ton, New Jersey, November 27, 1954. 

Rocky University of Wyoming, 
Laramie, Spring, 1955. 
SOUTHEASTERN, Tennessee Polytechnic Insti- 
tute, Cookeville, March 11-12, 1955. 
SouTHERN CALIFORNIA, Santa Monica City 
College, March 12, 1955. 

SOUTHWESTERN, University of New Mexico, 
Albuquerque, Spring, 1955. 

Texas, Abilene Christian College, Abilene, 
April, 1955. 

Upper NEw York University of Buf- 
falo, May 14, 1955. 

Wisconsin, Cardinal Stritch College, Milwau- 
kee, May, 1955. 
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opportunities 
for 
MATHEMATICIANS 


These are excellent positions 
requiring top level work as 
consultants and educators in the 
application of IBM electronic 
computing equipment to the most 
challenging problems of business 
and science. 

The opportunities for growth in 
this vast new field are exciting. 


with 
| 5 M Employment assignments are 
available in major cities throughout 
the United States. 


Good salaries, excellent working 
conditions, exceptional employee 
benefits. 


Representatives 


Minimum requirements: Major 

or graduate degree in Mathematics, 
Physics, or Engineering with 
Applied Mathematics equivalent. 
Desirable, but not required: 

Previous experience in teaching 
applied mathematics and use of 
automatic computing equipment. 
Write, giving full details, 

including experience and education, 
to 


Dr. C. C. Hurd, Director 
Applied Science Division 
Room 102 

INTERNATIONAL BUSINESS MACHINES 
590 Madison Avenue 
New York 22, N. Y. 
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Up-to-date mathematics texts from Ginn 


ADVANCED CALCULUS (in press) 
‘by Angus E. Taylor, 
University of California at Los Angeles 
ANALYTIC GEOMETRY AND CALCULUS 
by William R. Longley, Percey F. Smith, 
and Wallace A. Wilson 
COLLEGE ALGEBRA AND PLANE TRIGONOMETRY 


by James H. Zant, 
Oklahoma A. and M. College 


Scics Offices: New York Chicago 16 Atlanta 3 
COMPANY Dallas | Columbus 16 San Francisco 34 ~—- Toronto 7 


TEXTBOOKS OF 
DISTINCTION 


Home Office: 
Boston 


MATHEMATICS AND PLAUSIBLE 
REASONING 


By GEORGE PO6LYA 


A two-volume guide to the practical art of plausible reasoning 
by the famous author of How to Solve It. Professor Pélya uses mathe- 
matics to show how hunches and guesses play an important part in 
even the most rigorously deductive sciences. 


Volume I. Inductions and Analogy in Mathematics 336 pages. $5.50 
Volume II. Patterns of Plausible Inference 240 pages. $4.50 
The set $9. 


Order from your bookstore, or 


PRINCETON UNIVERSITY PRESS, Princeton, New Jersey 
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Special Sale 
OF VOLUMES OF 


THE AMERICAN MATHEMATICAL 
MONTHLY 


Complete volumes from volume 23 (1916) to volume 52 (1945) 


inclusive, are now available at the following rates: 


Single volumes: $10 per volume 


Five or more volumes (any years) : $ 5 per volume 


We pay transportation charges if payment accompanies order. 
* 


Harry M. Gehman, Secretary-Treasurer 
Mathematical Association of America 
University of Buffalo 
Buffalo 14, New York 
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Macmillan 


DIFFERENTIAL AND INTEGRAL CALCULUS 


fifth edition 
by Clyde E. Love and Earl D. Rainville 


—An important revision of many of the old topics with the follow- 
ing new topics added: 


(1) Work (2) Circle of curvature (8) Integral test for infinite series 
(4) Summation of power series (5) Oblique and curvilinear asymp- 
totes (6) Evaluation of iterated integrals by inversion of order (7) 
Evaluation of iterated integrals by change of coordinate system 
(8) Short appendix on rigorous presentation of limits, 


1954 526 pages $5.75 


TRIGONOMETRY, sohn 


Using the “general angle” approach to the subject, Dr. Randolph 
has written a trigonometry text adaptable to courses with emphasis 
ranging from the simplest numerical work to modem stress on 
analytical trigonometry. He presents the subject in such a way that 
the student not only learns facts about trigonometry but also ob- 
tains experience in sound thought processes. 


1953 220 pages $3.00 


The Macmillan 


60 FIFTH AVENUE, NEW YORK 33 
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PLANE TRIGONOMETRY 


By Gorpon Futter, Alabama Polytechnic Institute. 270 pages, (with tables), 
$3.25: (text alone), $2.75 


Here is a work designed as a standard text for courses in all college mathematics 
departments. The author covers both analytical and numerical trigonometry and 
avoids explanations that are too brief. Each new topic presented is illustrated with 
problems worked out in detail. 


GENERAL COLLEGE MATHEMATICS 


By W. L. Ayres, C. G. Fry, and H. F. S. Jonan, Purdue University. 288 
pages, $3.75 


A new approach to freshman mathematics, designed to help the student in his ap- 
plication of mathematical topics to the biological and social sciences. Keyed pri- 
marily for those students who will take only one year of college mathematics, this 
text examines such subjects as ratios and percentages, linear and quadratic equa- 
tions, trigonometry, interest and its application to installment buying, laws of 
growth, statistics, etc. 


CALCULUS AND ANALYTIC GEOMETRY 
By C. T. Hotmes, Bowdoin College. 416 pages, $5.00 


Designed for a «ombination course in which the concepts and techniques of the 
calculus are the aain objectives. Although calculus is emphasized, the essentials 
of analytic geometry are presented in sufficient detail for a subsequent major. The 
concept of integration is introduced early in the text. 


DIFFERENTIAL EQUATIONS 
By Rosert C. Yates, United States Military Academy. 215 pages, $3.75 


This text is designed to prepare the student for work in modern engineering prac- 
tice and theory and to present the basic mathematical tools necessary for analysis 
and solution of problems leading to differential equations. The book ranges from 
first order equations through general linear ordinary equations with constant co- 
efficients, special equations mostly of a nonlinear character, numerical solutions 
and solutions by series, and then into the solution of the wave equation by separa- 
tion of variables. 


Send for copies on approval 


McGRAW-HILL BOOK COMPANY 
330 West 42nd Street ° New York 36, N. Y. 


GEORGE BANTA PUBLISHING COMPANY, MENASHA, WISCONSIN 
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INTRODUCTION AND CRITIQUE 


F. J. WEYL, Investigator, NRC Committee on Training and Research 
in Applied Mathematics 


1. Sponsorship and Organization. The papers appearing here under a single 
cover were presented in the course of a somewhat novel type of symposium on 
special topics in applied mathematics which was organized by the Mathematics 
Division of the National Research Council and accommodated by the American 
Mathematical Society as an integral part of its regular meeting at Northwestern 
University, Evanston, Illinois, on November 27—28, 1953. The character of this 
symposium was largely set by the fact that its organization constituted one of 
the phases of a Survey of Training and Research in Applied Mathematics which 
the National Research Council has been conducting during the past year under 
provisions made by the National Science Foundation jointly with the three 
military Offices of Research in response to a recognized need for a comprehensive 
stock-taking of this country’s resources in the way of mathematics and mathe- 
maticians likely to play a significant part in the future advances of our science 
and technology. 

The arrangement of one or several conferences on the various phases of train- 
ing and research in applied mathematics formed part of the plans for the Survey 
from the beginning. As subsequent work brought the picture better into focus, 
it became clear that the modern range and depth of the subject matter of applied 
mathematics was as yet insufficiently understood by many mathematicians. To- 
day applied mathematics is emphatically no longer just the art of solving by 
dull methods repulsive differential equations whose coefficients carry physical or 
engineering names! Rather than merely pointing out the necessity of a systematic 
effort to dispel this deadly myth, the NRC Committee on Training and Research 
in Applied Mathematics which was in charge of the Survey decided to take at 
once an active step in this direction. Hence one of the planned conferences was 
set aside for this particular purpose and developed into the Symposium on Spe- 
cial Topics in Applied Mathematics here reported. 


2. Aim and Impact. The Symposium was designed to present characteristic 
current research to a large audience of mathematicians, illustrating particularly 
active sectors of the front along which mathematics interacts today with other 
scientific disciplines. Its philosophy differed therefore markedly from the Ap- 
plied Mathematics Symposia convened annually by the Mathematical Society 
or the younger series held under Office of Ordnance Research sponsorship. No 
attempt was made, as in the latter, to aim the program at the workers and 
specialists in some unified area, however broadly defined, and the interchange of 
information among them. Instead, the emphasis was on diversity of topics, on 
highlighting some exciting developments of the day and doing it in a manner to 
induce broad understanding among the mathematically literate. What we 
wanted was a collection of significant examples, not brought as yet too often or 
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accessibly to the attention of the mathematicians, which would show just how 
varied and substantial the contributions of the genuinely secular mathematician 
(as distinct from the monastic one) can be in this interaction of his modes of 
thought and analysis with scientific endeavors in other fields. 

The Symposium was therefore planned in three 90-minute sessions of three 
half-hour papers each. The first of these sessions was to be concerned with new 
mathematical theories and techniques which are currently being fashioned for 
applied purposes; the second one was to be dedicated to the characteristic mix- 
ture of combinatorial analysis and probability theory now being developed for 
the sake of constructing analytic models and the numerical simulation of con- 
trol and communication systems; the third one, finally, was to deal with situa- 
tions where the interplay between physical ideas and mathematical construction 
is at present the critical element. This structure of the sessions is still discernible 
in the Symposium as it took place in the end, although an unusually large num- 
ber of last minute changes beyond the planners’ control compromised it to some 
extent. It proved helpful, however, in assuring the diversity and range which was 
in fact achieved. 

The theory of partial differential equations, and especially of the physically 
relevant non-linear types, shows no signs of exhaustion as a primary source of 
important and stimulating problems which fascinate the mathematicians of 
applied interests. There are the basic questions revolving around the idea of 
what constitutes a well-formulated problem in this area which J. Leray touches 
in his paper and which have naturally led to the viewpoint, pioneered by Leray 
and others, of looking at solutions not individually but as families with topologi- 
cal and algebraic structure, capable of being supplemented by ideal elements, 
etc. There are the powerful function theoretical arguments which are being de- 
veloped for the construction and analysis of important particular solutions, here 
represented by the methods with which P. Garabedian and his co-workers have 
recently attacked hydrodynamic flows possessing free surfaces. There are finally 
the fascinating problems posed by the nature of the singularities which can 
arise in the solution of these non-linear systems and the structure of the loci on 
which they can occur: A. H. Taub discusses such a problem using methods of 
differential geometry to analyze the configuration of interacting shock fronts in 
compressible flow. All of these are vital areas of substantial mathematical re- 
search today. 

Ultimately destined, perhaps, to be of comparable importance in its mutually 
stimulating interaction with modern analysis is the area of signal and noise prob- 
lems, illustrated by those presented in the paper of M. Kac. Here again the quest 
for analytic and measure-theoretical information on certain families of functions, 
containing individually observed time series as members, leads naturally to 
problems of substantial interest also to current front-line thinking in pure analy- 
sis. Another fertile field for mathematical ingenuity in applied contexts is furn- 
ished by what, in effect, are elementary problems but which are posed with 
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such a massive complexity of detail as to render the solution totally inaccessible 
to the traditional class-room approach. S. Chandrasekhar’s discussion of eigen- 
value problems of high order furnishes a good example how such an obstacle 
can be overcome by an ingenious mathematical trick, perhaps even, as in this 
case, without knowing in a rigorous way what accounts for the success of the 
trick. 

Probably the least familiar aspect of the “new look” in applied mathematics 
is the steadily growing interaction between the combinatorial and algebraic 
mathematics of discrete structures on the one hand, and the use and operation 
of large-scale digital computers on the other. The use of Boolean algebra in the 
analysis of switching circuits is of course no longer new, but the idea of carrying 
forward basic research in Boolean algebras, as reported in the sequel by D. E. 
Muller, with the hope of anticipating the future needs of switching circuit design 
rather than being merely concerned with supplying the present ones, and to 
make use, moreover, of high-speed computing equipment in doing it—this 
gives evidence of the substantial development of a major new province of 
applied mathematics. In the same general area lies the work of Project SCAMP 
at the Institute for Numerical Analysis, whose purpose is that of exploring 
the use of high-speed automatic computers in the analysis of discrete structures. 
Some of it has been reported for the first time in any completeness at the Sym- 
posium by S. S. Cairns. Since summer after summer a sizable group of our ablest 
mathematicians participate in this work, it cannot fail in the long run to have a 
noticeable impact on American mathematical research. 

Another important development, once more in the traditional vein of applied 
mathematics, is illustrated by the investigations, summarized by E. Montroll, 
concerning the energy distribution over the frequencies at which a molecular 
crystal lattice is capable of oscillating. It turns out that this analysis depends in 
an essential way on the relations known to exist between the topological struc- 
ture of a closed manifold, here a certain kind of phase space, and the analytic 
properties of functions and integrals defined thereon. These problems are among 
the most intensively investigated ones in modern analysis, and if they should 
furnish the appropriate terms for analyzing significant aspects of current theo- 
ries of molecular physics, we might be faced with their proliferous growth similar 
to the development of Hilbert space methods under the impact of quantum 
mechanics. 

Everyone who scans this list will miss certain topics which he would consider 
as particularly interesting or important for the development of modern applied 
mathematics. Let him be assured that but few, if any, of these omissions were 
inadvertent. It is perhaps indicative of how much needs doing along these lines, 
that a first attempt could achieve as much and still leave so much untouched. 
Probably most glaring is the absence of a report, aimed at mathematicians, on 
the current problems and difficulties faced by quantum mechanical theories of 
fields. All endeavors of the Committee to provide such a report were unsuccess- 
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ful, and in hindsight this may have been all to the good, for the matter is serious 
enough to warrant a special effort. 


3. Experimental Features. In arranging for a novel type of meeting, some 
innovations as regards organization are inevitable. The Symposium on Special 
Topics in Applied Mathematics introduced two such departures from standard 
practice. For one, the sessions of the Symposium were interspersed among the 
regular sessions of the Mathematical Society meeting, rather than grouped into 
a separate conference preceding or following the former and thus too easily 
avoided by the non-specialist. The example which was followed here is that of the 
American Physical Society whose meetings provide for a sizable number of in- 
vited half-hour papers, mixed in among the contributed papers either singly or in 
groups. Thus, the speakers were permitted some expository breadth in order to 
bring those of their colleagues not immediately concerned with the subject on 
hand up to date. The second innovation consisted in inviting representatives of 
other scientific fields to speak to mathematicians. The ones selected for this pur- 
pose naturally have extensive contact with mathematics in their daily scientific 
lives, yet they do not think of themselves as mathematicians. It is from such 
persons that a mathematician will generally get his most stimulating insight 
into the mathematical problems which arise in other scientific disciplines. Both of 
these experiments proved successful at the Symposium under discussion and 
their adoption, also in other than applied mathematical contexts, is recom- 
mended. 
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THE PHYSICAL FACTS AND THE DIFFERENTIAL EQUATIONS 
JEAN LERAY, College de France 


1. Introduction. A differential equation raises not only the problem of prov- 
ing: 
(A) the existence and regularity of its solutions; their uniqueness and their 
continuity with respect to the datum or at least their compactness for the datum 
running on a compact set. 

Indeed, when those properties cannot be established, then the following 
alternatives arise: 

(B) allowing the solution to have convenient singularities, to obtain exist- 
ence- and compactness-theorems without the regularity-theorem (and indeed with- 
out the uniqueness-theorem even when uniqueness should be expected), 

(C) to construct counter-examples proving that the compactness-theorem does 
not hold. 

Assume that the equation is the conclusion of a physical theory expecting 
the existence, the regularity (and maybe the uniqueness) of its solutions. If (C) 
is possible, then the properties of the equation contradict the theory. On the 
contrary, if (A) holds, then the physical phenomenon is actually ruled by the 
equation. If not, and if (B) is possible, that is if the expected regularity (and 
uniqueness) is not obtained, then the physical theory neglected something 
which appears in the end to be not entirely negligible. The theory does not com- 
pletely describe the phenomenon, which eludes perhaps any complete descrip- 
tion. This often happens when the equation is non-linear, especially in the me- 
chanics of fluids; their equations are based on the assumption that the flow is 
laminar, but the experiment shows turbulence. 

Many papers have been devoted to problem (A); few to (B) and (C). The 
processes required for the establishment of (A) are often difficult; wherever (B) 
and (C) have been solved, it was by means of very simple processes, although at 
times somewhat concealed by technical details. We shall describe some of these 
here. 


2. An equation which rules an evolution and for which (B) holds. Consider 
the Navier-Stokes system ruling incompressible viscous flows; it is parabolic 
and non-linear; the velocity is given at the initial time ¢=0 and has to be ob- 
tained for any t>0. 

By successive approximations the flow can be calculated from t=0 to t= 7), 
then from t= 7, to t=T>, etc. Let T=lim T,,; the problem is solved and its solu- 
tion is regular if T=-+ ©. That occurs in some special cases: two-dimensional 
flows without boundary, sufficiently small data, etc. The proof is based on well 
known properties of the vortex; it can also be based on the relation expressing 
the conservation of the energy. 
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But the homogeneity of the Navier-Stokes system is such that the conserva- 
tion of the energy and the smoothing action of the viscosity seem to be insuffi- 
cient for maintaining the flow so regular that T=+ 0. 

However, the experiment shows that the flow exists from t=0 to t=+o, 
Furthermore the Navier-Stokes system can be slightly modified in such a way 
that: 

(a) the energy relation is preserved; 

(b) the smoothing properties of the viscosity are so increased that the flow 
cannot become irregular and hence exists from t=0 tot=+. 

Let this modified system tend to the Navier-Stokes system; its solution con- 
verges weakly to one (or several) solutions of the Navier-Stokes system. This 
solution is not regular; it is only square integrable at any moment of time; it 
satisfies the Navier-Stokes system in the sense of L. Schwartz’s distribution 
theory. Problem (B) is solved. 

In some cases (three-dimensional flow without boundary, two-dimensional 
flow inside a convex boundary) a detailed study shows that such an irregular 
solution is regular except for a compact subset of values of ¢ of measure zero. 

The uniqueness of the flow and the continuity of its kinetic energy at the 
moments of time where it is not regular are two open problems; they are related. 

No example of irregular flow has been obtained, but the above-mentioned 
homogeneity of the Navier-Stokes system suggests a process which probably 
leads to such an example. 


3. An equation ruling an equilibrium and for which (A) holds in some cases, 
(C) in the other cases. Consider the Dirichlet problem for non-linear, second 
order, elliptic differential equations in two independent variables. The fixed point 
theory and the theory of the linear elliptic equations reduce the existence theo- 
rem to the compactness theorem, more precisely to the obtaining of a@ priori 
bounds for the solution and its first and second derivatives. These bounds fol- 
low from an extension of the maximum principle to the solutions of the studied 
equation, to their first and to their second derivatives. Assumptions about the 
equation and the Dirichlet data have to be made; they mean that the equation 
has a convenient behavior on the vertical cylinders (especially on the cylinder 
containing the graph of Dirichlet’s data) or, if the equation belongs to Monge- 
Ampére’s type, on the curves (especially on this graph). 

The proof that those sufficient assumptions are necessary (when the equa- 
tion behaves regularly on those cylinders and curves) is quite simple. It consists 
in the construction of a solution of the equation such that: 

(a) this solution has at a point a vertical tangent plane or an infinite curva- 
ture; 

(b) this solution is the limit of regular solutions defined on the same domain. 

This is easily done by applying the Cauchy-Kowalewski theorem to the equa- 
tion transformed by a change of coordinates or by a contact transformation. 
The non-linear Dirichlet problem in three or more independent variables is 
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open. A priori bounds cannot be obtained again by a maximum principle; no 
existence or compactness theorems are known; problem (A) is not solved and 
probably problem (B) should be considered. However, the construction of the 
counter-examples showing that (C) holds could easily be extended. 
Bibliography 
Navier-Stokes system 
Its origin, the linear Navier-Stokes system: 
C. Oseen, Hydromechanik, (Leipzig, 1927). 


Three-dimensional flow without boundary: 
J. Leray, Acta math., vol. 63, 1934, pp. 193-248. 


Two-dimensional flow inside a convex boundary: 
J. Leray, Journal de math., vol. 13, 1934, pp. 331-418. 
Regular flows: 


J. Leray, Journal de math., vol. 12, 1933, ch. iv, pp. 64-82. 

J. Leray, Comptes rendus Acad., vol. 194, 1932, pp. 1892. 

W. Wolibner, Math. Zeits., vol. 37, 1933, pp. 698-726. 

E. Hélder, Math. Zeits., vol. 37, 1933, pp. 727-738. 

J. Leray et L. Robin, Comptes rendus Acad., vol. 205, 1937, p. 18. 
P. Casal, Congrés intern. Mécan., 1952. 


The general case: 
E. Hopf, Math. Nachrichten, vol. 4, 1950, pp. 213-231. 
Non-linear Dirichlet's problem 
First studied by: 
S. Bernstein, Math. Ann., vol. 69, 1910, pp. 82-136; vol. 95, 1926, pp. 585-594; vol. 96, 1927, 
pp. 633-647; Annales Ecole norm., vol. 27, 1910, pp. 233-256; vol. 29, 192, pp. 431-485. 
His preliminaries about linear equations have been largely simplified by 
J. Schauder, Math. Zeits., vol. 37, 1933, pp. 623-634; Studia Math., vol. 5, 1935, pp. 34-42. 
L. Bers, L. Nirenberg (to be published) 
The uniqueness assumption has been removed by 
J. Leray, J. Schauder, Annales Ecole Norm. Sup., vol. 51, 1934. 
H. Lewy and J. Schauder found a mistake in S. Bernstein’s a priori bound of the first derivative 
of a quasi-linear equation (Math. Ann., vol. 69 p. 129: the expression of K is obviously wrong); 


furthermore, his conclusions are not invariant with respect toa change of coordinates; his study of 
the Monge-Ampére equation is incomplete. His results have been corrected and completed by 


J. Leray, Journal de Math., vol. 17, 1938, pp. 89-104; vol. 18, 1939, pp. 249-284. 


y 
V 

t 
n 
| 
d 
t 
4 
d 
e 
n 
l- 
1. 
is 


RECENT ADVANCES AT STANFORD IN THE APPLICATION 
OF CONFORMAL MAPPING TO HYDRODYNAMICS 


P. R. GARABEDIAN, EDWARD McLEOD, JR., and MARTIN VITOUSEK, 
Stanford University 


The present paper forms an expository report of researches carried out re- 
cently at Stanford University by a group of mathematicians interested in free 
boundary problems of hydrodynamics. The discussion centers about a unified 
technique in geometric function theory. 

The steady, irrotational motion of an incompressible, inviscid fluid in the 
z-plane is governed by an analytic complex potential ¢(z) =¢+7. For a fluid of 
unit density, the pressure p is given by Bernoulli’s equation 


(1) 4| ¢’(2) |? + gy + p = const., 


where g is the gravitational constant. This law, in general, determines an addi- 
tional boundary condition along free streamlines. 

Following [3], we consider the uniform plane flow of a liquid past a small 
bubble of gas. We assume that the bubble is so small that the forces due to 
surface tension surpass significantly those due to gravity. Setting g=0 in (1) 
and resolving all forces along the normal to the free streamline, we obtain on 
the boundary C of the bubble a condition of the form 


(2) | ¢’(2) |? = Tx, 


where x is the curvature of the free boundary C and T represents the magnitude 
of the surface tension. 

Denoting differentiation with respect to arc length s along C by a dot, we 
throw (2) into the form 


(3) = 

and we obtain along C upon integration 

(4) = f = 

where g(z) is an analytic function in the flow region D. Hence g(z) maps D con- 
formally upon a Riemann surface bounded by a circle of radius T, while ¢ maps 


D conformally upon the exterior of a horizontal slit. If w(z) maps D upon the 
exterior of the unit circle, we can therefore write 


1 
(5) f=w+—)> 
w 


(6) q = Tw ———_ A?> 1. 
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1 2 
(1 =) dw 
7 1 f w? 
@) tT d [ w? + A? ] 
w 
dw 1+ A*w? 
where A is to be determined so that the numerator and the denominator of the 


integrand vanish simultaneously. We obtain the explicit formula 


2 1 
(8) s=w-—- 


But since dgdz =df{?, we find 


by integrating (7), and, together with (5), this determines completely the shape 
of the bubble and the corresponding flow. 

In [4] a treatment is given of a free surface flow representing a water foun- 
tain in which gravity, rather than surface tension, is predominant. Taking 
g=1/2 and choosing suitable coordinates, we can put (1) in the form 


(9) ——+—=0 


along the free streamline. Following [2], we denote by \(f) an analytic function 
in the flow region D which coincides with y on the free surface, and we substitute 
into (9) to obtain by analytic continuation the ordinary differential equation 


1 
(10) z'(z’ — 2id’) + 0, 


where the prime indicates differentiation with respect to ¢. Along a vertical 
linear streamline, this differential equation in the unknown ) is real, and we can 
expect to obtain a real solution for real initial data. 

The desired flow is required to rise along the y-axis until it reaches a free 
surface at the origin which spreads smoothly to either side and falls, joining two 
fixed vertical boundaries which descend to infinity. The continuation of the 
free streamline is linear and vertical, so that we can expect A(z) to map it onto 
the real axis. Since {(z) has the same mapping property, we are led to set 
A= in (10). Integration then gives 


(11) 
0 


a formula which does, indeed, represent a gravity flow of the type described. 
These techniques lead to further interesting free surface flows involving, 

respectively, surface tension and gravity. However, such examples do not ex- 

haust the method. A successful treatment of cavitation in axially symmetric 
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flow yields to the present technique. The essential difficulty for this case stems 
from the more complicated partial differential equation 


1 ay 1 oy 
O202* 2(2—2*) dz 2(¢ —2*) dz* 


(12) =0, s=x+iy, 2* = x — iy, 


satisfied in the meridian plane by the stream function y governing the flow, 


but it can be largely overcome by appropriate exploitation of the Riemann 
function 


(13) R(z, 2*; t, *) = 


(z — t*)(2* — 2) 


associated with this equation in the domain of two independent complex vari- 
ables z and z*. Here F[w]=F(—1/2, —1/2, 1, w) is the hypergeometric series. 
Indeed, the stream function y of an axially symmetric free surface flow for which 


the analytic curve 2 =g(z) appears as the free boundary has an explicit represen- 
tation 


where 2p is a fixed point on the curve. Analysis of the conformal mapping prop- 
erties of the analytic function g(z) analogous to the study of g(z) and A(z) carried 
through above yields explicit examples based on (14) of axially symmetric cav- 
ity flow in the large [1]. 
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SINGULARITIES ON SHOCKS 
A. H. TAUB, University of Illinois 


1. Introduction. It is our purpose to discuss the situation which arises when 
a shock front is such that there are points on it at which the tangent to the 
shock, the curvature of the shock or some derivative of the curvature is discon- 
tinuous. Such shocks will be called shocks with singularities. They occur in many 
shock interaction problems; for example, the shock MTRO, in Figure 1, which 
illustrates the shock configuration arising in the Mach reflection of a plane 
shock from a rigid wall [1], may be considered as a shock with a singularity at 
the point T at which point the tangent to this shock is discontinuous. 


Fic. 1. Incident shock is TI; reflected shock is TRO; Mach shock is TM; 
slip stream is TS; wall is CW. 


We shall restrict the discussion to two dimensional pseudo-stationary and 
two dimensional stationary flows. It will be sufficient to discuss in detail only 
one of the cases which will be chosen to be the first one, for the key equations 
needed in the discussion have been shown, in a previous publication [2], to be 
identical. 

The shock configuration in a two dimensional pseudo-stationary flow at 
time ¢ may be represented by the equations 


x = ta,(s) 

ta,(s), 

where x and y are fixed Cartesian coordinates in the plane in which the shock 
11 


(1.1) 
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is moving, and s is a parameter along the shock which may be chosen so that 


at regular points. At a singular point of the shock some derivative of the func- 
tions a,(s) is discontinuous. 

We shall assume that at time ¢ each point of the region behind the shock 
with a singularity is occupied by a single fluid particle which has crossed the 
shock at some time fo earlier than ¢. 

Let us denote by so the value of s corresponding to a singular point of the 
shock and refer to the point on the shock as the point so. The locus of particles 
which have crossed the shock at so at times fp <¢ must be a unique locus separat- 
ing the region behind the shock with a singularity into just two regions in ac- 
cordance with the assumption made above. However, this locus will be a locus 
of discontinuities in the flow behind the shock and therefore must be either a 
slip-stream or a shock. Only the former possibility is allowable as will now be 
shown. 

The position at time ¢ of a particle which has crossed a shock at a point s 
at time fo is given by 


= } 
(1.3) x = ta,(s, 7) 
y = ta,(s, 7), 
where 
t 
(1.4) rt = log — 
to 
and a,(s, 7) are solutions of the equations 
da; 
dr 


where u; is the particle velocity relative to the fixed Cartesian coordinate sys- 
tem; the a,(s, 7) satisfy 


(1.6) a,(s, 0) a,(s), 


where the functions on the right hand side of this equation are the functions 
entering into equations (1.1), the equations specifying the shock. 

Equations (1.3) or the functions a;(s, r) for fixed s and variable to give the 
locus of particles which have crossed the shock at a point s at times fp St. We 
shall call these loci streak-lines. Our fundamental assumption implies that for a 
shock with a singularity at so that the locus given by 


(1.7) lim a;(so + 7) = az (7) 
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be identical with the locus given by 


(1.8) lim a;(so — €, 7) = a; (r) 
0 


where 7 is the parameter along this locus, defined by equation (1.4) with ¢ fixed 
and to variable. The two functions obtained by these two limiting procedures 
need not be identical in r but they must represent the same curve. In Figure 1 
this curve is represented by T'S and called the slip-stream. 

From equation (1.5) it follows that the tangent vector to the locus of particles 
that have crossed the shock at s at times t9<¢ has components proportional to 
U;. We denote by Uj‘ (ai, az) the values of U; at the point a; on or above the 
singular locus described by (1.7) (in Figure 1 this region is STR) and by 
U;-(a1, a2) the corresponding quantity at the point a; on or below the singular 
locus described by (1.8) (in Figure 1 this region is STM). Since these two loci 
must be identical we must have 


(1.9) (a7 Ua (ae = 0. 


That is, their tangent vectors at the same points are equal: 


_ 


(1.10) 
where 
(1.11) (0°)? = (U3)? + (Ua)? 


Equation (1.9) states that there is no flow across the singular locus (1.7) (or 
(1.8)) and hence this locus must be a slip-stream. It is a consequence of (1.9) 
and the Rankine-Hugoniot equations given below that 


(1.12) (ai(r)) = p (ai 


where p*(a,) is the pressure at a point a; in the flow on or above the singular 
locus and p~(a,) is the pressure at a point a; in the flow on or below the singular 
locus. Equations (1.9) and (1.12) are identities in r, and af (r) represents a point 
on the singular locus. 

In the subsequent sections of this paper we shall show how equations (1.10) 
and (1.12) and the derivatives of these equations may be used to determine 
geometric properties of the shock configuration at a singular point. 

Since equations (1.9) and (1.12) are identities in r, the derivatives of these 
equations with respect to r must also hold. It follows from equation (1.9) that 


Us + = 0, 


where the comma denotes the ordinary derivative with respect to variable a 
and the argument of the functions Uj and U; is af(r). In view of equation 
(1.10), this equation may be written as 
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However, the curvature of a curve whose tangent vector is proportional to 
U; is given by 


€%2U;,;U Up 


(1.14) 


Hence it follows from (1.13) that 
(1.15) K+(at(r)) = K-(a*(r)), 


where Kt is the curvature of the locus described by (1.7) and K- is the curva- 
ture of the locus described by (1.8). Differentiating both sides of equation (1.15) 
we obtain 


(1.16) KY = K~’ 


where the prime denotes the derivative with respect to the arc length along the 
singular locus. Similarly we may show that 


(1.17) = 


where the superscript m indicates the mth derivative with respect to the arc 
length of the singular locus. 


We may also show that equation (1.12) implies that 
(1.18) = p-™, 


Equations (1.17) and (1.10) are merely the formal consequences of the 
requirement that equations (1.7) and (1.8) represent the same locus, for the locus 
is intrinsically specified in terms of the tangent vector, the curvature, and the 
derivatives of the curvature. 


2. The equations describing the flow. In order to determine the geometric 
properties of a shock at a singular point we use equations (1.17), (1.10), and 
(1.18) evaluated at the shock (i.e., at r=0) and the relations obtained by 
Thomas [3] relating the curvature and the derivatives of the curvature with 
respect to arc length of a streak-line in stationary flow with the corresponding 
quantities of the shock from which the streak-line emanates. The identical rela- 
tions may be used in pseudo-stationary flows as follows from the results of 
Taub [2]. Thus we assume analyticity of the flow in regions behind the shock 
front with a singularity and replace equations (1.9) and (1.12) by sequences of 
conditions at the singular point of the shock, namely, by equations (1.10), 
(1.17) and (1.18) evaluated at r =0. These equations are then replaced by condi- 
tions involving the quantities characterizing the shock by using results from 
consideration involving the differential equations that must be satisfied by the 
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flow variables and the Rankine-Hugoniot equations that must hold across a 
shock. 
For pseudo-stationary flows the former equations may be written as: the 
ad conservation of mass, 


(2.1) (eU;),« + 2p = 0, 


the conservation of momentum, 
1 

(2.2) + Uit — = 9, 
p 


_ and the conservation of energy, 
(2.3) U.U;U;,; + 0? + 2) 0, 
where p is the density and 


(2.4) c= =.. 
p 


The Rankine-Hugoniot equations may be written as 


(2.5) U, = Ud = = Vu 

(2.6) 
. y 
(vy + 1)o? 

2:7 = 

is (2.8) Un = = ((y — 1)o? + 2), 


(y+ De 
where \; are the components of the tangent vector to the shock, 


Urn 


id (2.9) o 
C1 C1 
th p 
1g (2 10) 7= re 
a- 
of (2.11) 
*k P1 
f ‘ re 
) and we have used the subscript 1 to denote quantities evaluated on one side of 
‘ the shock, a side which we will assume to be uniform. Thus in Figure 1 the sub- 
a script 1 may refer to the region TRI or to the region to the right of TM. 
a In the work of Thomas referred to above, equations (2.1) to (2.3) and the 


derivatives of equations (2.5) to (2.8) along the shock when evaluated at the 
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shock were considered as algebraic equations for the quantities U;,;, p,; and p,;. 
The solutions of these equations were then used to give explicit expressions for 
the derivatives of the flow variables behind a shock in terms of properties of the 
shock. Similar methods were used to obtain relations between higher derivatives 
of the quantities U;,p and p at a point behind a shock in terms of the derivatives 
of the equations characterizing the shock. 


3. Thomas’ results. We shall describe a method for obtaining Thomas’ 
results and the addition equations needed for the carrying out of the procedure 
outlined above. We consider the functions 


(3.1) a; = a,(s, 7) 
which are solutions of equations (1.5) subject to the initial conditions (1.6) as 


a transformation from the variables a; to a new set of independent variables s 
and r. Then 


(3.2) 
Or 
(3.3) =, 
Os 
where X, is defined by equations (3.3) and is such that 
It then follows that 
and 
where 
(3.7) Un = 


We now have 
OU; ar OU; as 


Or 0a; Os 0a; 
Om 


(3.8) 
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Hence 
log U, 
(3.9) Ui4 = 
Or 
and 
1 ou; 
(3.10) = 2 — U;- ) 
Os Or 
where 
(3.11) U, = 


Equations (3.9) enable us to obtain a first integral of the conservation of 
mass equation. Equation (2.1) may be written as 


op 
— + pUi + 2p = 0 
Or 


or as 


8 log p 0 log U, 
0g 
Or Or 


+2=0. 


Hence we must have 
(3.12) n = po one, 


where the arguments of the functions on the left hand side of this equation are 
s and 7 and we have introduced the convention that 


(3.13) fo = f(s, 9). 
The conservation of momentum equations may be written as 


ou; 


(3.14) 


1 
-— fs 
p 


Multiplying these equations by U; and summing we have 


1 dp 

(3.15) 
or 

(3.16) = — i 

2 Or p Or 


and multiplying by A; and summing we have 
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(3.17) au; ae 1 ap 1 dv 
ar p as 2 as 


The conservation of energy equation may be written as 
= 0 
Or Pe 


and admits the first integral : 
p 
(3.18) . 


Po. 


Hence in the s, tr coordinate system, in which the curves s =constant are the 
loci of particles which have crossed the shock at earlier times, the streak-lines, 
and the curves r =constant are a family of curves which include the shock (for 
the value r=0), the conservation equations admit two first integrals and take 
on the relatively simple form (3.16) and (3.17) in addition. The use of these 
coordinates in the integration of the equations of pseudo-stationary flows will 
be discussed elsewhere. 

Equation (2.3) may be written as 


ou; c 
U; — - — — = 2c — UW; 
Or U, Or 
If we let 
= 
then 
OU, ou; 
NU; = 0, + - 
T Or 
Hence 
and 
ou, 0U Ou; OU, 
Or Or Os U, Or 


Substituting in this equation from equation = we obtain 


aU. Us ap 1 
Os p Os 2 as 


ar U2— cr 
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or 


EGS —U.) + +t 
AY 


— 


At a point on the shock, equations (3.17) and (3.20) express the rate of 
change of the pseudo-velocity vector U along the streak-lings in terms of the 
derivatives of the flow variables along the shock. That is, the right hand sides 
of these equations may be computed from the functions a; occurring in equa- 
tions (1.1) and from equations obtained by differentiating the Rankine-Hugo- 
niot equations (2.5) to (2.8). We may also compute 0p/d7r and 0p/dr behind 
a shock from the right hand sides of equations (3.17) and (3.20), that is, in terms 
of the derivatives of the flow variables along a shock. 


(3.20) 


Thus 
Op p 0U, 
ar U, ar 
Substituting from (3.20) we obtain 
dp p [ ( (4 ap )] 
3.21 — = — —— |] U.{— — JU; }}. 
It follows from the equation preceding equation (3.18) that 
0 — pc? 0U; 
= ms Ae =) 
ar dr U2 — as 
(3.22) 


The curvature of a streak-line may also be expressed in terms of the flow 
variables and their derivatives with respect to s. From the definition of curvature 
we have 


Or 


—vK = €ik U; 


= — — 
p 


as follows from the conservation of momentum equations. However 


1 


Or 


| 
= 
+ U; ( as v as U; 
1 


20 SYMPOSIUM ON SPECIAL TOPICS IN APPLIED MATHEMATICS 
Hence 
1 /d 0 
(3.23) vK = (2 u.- 
pU,, Os 


where we may substitute from the equation (3.22) for 0p/dr. 

Equations (3.22) and (3.23) when evaluated for r=0, that is, at the shock, 
give explicit expressions for p+“ and K+ needed in the determination of the 
geometric properties of a shock with a singularity by means of equations (1.17) 
and (1.18). Explicit expressions for p+™ and K+™ may also be obtained from 
equations (3.22) and (3.23) respectively as follows: Either of these equations 
may be differentiated m times with respect to r. The right hand sides may be 
reduced to expressions involving only derivatives with respect to s. The deriva- 
tives with respect to s may be evaluated by differentiating the Rankine-Hugo- 
niot equations. 

We shall demonstrate the method outlined by evaluating equations (3.22) 
and (3.23) for r=0. It follows from equations (2.5) to (2.8) that if the flow 
ahead of the shock is uniform, then (cf. [2] eq. (4.2) to (4.6)) 


i; 
Os 
(3 24) OU 1n kU do 
OU ou; 
= kU, = —1— keyo = — 
Os Os 


where & is the curvature of the shock and \, in these equations represents the 
tangent vector to the shock. From the definition of s in equations (1.1) it follows 
that 


(3.25) (ss, 0) = AA; = 1. 
Further it may be shown that at r=0 
Op 4kU Unp 


Os y+i1 

Op 4kU on 
3.26 — = (2'= (7 = — 1) oa — 
(3.26) G+ | (y — 1)@— 1)) 

ou, 2 

ds y+1 n 
and 


ou; 
(3.27a) 1+ = — k(n —1)U, 
AY 


ay 
d 
: 
@ 
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. & n . 
( 


Substituting from these equations into (3.22) and (3.23) we obtain at r=0 


(3.28) : +~—*)ut - ut) 


Or 
and 
U; ( 8 2 
3.29 K = k— 1)U, 
respectively. 


It has been shown by Thomas [3] for stationary flows and by Taub [2] for 
pseudo-stationary flows that at r=0 


(3.30) K™ = G,(n, U./U,)k™ + H, 


where K™ is the quantity defined in section 1, k™ is the mth derivative of the 
curvature of the shock with respect to s, and H, is a polynomial in & and its 
derivatives with respect to s of order less than n. 

It follows from the arguments given above that at r=0 


(3.31) = fa(n, Ur/Un, RY, ROY), 


That is, behind the shock p™ is a function of the strength of the shock, its in- 
clination with respect to the flow incident upon it, the curvature of the shock and 
the first —1 derivatives of this curvature with respect to s. 

4. Conditions at singular points of a shock. The formulas of the preceding 
section enable us to translate the conditions which must obtain at a singular 
point of a shock, namely, equations (1.10), (1.17) and (1.18), into sequences of 
algebraic relations between the variables n+, Ut /U7, k+ and k*“ (¢=1, 2, 

A shock which has a straight portion and a singular point at which the 
curvature or some derivative of the curvature jumps to a non-vanishing value 
cannot have a uniform flow behind the straight portion of the shock at the singu- 
lar point. This follows from the fact that if the flow were uniform, equations 
(3.30) would apply for the calculation of K+™ and hence these quantities would 
vanish. However, the quantities K-™ would not vanish and it would be impos- 
sible to satisfy equations (1.17). The discussion of a special case of a shock 
with such a singularity is contained in [2], wherea shock with a point at which 
the curvature jumps from zero to a non-vanishing value is discussed. 

In a shock configuration such as shown in Figure 1, a shock with uniform 
flows incident on it and with a singularity at which the tangent is discontinuous 
may exist. The sequence of algebraic conditions which must hold may then be 
paired as follows: 


- 
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=0U; 
(4.1) 


the first of these being (1.10) and the second being (1.18) for »=0; 
K+ = K- 

(4.2) 


= 


where we have used (1.17) and (1.18); and the remaining equations may be 
written as 


Kt+™ = 


(4.3) poor) n= 1, 


Equations (4.1) are the fundamental equations of the three-shock theory 
(cf. [1]). They may be solved for n+, U} /Ux, n~ and U; /U, as functions of the 
strength and angle between the flow incident on and the normal to the shock 
TI of Figure 1. The solutions are not unique and only very special ones exist for 
some values of these two parameters characterizing the shock TY. 

Equations (4.1) and (4.2) have been solved by Clutterham and Taub [4]. 
Because the function Ho of (3.30) vanishes and because p“ depends linearly on 
k (cf. (3.28)) these equations are independent of k+ and imply a relation between 
the two variables characterizing the shock TJ. This in turn determines possible 
world lines of the point T in Figure 1. However, if either k+ or k~ is given, then 
the remaining curvatures, k~ or k+ and K+ and K-, are determined. 

The sequence of equations (4.3) may be used to determine k+) and hence 
K+ as functions of kt or k-. 

The fact that k~ cannot be determined from the equations (4.1), (4.2) and 
(4.3) is to be expected; for in a shock configuration such as in Figure 1, the 
shock TM must satisfy an additional condition, namely, it must be orthogonal 
to the wall CW. Such a requirement serves to determine k- (cf. [4]). 

Thus a shock configuration such as that occurring in Figure 1 may be deter- 
mined to an arbitrary degree of accuracy by solving equations (4.1), (4.2) and 
(4.3) with in the latter equations chosen sufficiently great and satisfying the 
boundary condition on CW. However, the solution will not be unique for a given 
shock TY. 
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SIGNAL AND NOISE PROBLEMS 
M. KAC, Cornell University 


In this brief report I shall limit myself to a cursory discussion of several 
problems chosen on grounds of personal predilection and with the view of show- 
ing how “pure” and “applied” mathematics can intermingle with profit to both. 


1. The problem of detection. This is a purely statistical problem which can 
be formulated as follows: Let x(t) be a stationary random process (noise) and 
s(t) a signal. An observer receives a record y(t) which is either x(t) or x(t) +5(t) 
and is to decide whether the signal is present or not. 

In the simplest case, when s(t) is periodic with period @ and a finite number 
of observations is made at times 4, 4: +0, - - - , 4+(n—1)@ and if furthermore 6 
is long compared with the correlation time of the noise (i.e., x(t: +70), 7=1, 2, 

++, m—1, are independent), the problem can be solved by the use of the 
Neyman-Pearson theory. A detailed account can be found in [1]. More recently 
various authors [2] considered the sequential approach to the problem. In its 
general form the problem clearly belongs to decision theory and it may be 
hoped that here the theories of the late A. Wald will play an important part. 


2. Spectra and correlation. The relation between the power spectrum 
A(w) and the correlation function p(t), 
(2.1) p(t) -f A(w) cos wt dw, 
is a standard tool in the theory of noise. It is often desirable to consider proc- 
esses for which . 
(2.2) p(t) = 0 |t| > A. 


The relations (2.1) and (2.2) combined with the fact that A(w) 20 induce strong 
restrictions on p(t). One can show, for instance, [3] that if m is an integer 


the constant being the best possible. 


(2.3) 


3. Integral equations. In the theory of radio receivers with square-law de- 
tectors one is led to the problem of finding the distribution function of an ex- 
pression 


(3.1) — 1) + y(t — 1) } dt, 


where x(#) and y(t) are independent, stationary Gaussian processes with the 
same correlation function p(t). 
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The characteristic function of the distribution function of (3.1) is given by 
(see [4]) 
D-“(ié), 


where D(A) is the Fredholm determinant of the integral equation 
(3.2) ols DK = 
0 


Integral equations of the form (3.1) appear in other branches of pure and ap- 
plied mathematics and it is at least amusing to contemplate solving them by 
building a corresponding receiver and determining the distribution function ex- 
perimentally. 


4. Zeros of random functions. This is an extremely interesting and difficult 
problem where much further work needs to be done. If x(t¢) is a stationary 
Gaussian process with power spectrum A(w), the average number of zeros per 
unit time is given by Rice’s formula [5]: 

1/2 


f (w) dw 


f (w)dw 


The fact that by counting zeros one can obtain information about the spectrum 
is in itself of great practical interest. It has been, for instance, applied to 
turbulence by H. W. Liepmann and his group at the California Institute of 
Technology. A more detailed study of the distribution of zeros of random func- 
tions for even the simplest processes encounters great analytical difficulties. A 
closely related problem is the following: 

Let 


(4.1) 


(4.2) = ay COS + ox) 


and assume that the frequencies \, are rationally independent. Let N(T, a) be 
the number of roots of 


=a 


in 0OStST. 
It can then be demonstrated [6] that 


N(T, 
(4.3) lim lim f = {TI II Jo(axn) I + dan?) dédn. 
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This is the counterpart of (4.1) and can be interpreted by saying that the aver- 
age distance between consecutive a-values of f(t) is given by the inverse of the 
expression on the righthand side of (4.3). If one asks for the average of the 
square of the distance between consecutive a-values one runs into difficulties 
which, at least at present, appear insurmountable. This is already true for the 
simplest case 


= a, COS Al + COS 


(except for special values of a;, a2 and a). An interesting application of (4.3) to 
the theory of unimolecular reaction rates was made recently by N. B. Slater 
[7]. In conclusion let me mention another related problem, this time with no 
practical implications. 

Consider a polynomial of degree n 


(4.4) xu 
k=0 


whose coefficients are independent, normally distributed, random variables each 
having mean 0 and variance 1. It is then easy to show [8] that the average num- 
ber of real roots of (4.3) is asymptotically 


(4.5) = log n. 
Moreover, a tedious but rather elementary calculation shows that the standard 
deviation about the mean is of lower order and consequently it is very rarely 
that a random algebraic equation of high degree wil! have a number of real roots 
which is significantly different from (4.5). 
These conclusions remain valid for a much wider class of independent ran- 


dom variables [9] but proofs become enormously more tedious. For the simplest 
case 


Prob. {X, = 1} = Prob. {X, = — 1} =} 


the proof that (4.5) is still asymptotically the average number of real roots is 
lacking! 

The present exposition was naturally permeated with probabilistic consider- 
ations. But perhaps it is not too idle and inappropriate to contemplate here the 
possibility of a statistical approach to various questions in pure mathematics. 
As an example let me consider the following question: how good is the classical 
Descartes’ rule of signs? 

As applied to an individual algebraic equation the question is largely mean- 
ingless. Interpreted statistically it can be properly formulated and answered. 
The average number of changes of sign in (4.4) is 2/2 (this is to be compared 
with 2! log m which is the average number of positive real roots of (4.4)). 
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y However, one can do better. If one considers the polynomial 


the number of real roots is the same as for (4.4); the average number of changes 
of sign can now be shown to be of the order C+/n. This is about the best one can 
do, and yet C-+/n is still so far from the correct order z—! log n that we must con- 
clude that Descartes’ rule of sign is extremely unlikely to give a good estimate 
for equations of high degree. 
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BOOLEAN ALGEBRAS IN ELECTRIC CIRCUIT DESIGN 
D. E. MULLER, University of Illinois 


A brief description of how Boolean algebras may be used in the design of 
computer circuits is followed by an account of studies in the field of Boolean 
algebras which were carried out with the aid of the ILLIAC computer at the 
University of Illinois. The computer was used to investigate relationships be- 
tween functions in Boolean algebras which uncovered symmetries leading to the 
development of certain theorems. These theorems turned out on later expansion 
to have a number of applications to computer design. Hence, the work under 
discussion is quite distinct from the customary use of the computer in carrying 
out massive algebraic manipulations with the objective of settling some specific 
design problems. 

One of the main motives for the introduction of algebraic techniques into 
the design of switching circuits is the problem of devising the simplest possible 
circuit toward a specified end under a given criterion of simplicity. In principle 
this can, of course, always be achieved by brute force since it can be reduced to 
the trial of a finite, although generally very large number of combinations. A 
mathematically satisfying solution, however, has not yet been found and will 
certainly require a deeper insight into the basic structure of Boolean functions 
than available at present. 

The investigations are concerned with this latter type of functions and 
therefore apply only to circuits of the so-called static type in which the outputs 
are expressible as Boolean functions of the inputs. This excludes, in particular, 
digital circuits with feed-back loops. Any Boolean function of p independent 
variables can be expanded into a canonical representation consisting of a linear 
expression in 2? standard functions with coefficients which are either zero or 
one. The set of coefficients serves to completely define the function being repre- 
sented, and conversely any ordered set of 2” zeros and ones represents a Boolean 
function. Hence every such function can be regarded as corresponding to a 
vertex of the 2?-dimensional unit cube, and the geometric properties of the 
latter should be reflected in algebraic properties of the corresponding totality of 
Boolean functions. A particular orthogonal coordinate system, with its origin 
at the center of the cube, is introduced by passing axes through the vertices cor- 
responding to the independent variables themselves and through the points cor- 
responding to all functions which can be obtained from the latter, using only 
the operation “exclusive or.” 

Coordinates of points with respect to these axes for the case p=4 were 
studied by the computer. For each point the number of coordinates of absolute 
magnitudes 0, 1, 2, - - - , 8 were obtained, and all points for which these nine 
numbers agreed were put into the same equivalence class. Such equivalence 
classes represent sets of functions bearing similar relationships to the axes and 
therefore, presumably, similar algebraic form. There were’ 8 such classes of 
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which 3 attracted particular attention, because the number of functions con- 
tained in them was 2048, an exact power of two. Jointly they contained all points 
corresponding to functions, each one of which can be expressed as a sum of 
products of no more than two of the independent variables, and it turned out 
that the squared distance between any pair of such points is never less than 4. 
This suggested giving a special name to the set of all points corresponding to 
functions which are the sum of products of no more than p—g of the inde- 
pendent variables: this is a net of order q. It can then be proved that the distance 
between any pair of points in a net of order q is at least 2%. 

An unexpected application of this concept was found in the field of error 
detection which at the same time suggested a number of problems. In this case 
one wishes to find the most numerous set of vertices of an m-dimensional cube, 
any pair of which is separated by at least the distance d. The question whether 
for n=2?, d=2%, the net of order qg furnishes a solution has not been com- 
pletely answered. 

The simpler question of whether additional vertices can be found such that 
the distance property remains valid for the augmented set has been answered 
affirmatively for large enough p except when g=0, 1, p—1, and p. For these latter 
values of g, the net cannot be augmented for any value of ». The addition of 
further points is known to be possible when p=7, g=4. For p=5, g=3, a com- 
ot puter study has shown that no augmentation of the net is possible. The case 
=6 remains open. 
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COMPUTATIONAL ATTACKS ON DISCRETE PROBLEMS 
S. S. CAIRNS, University of Illinois 


This is a summary of some recent work done at the University of California 
at Los Angeles under the so-called SCAMP project, of which I was chairman 
during the summers of 1952 and 1953. This project represents the first organized 
effort to attack discrete problems with the aid of large-scale computing equip- 
ment. It is a pioneer effort and only in the early stages of its expected develop- 
ment. The objectives of the computational procedures are (1) to obtain numer- 
ical information useful in the development of general theories, (2) to verify 
hypotheses which are beyond the reach of hand computations, and (3) to ac- 
cumulate data on the basis of which to make conjectures as foundations for 
further research. 

Thus far much, and perhaps too much, attention has been devoted to the 
existence problem for finite projective planes of order m. This can be inter- 
preted as the problem of finding the incidence matrix A of such a plane; that is, 
a matrix A of order n=m*+m-+1, whose elements are all zeros or ones, there 
being exactly m+1 ones in each row and exactly one column in which both 
members of any given pair of rows have a one. The current state of this problem 
is the following: Finite projective planes are known to exist whenever m is the 
power of a prime and are known not to exist provided that m=1 or 2 mod 4 and 
that the square-free part of m has a prime factor congruent to 3 mod 4. In the 
particular case m =6 this had previously been shown by a direct attack, relying 
partly on numerical methods.* The question remains therefore unsettled for a 
sequence of values of m beginning with 


m = 10, 12,15, 20---. 


The smallest of these is sufficient to present a formidable challenge, yet not so 
terrifying as to discourage an attack which combines mathematical ingenuity 
with the capabilities of high-speed automatic computing equipment. 

The three basic methods of approach which have been tried are character- 
istically available in dealing with any problem of this kind, an ever-increasing 
number of which are presenting themselves in such connections as the theory of 
switching circuits and communication networks, the scheduling of production 
and transportation, as well as other assignment and allocation problems. Identi- 
fying them as the combinatorial, the algebraic, and the analytic approach, we 
shall illustrate each of them by summarizing how far it has carried towards 
solving our problem. 


The combinatorial approach. The aim will be that of laying down principles 
or carrying through constructions which permit on the basis of succinct criteria 


* G. Tarry, “Le Probléme des 36 Officiers” Comptes Rendus, vol. 1, 1900, pp. 122-123 and 
vol. 2, 1901, pp. 170-203. 
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(or a relatively small number of trials) the elimination of large classes of candi- 
dates from the frequently astronomical number of possible ones. In the present 
instance, a matter of proper labelling makes it possible to reduce the problem- 
atical portion of A to a right lower corner which consists of an m Xm array of 
mth order permutation matrices ¢,, (p,¢=1, 2, - - - , m), subject to certain con- 
ditions resulting from the definition of a finite projective plane. It is possible, 
moreover, to arrange it so that this m Xm array of matrices will be bordered on 
top and on the left by unit matrices. There now exists a suitable way of mapping 
mth order permutation matrices on column vectors. On assembling the column 
vectors thus corresponding to Gy, Cp2,°°**, Com there results for each 
p=2, 3,---,m, a latin square. A set of (m—1) latin squares arising in this 
manner from the incidence matrix of a finite projective plane consists of mu- 
tually orthogonal ones; that is, no ordered pair of elements occurs in identical 
positions more than once in any ordered pair of squares. This reduces our prob- 
lem to the existence problem for a set of (m—1) mutually orthogonal latin squares of 
order m, thus immersing it into the main stream of a large and important family 
of combinatorial problems which depend for their solution on the construction 
of permutation matrices with particular properties (e.g., that of maximizing a 
suitable function) and which are encountered in the design of experiment as 
well as in assignment problems and their variants. 

Tarry (loc. cit.) showed that for m=6 there is not even a pair of such latin 
squares; Euler had already conjectured the same for any m=2 mod 4, but this 
is still an open question. Even the existence problem for a pair of mutually 
orthogonal ones of order 10 has resisted all efforts to date: programming the 
method of exhaustive search, which Tarry carried out by hand in the case m=6, 
on a modern electronic computer shows it to be beyond the scope of current 
machines in spite of all improvements thus far devised; and the systematic 
construction of latin squares most prone to possess orthogonal companions has 
not as yet produced one which actually does. 


The algebraic approach. An incidence matrix A can be characterized as a 
solution with non-negative integral elements of 


(1) AAT = mI, + Un, 


where J, is the unit matrix and U, has all its elements equal to 1. Relaxing the 
requirement that the elements of A be integral, one can find solutions of (1) with 
rational elements, satisfying at least some of the conditions which define an in- 
cidence matrix. For example, this definition can be satisfied as far as the first 
(m+1) rows are concerned. As yet, however, there is little one can do to exploit 
these rational solutions in order to come nearer to the ultimate goal. Experi- 
ments were conducted, using the machine, in order to get information on the 
structure of a matrix G which would transform a rational solution S into an 
incidence matrix. If 


T = SG 


= 
a 
A 
é 
di 
| 
9 
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where S solves (1) and is normal, then T is also a normal solution of (1) provided 
G is orthogonal and its row and column sums are either all +1 or all —1. If T 
is integral it is an incidence matrix. Matrices G were therefore studied in cases 
where S and T are known, in the hope of finding ways of characterizing a man- 
ageably small class of possible G’s applying also in unknown cases. So far this, 
too, has proved fruitless. 


The analytic approach. This mode of attack consists in imbedding the dis- 
crete set of alternatives into a space of continuous variables. What is generally 
wanted is a procedure which, starting from an arbitrary point, lays down a 
“search path” on the now continuous manifold of candidates, leading to the 
desired object. For a computational approach there is required, in addition, a 
decision procedure when to break off such a search and from where to start 
anew, in case the original choice of a starting point has turned out to be a very 
bad one. In the case of our incidence matrix for the finite projective plane of 
order m, the conditions which define it characterize certain subspaces in the n?- 
dimensional space of all matrices of order » =m*+m-+1; and the task becomes 
that of determining whether or not the intersection of these subspaces is empty. 
Hence the type of search procedures required is of the sort that, if it has failed 
to yield a common element at a certain stage, then it follows that there can be 
no such element. Computationally it has turned out more convenient to con- 
sider the existence of 


(A = (m+ V+ | 


in place of A, where Y is characterized as follows, if A is an incidence matrix: 


, 


(a) Yu = = 1, 

1 Vm +1 
(c) Y is orthogonal. 


As yet it has not proved possible to determine whether or not these three sub- 
spaces have a non-empty intersection for m = 10. 


ON CHARACTERISTIC VALUE PROBLEMS IN HIGH ORDER 
DIFFERENTIAL EQUATIONS WHICH ARISE IN STUDIES ON 
HYDRODYNAMIC AND HYDROMAGNETIC STABILITY 


S. CHANDRASEKHAR, University of Chicago 


1. Introduction. Recent studies in hydrodynamic and hydromagnetic 
stability have disclosed the existence of a class of characteristic value problems 
in differential equations of high order—orders as high as twenty-four have been 
encountered—which appear to have a genuine mathematical interest. A partial 
list of these problems will be found in the Appendix. These problems have been 
solved, singly, as they have arisen in the physical connections. But one feels 
that there must be a general theory which embraces them all though such a 
theory is lacking at the present time. It is the object of this paper to bring these 
problems and the methods which have been developed for their solution to the 
attention of the mathematicians in the hope of stimulating their interest. 

In the solution of the problems listed in the Appendix, methods of two differ- 
ent kinds have been found useful. In the first of these methods a variational 
procedure is developed which in the manner of its application is different from 
the usual ones. The second of these methods is based on expansion in orthogonal 
functions but again applied in an unusual manner. We shall illustrate the prin- 
ciples underlying these methods by considering one of the simpler problems 
listed in the Appendix. 

2. A typical problem. G. I. Taylor's [1] investigation in 1923 on the stability 
of viscous flow between two concentric rotating cylinders provided the first 
example of a case of hydrodynamic instability for which the criterion was theo- 
retically predicted and experimentally verified. The mathematical problem un- 
derlying this classic investigation in hydrodynamic stability is the following: 

With certain simplifications suitable to the circumstances under which the 
experiments were performed, the problem requires the solution of the sixth 
order equation 


(1) (D? — a*)*»y = — a®T(1 + az)o, 
with the boundary conditions 
(2) v = (D? — a’)v = D(D*? — a*)v = 0 


for z=0 and 1, where T is the characteristic value parameter, D=d/dz and a 
and a are assigned (real) constants. In the physical problem a is negative and 
one is interested in the range 0D2a2= —3.0. 

Equation (2) provides six boundary conditions (three at z=0 and three at 
z=1) and the requirement that a non-trivial solution of equation (1) satisfy 
these conditions will lead to a determinate sequence of possible values for T 
(for given a? and a). Among these possible values of T (for given a and varying 
a’) there will be a smallest (positive) value; and in the physical problem particu- 
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lar interest is attached to the minimum of these smallest positive values of T 
as a function of a? (for fixed a). 

In discussing the characteristic value problem presented by equations (1) 
and (2), we shall distinguish the cases a=0 and a0. The problem has basically 
different characters in the two cases: in the former case (a =0) it is “self-adjoint” 
in some sense (yet to be defined!) while in the latter case it is not. Because of 
this difference the solution of the problem, while it can be made to depend on an 
extremal principle in the case a=0, it cannot be so accomplished in the case 
a0. As far as the general solution of equations (1) and (2) is concerned, a 
separate detailed discussion of the case a=0 would hardly seem justified; but 
it does happen that the case a=0 is the simplest proto-type of a wide class of 
problems to which similar methods of solution can be applied. For this reason 
we shall treat the case a=0 at some length. 


3. The case a=0: example of the variational procedure. When a=0, the 
solution of the characteristic value problem can, in principle, be achieved very 
simply. For, in this case, the solution of the equation 


(3) (D* — a*)*» = — a*Tv 


must clearly be of the form 
6 

(4) = Aver, 
t=1 


where the q,’s, occurring in pairs, are the roots of the characteristic equation 
(5) = — 07, 


and the A,’s are constants of integration. The requirement that the solution 
represented by equation (4) satisfy the boundary conditions (2) will lead to a 
system of six linear homogeneous equations; and the determinant of this system 
must vanish if we are not to have the trivial solution A;=0,i=1, - - - , 6. And 
the condition that the determinant vanish will provide an equation for deter- 
mining T. There is, of course, no difficulty of principle in carrying out this pro- 
cedure. Indeed, it has been carried out by Pellew and Southwell [2] who find for 
example that 


(6) T = 1707.8 for a = 3.12. 


(This is the value of a at which T, as a function of a, attains its minimum.) 

For the case on hand the carrying out of the direct method of solution de- 
scribed in the preceding paragraph is particularly simple as the roots of equation 
(5) can be explicitly written down in terms of the known cube roots of —1; 
thus 


(7) = a + w¥(a*T), = — 1). 


; 
5 
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But this simplification does not arise in other problems which have to be con- 
sidered ; also, the characteristic equations which have to be solved are of higher 
order. It would therefore be useful if we can devise a method of solution which 
will avoid the necessity of solving complicated algebraic equations of high 
order by laborious methods of trial and error. This would be specially important 
if, as it often happens, we are required to establish the dependence of the char- 
Fe, acteristic values on two or more parameters rather than on a single parameter, 
20 such as a, in the present problem (a=0). And it does appear that in many of 
the cases which arise, convenient variational methods can be devised which we 
shall now illustrate by considering the problem presented by equations (2) and 


(3). 
Letting 
(8) G = (D? — a*)v and F = (D? — a*)G = (D* — a?*)2, 
we can rewrite equation (3) in the form 
(9) (D? — a*)F = — a°T», 
while the boundary conditions (2) require that 
(10) »v=G=DG=0 for z = 0 and 1. 


Let 7; denote a particular characteristic value and let the various functions 
derived from the solution v; belonging to 7; be distinguished by a subscript 7. 
Multiplying the equation governing v; by G; belonging to a different character- 
istic value 7; and integrating the resulting equation over the range of z, we ob- 
tain 


(11) ‘GAD? — a )Fydz = — f — a’)v;dz. 


The integral occurring on the left-hand side of this equation can be reduced by 
two successive integration by parts; thus, 


(12) f "GAD" — a*)Fids = [G;DF; — (DG))F;). + f "RAD? — dz. 


The integrated parts vanish in virtue of the boundary conditions (cf. equations 
(10)) and we are left with 


(13) — dz = f ids. 


Similarly, after an integration by parts the right-hand side of equation (11) be- 
comes 


(14) a°T; [(Dv;)(Dv;) + 


J 
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Thus, 
1 1 
= aT; i i 
(15) FF jdz = [(Dv,)(Dv;) + 


Noticing that interchanging 7 and j in this equation replaces T; by T; but leaves 
the equation otherwise unaffected, we conclude that 


1 
(16) f FF jdz = 0 ij; 
0 


and, further, that when i=j7 the corresponding characteristic value, T, can be 
expressed as the ratio of two positive definite integrals in the form: 


1 
f 
0 


a? f [(Dv)? + ]dz 


(17) T= 


These facts clearly imply the existence of an extremal principle which can be 
made the basis of a variational method of solving the underlying characteristic 
value problem. For comparison with the method we shall describe in §4 for the 
case a0, we shall formulate the extremal principle in the context of a slightly 
transformed equation. 

Operating on equation (3) by (D?—a*) we get 


(18) (D? — = — 
as the equation governing 
(19) W = (D? — a?)o. 


According to equations (2) and (3) the corresponding boundary conditions on W 
are: 


(20) W=DW=0 @)G = (D?— forz = Oand i. 
Again rewriting equation (18) in the form 
(21) (D? — = — a®TW, 


and multiplying the equation governing W; (belonging to T;) by F; (belonging 
to T;) we obtain 


(22) f 'F — a*)Fydz = — a®T; f W(D*? — a*)G;dz. 
0 


Making use of the boundary conditions (20) we can reduce equation (22) by 
one or more integrations by parts to the form 


| 


a 
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(23) [(DF;)(DF;) + a°F F;]dz = 

From this equation it follows that 

(24) = 0, ix 
and that when 1=j we can express T as: 


f [(DF)? + a?F?]dz 
0 


1 
a* f 
0 


This last formula, like (17), expresses T as the ratio of two positive cefinite 
integrals. 

Consider now the effect on T (evaluated according to equation (25)) of an 
arbitrary variation, 5W, in W compatible only with the boundary conditions 
on W. We find in a straightforward manner that 


(25) T 


1 
a? f G*dz 


(26) f {(D* — a°)F + dz, 


where it should be noted that in the reductions leading to equation (26) the rela- 
tions implied in the definitions of F and G (cf. equation (20)) have been used; 
in particular 


(27) = (D* — 


Hence, to the first order, 57=0 for all small arbitrary variations in W which 
satisfy the boundary conditions, provided 


(28) (D? — a*)F = — @TW, 


t.e., if the differential equation governing W is satisfied. It is evident that the 
converse of this proposition is also true. Further, it follows from (25) that the 
true solution of the problem (belonging to the lowest characteristic value T) 
leads to the minimum value (in the sense of the calculus of variations) for T 
when evaluated according to (25). 

While the minimal principle formulated as above is valid for all small arbi- 
trary variations 5W compatible only with the boundary conditions on VW, it is 
true also for all small arbitrary variations 6F, the variation, 5W, in W being de- 
termined in terms of 6F by the equation 


(29) (D? — = 


= 
- 
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and the boundary conditions 
(30) = DiW = 0 forz = Oand 1. 


The validity of this minimal principle in this more restricted form actually 
provides a more effective basis for a variational procedure for determining T. 
For on this latter basis the procedure would be the following: 

Assume for F an expression involving one or more parameters, A,, and which 
vanishes at z=0 and 1. With the chosen form of F solve the equation 


(31) (D? — a*)*W =F 


for W and arrange that the solution satisfies the boundary conditions W=DW 
=0 for z=0 and 1; since equation (31) is of the fourth order there will be just 
enough constants of integration to do this. With W determined in this fashion 
evaluate T according to (25) and minimize it with respect to the parameters 
A,. In this way we shall obtain the “best” value of T for the chosen form of F. 

In practice it is found that even with the simplest trial functions, the varia- 
tional procedure in the foregoing form gives surprisingly high accuracy in the 
deduced values of T. Thus, with the trial function, F=sin 7z, with no variational 
parameter, the method gives T=1715.1 for a=3.12; while the function, F 
=sin rz+A sin 3rz, with one variational parameter leads to T =1707.9 for the 
same value of a; these values obtained in the “first” and the “second” approxi- 
mations should be compared with the value T = 1707.8 obtained from an “exact” 
solution of the problem. The origin of this high precision in the deduced values 
of T must clearly be traced to the fact that in satisfying four of the six boundary 
conditions of the problem (namely W=DW=0 for z=0 and 1) we have exactly 
satisfied the second of the pair of differential equations, 


(32) (D? — = —a?TW and (D? — a*)*W =F, 


which governs the problem. 

[It may be noticed here that had we used equation (17) (instead of (25)) as 
the basis of the variational method, we should have had to assume a form for 
DG such that not only does it vanish at z=0 and 1, but also G, obtained after 
integration, vanishes at the same points; further, W will have to be obtained 
as the solution of (D?—a?)W=G which vanishes at z=0 and 1. Apart from 
these differences in detail the application of the variational method based on 
(17) proceeds along essentially the same lines. But the method based on (25) 
is preferable since it avoids the restrictions implied by the requirement that both 
G and DG vanish at z=0 and 1.] 


4. The case a~0: example of a method based on expansion in orthogonal 
functions when no extremal principle exists. Returning to equation (1) and the 
general case a0, we can readily verify that propositions similar to those em- 
bodied in equations (15), (16) and (17), for the case a=0, cannot be established 
now. Thus, if equation (1) and the boundary conditions (2) together represent 
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a system which might be described as “self-adjoint” when a=0, it cannot be so 
described when a0. It is important to note that this conclusion in no way 
depends on the sign of a; in particular it is independent of the circumstance 
that when a< —1 the operator whose characteristic values we are seeking be- 
comes singular on account of (1+az) having a zero in the range 0 S21. 
Since equations (1) and (2), in general, do not allow the formulation of an 
extremal principle on which a variational procedure might be based and since 
the general solution of equation (1) cannot also be written down in any con- 
venient form, it would appear that the only remaining course is to expand the 
unknown functions in Fourier series. In Taylor’s original investigation this was 
the method which was adopted: Taylor expanded » in a sine series of the form 


(33) v= V, sin 
n=1 


and obtained T as the characteristic root of an infinite matrix. But the process of 
determining the root was not a very convergent one; and he succeeded (with 
considerable effort) in determining the characteristic root for only one value of 
a< —1. However, we may expect that a more rapidly convergent process will be 
obtained if in using the method of expansion in Fourier series, we incorporate in 
the method the same basic idea which led to the high precision of the variational 
method for the case a=0. We shall now indicate how this can be accomplished 
(see Chandrasekhar [3]). 
First we transform equation (1) by rewriting it in the form 


(34) (D? — = — a*T», 
and operating on it by (D?—a?); in this way, we obtain the differential equation 
(35) (D? — a?) { (D? — owt = —a TW, 
1+ az 
for 
(36) W = (D? — a?)». 


The corresponding boundary conditions on W are the same as those given in 
(20). 
Equation (35) is equivalent to the pair of equations 


(37) (D? — a®)*W = (1+ az)y 
and 
(38) (D? — a*)y = — a®TW, 


while the boundary conditions (20) can be expressed alternatively in the forms 


1 
By 
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(39) W = DW=0 and y=0 for = 0 and 1. 


Since y has to vanish at z=0 and 1, we can expand it in a sine series of the 
form 


(40) y= > C, sin 


n=l 


Having chosen y in this manner, we next solve the equation 
(41) (D? — a*)?W = (1 + az) >. C, sin ns, 
n=1 


obtained by inserting (40) in (37), and arrange that the solution satisfies the 
four remaining conditions on W; since equation (41) is of the fourth order there 
will be just enough constants of integration to do this. With W determined in 
this fashion and y given by (40), equation (38) will lead to an infinite determi- 
nant which must be zero if all the C,’s are not to vanish. In this way we shall ob- 
tain a characteristic equation for determining T. 

When the details of the method described in the preceding paragraphs are 
carried out, one finds that the process of solving the infinite order characteristic 
equation for T, by setting the determinant formed by the first ” rows and col- 
umns equal to zero and letting m take increasingly larger values, converges very 
rapidly indeed. Thus for a= —2.5 and a=5.00, 5.05 and 5.10 the values of T 
obtained in the third and the fourth approximations (the “order” of the ap- 
proximation being the order of the determinant which is set equal to zero in the 
determination of T) are: 


5.00 4.607 X 10 4.626 X 10‘ 
(42) 45.05 T (3rdapp.) = {4.600 X 104 T (4th app.) = 14.619 104. 
5.10 4.604 X 4.623 X 10¢ 


It is seen that the values of T given in the third and the fourth approximations 
differ by only four parts in a thousand. The origin of this rapid convergence 
clearly lies in the splitting of the original equation of order six into a pair of 
order two and four respectively and satisfying the equation of order four 
exactly. This basic idea underlying the method is capable of extension and ap- 
plication to a wide class of problems. 


Appendix 
We shall list here some of the more important characteristic value problems 
in high order differential equations which have occurred in recent studies on 
hydrodynamic and hydromagnetic stability. The particular physical connections 


in which they arise are indicated. The references are to papers in which the solu- 
tions of the problems will be found. 


40 SYMPOSIUM ON SPECIAL TOPICS IN APPLIED MATHEMATICS 
I. The inhibition of convection by a magnetic field [4, 5]. 
1) To solve 
(D? — a?)[(D? — — QD?|]W = — a?RW, 
together with the boundary conditions 
W = [(D? — a’)? —QD*|W =0 on z= +}, 
and 
either DW =0 on z= +}, 
or DW=0 on z= +}, 
or DW =0 on and DW=0 on z= —}, 


where a and Q are assigned positive constants and R is the characteristic value 
parameter. 

The physical problem requires the minimum (with respect to a) of the lowest 
characteristic value R for various assigned values of Q. 

2) To solve 


(D? — a?)[(D* — — Q(D + ia)?|W = — a RW, 
together with the boundary conditions 
|W| =| (D? — — Q(D + ia)*W|=0 on z= +}, 
and 
either |DW|=0 on z= +}, 


or |[DW|=0 on s= +}, 

or |DW|=0 on and |D*W|=0 on z= —}, 
where a, @ and Q are assigned constants and R is the characteristic value 
parameter. 


The physical problem requires the minimum (with respect to a and @) of 
the lowest characteristic value R for various assigned values of Q. Since W is 
complex, we have here a genuine problem in an equation of order twelve. 

II. The instability of a layer of fluid heated below and subject to Coriolis ac- 


celeration [6]. 
3) To solve 
(D? — = — (= i) DW 
and 
(D? — a?) G — — i) vz| 


= [(D? — + TD?|W = — a*RW, 


j 
- 
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together with the boundary conditions 
22 
W= (= Dz = 0 on s= + }, 
v 
and 
eithr DW=Z=0 on z= +}, 
or DW=DZ=0 on z= +}, 
or DW =Z=0 on and DW=DZ=0 on —}, 
where a, v, d and 2 are assigned constants, T= 40%d‘/y? and R is the character- 
istic value parameter. 
The physical problem requires the minimum (with respect to a) of the low- 


est characteristic value R for various assigned values of T. 
3a) To solve 


(D? — a? — = — (= a) DW, 


22 
— a*)(D? — a? — ic)W — (= =F 
v 
and 
(D? — a? — iao)F = — Ra’W, 
together with the boundary conditions 
|F[=|W|=0 on +}, 
and 
eithr DW=Z=0 on s= +}, 
o DW =DZ=0 on +}, 
or DW=Z=0 on and DW=DZ=0 on «= —}, 
where @ is a further assigned constant and o (which is real) is a parameter to be 
determined (for given T=40?d‘/y? and a*) by the condition that R is real; and 
the physical problem requires the minimum (with respect to a?) of these real 
values of R for various assigned values of T and @. 
III. The instability of a layer of fluid heated below and subject simultaneously 


to a magnetic field and Coriolis acceleration [7]. 
4) To solve 


(D* — at) [{(D* — — + ]W = — Ra*[(D* QD*|W 


together with (several different) five pairs of boundary conditions on W at 


‘ 
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2= +4, where a, Qand T are assigned constants and R is the characteristic value 
parameter. 

IV. The stability of viscous flow between rotating cylinders [1, 3, 8, 9]. 

5) To solve 


(D*? — = — a*T(1 + az)o, 
together with the boundary conditions 
v = D*y = D(D*? — a?)) = 0 on 2 =0 and 1, 


where a and a (<0) are assigned constants and T is the characteristic value 

parameter. The physical problem requires the minimum (with respect to a) of 

the lowest (positive) characteristic value T for various assigned values of a(<0). 
6) To solve 


(DD* — = (4 + 


v? 
where D=d/dr and D*=D+1/r, together with the boundary conditions 
v = (DD* — )*)v = D*(DD* — = 0 for r = R, and 


\, v and B are assigned constants and A is the characteristic value parameter. 
V. The stability of viscous flow between rotating cylinders in the presence of a 
magnetic field [10]. 
7) To solve 


[(D? — a)? + Qa*]*y = — Ta*(D* — a*)y, 
together with the boundary conditions 
Dy = = on +4, 


where a and Q are assigned constants and T is the characteristic value param- 
eter. The physical problem requires the minimum (with respect to a) of the 
lowest characteristic value T for various assigned values of Q. 

VI. The stability of viscous flow between rotating cylinders in the presence of a 
radial temperature gradient [11]. 

8) To solve 


together with the boundary conditions 


dr dr? r dr r? a ‘ 


where m is an integer and S, is the characteristic value parameter. 


if 
hy 

dr? r dr fr vail 
a 
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9) To solve 
a? W 
W = — 
dr? dr 


together with the same boundary conditions as in (8) above. 
VII. The thermal instability of fluid spheres and spherical shells [12, 13, 14]. 
10) To solve 


)w = 


dr? r dr r? 
together with the boundary conditions 
da? 2d Kki+1)\? 
W =(—+-— —- W=0 
(a) (= r dr r? ) 
and 
dw aw 
either = 0 forr = 1, 
r dr? 
and 
W = as +0; 
2d Ki+1)7? 
(b) W = lw =o for r = 1 and n(<1), 
dr? dr? r dr r? 


where / is an integer and C; is the characteristic value parameter. 
VIII. The stability of superposed fluids (plane problem) [15, 16, 17, 18, 19]. 
11) To solve 


+ + Dw + — + Dw} 


together with the boundary conditions 


w=0 for and 


and 
either Dw=0 for -=0 and J, 
or for and I, 
or Dw=0 for and D*'w=0 for z=1, 


or D’w=0 for and Dw=0 for 


i+ 
q 
j 
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where p=p(z) and u»=,(z) are given functions of z, k is an assigned (real) con- 
stant and ¢ is the characteristic value parameter. (Note that o can be complex). 
IX. Tie stability of superposed fluids (spherical problem) [20]. 
12) To solve 


Ad, Som} + (i+ 2)(0 1) 


dr o dr drt dr? r? 
d du (W 
r o o dr dr dr\r 


where 


y, 


dr? or dr 
together with the boundary conditions 


dw aw 
W=0 and either — or 
dr dr? 


and 
W = O(r') as r—0. 


In the foregoing p, u and ¥ are given functions of r, / is an integer and o (which 
may be complex) is the characteristic value parameter. 
The case when 


= = constant, p = p1 = constant forO Sr <7; 


and 
= pe = constant, p = p2 = constant forr; <r ZR 


is of particular importance. In this case the equation reduces to the form 


ad’? (i+ 1 

dr? po r po 
in each of the two regions of constant p and wu. And solutions are sought which 
satisfy the boundary conditions 


vw 
= forr=R 
dr? 


W=0 and 


W = O(r') forr—0 


| 
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and 
aw 
dr? 


dw 
W, and + 2)0 nw} 


are continuous on r=7;. 
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FREQUENCY SPECTRUM OF VIBRATIONS 
OF A CRYSTAL LATTICE 


ELLIOTT W. MONTROLL, University of Maryland,* and Office of Naval Research 


1. Introduction. The heat capacity of a system of N interacting particles 
under an external constraint x is defined as C,= (0Q/0T),, the ratio of the added 
heat to the temperature rise of the system. Experimentally it is easiest to meas- 
ure the heat capacity of solids at constant pressure, C,, and of gases at constant 
volume (in a closed cor-tainer), C,. The quantities C, and C, are related through 
the thermodynamic equation C,—C,=9a?VT/x, a being the coefficient of 
thermal expansion, V the volume, and x the compressibility. Since statistical 
mechanical formulae lead naturally to C,, one must include the correction due 
to lattice expansion in applying theoretical results to the interpretation of ex- 
perimental C, measurements in solids. 

The modern theory of the heat capacity of solids was originated by Einstein 
[1] in 1907 through the application of the then new Planck quantum theory. 
The appropriate basic f»rmulae are those for the internal energy and heat ca- 
pacity of N independent harmonic oscillators of frequency v: 


1 1 
(1.2) C, = =) / sink? =). 


Einstein assumed that each atom in a solid vibrates about its equilibrium posi- 
tion independently of its neighbors and used equation (1.2) to fit the then avail- 
able heat capacity data. Equation (1.2) is a universal function of the variable 
kT /hv. Hence it was predicted that all experimental C, curves, when plotted on 
an appropriate scale, could be superimposed. The frequency »v determined the 
scale. The anomolously low value of the atomic heat of diamond was explained 
by stating that its frequency » was larger than that of most materials. 

The general form of experimental C, curves is given in Figure 1. Experimental 
deviations from the Einstein formula are to be expected at low temperatures 
(and indeed have been observed for years) because atoms in a solid do not vi- 
brate independently of each other. Interatomic forces are such that the oscilla- 
tions of a given atom are coupled with those of its neighbors. Hence, the 3NV 
equal frequencies of equation (1.2) should be replaced by the set of frequencies 
of the normal modes of the coupled oscillating atoms. 


* The preparation of this article was partially supported by the Office of Scientific Research 
of the U. S. Air Force. 
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It is useful to introduce a distribution function g(v) which has the property 
that g(v)dv is the number of normal modes between v and »y+dy. Then equations 
(1.1) and (1.2) become 


and 


Here vz is the largest normal mode frequency. Optical as well as thermodynamic 
properties of crystals can be discussed in terms of g(v). 


6 


5 
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Fic. 1. The variation of the specific heat (in calories per gram mol) of a typical solid as a 
function of temperature. 


The first calculation of the distribution of frequencies of normal modes was 
made by Debye in 1912. He postulated a solid to be an elastic continuum. Since 
a continuum has an infinite number of normal modes, Debye cut the frequency 
spectrum off at a frequency vy, such that the total number of normal modes was 
equal to the number of degrees of freedom of the solid. In the Debye [2] theory 
g(v) =9N(v/v1)*/v_ where vz, the largest frequency, is a function of the elastic 
constants of a solid; g(v) =0 when »y>vz,. Debye’s work was one of the great 
successes of the early quantum theory. The theoretical heat capacities based on 
his frequency spectrum are in good agreement with experimental results. The 
simplicity of the Debye theory combined with this fact has given it a long life 
as the dominant theory of the heat capacity of solids. 

The continuum model can be expected to give good results at very low tem- 
peratures where only long wave length oscillations are important. However, in 
the steep part of the C, curve it is not sufficiently accurate. A discussion of the 
experimental situation is given in a review by Blackman [3]. Born and von 
Karman [4] were the first to investigate the dynamics of periodically spaced 
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atoms. They introduced a model of a solid which is mechanically equivalent to 
a set of coupled springs and masses. 

Our main purpose here is to review the general characteristics of the fre- 
quency distribution, g(v), of the Born-Karman model. A striking feature of g(v) 
curves is the existence of certain singularities, logarithmic ones in two dimen- 
sional lattices and inverse square root singularities in the derivative of g(v) in 
three-dimensional lattices. The existence of peaks was first noticed by Blackman. 
The analytical nature of these peaks, our main subject of interest here, was first 
discussed in an investigation [5] of the vibrations in the plane of a two dimen- 
sional square lattice with interactions between atoms which were nearest and 
next nearest neighbors. The logarithmic peaks which occurred in the g(v) of this 
case were also found by Rosenstock and Bowers [6] for the transverse vibrations 
of atoms in a square lattice and finally by Smollett [7] for the in plane vibra- 
tions of the two dimensional lattices with long range Coulomb interactions. 
Similar results were obtained for the transverse vibrations in the two dimen- 
sional hexagonal graphite structure [8, 9]. 

Van Hove [10] then showed by an application of the critical point theory 
of M. Morse [11] that the existence of these singularities is a consequence of the 
periodic lattice structure and not of the detailed character of the atomic force 
constants. He discussed both two and three dimensional lattices. His theory of 
the nature of the singularities in g(v) for the three dimensional case was verified 
in specific examples by Newell [12] and Rosenstock [13, 14]. 

We shall use the abbreviations 2D and 3D in their usual connotation and 
FS to mean frequency spectrum. RL will denote reciprocal lattice. 

Many of the results to be reviewed here are applicable to the general prob- 
lem of wave propagation in 2D and 3D periodic structures. For example the 
electron energy level distribution in metals has the same behavior as the g(v) of 
lattice vibrations. 


2. Fundamental formulae for analysis of doubly periodic lattices. We shall 
introduce the discussion of the vibration of periodic lattices by analyzing some 
2D examples and then generalizing to 3D systems. 

Let us consider a square lattice which contains N rows and N columns of 
identical particles of mass M. The lattice points are identified by (l,m) (where 
both / and m range through the integers 1, 2, - - - , N). The components of the 
displacement of the particle (/,m) are designated by u:,m and ¥;,m. Our notation 
is summarized in Figure 2. 

We shall assume that our lattice is formed on a torus so that we can apply 
the simplifying Born-Karman boundary conditions: 


(2.1a) = my Vim = 


Physically one expects thermodynamic functions which are proportional to the 
number of particles in a system to be independent of surface effects in three 
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dimensions and of boundary effects in two dimensions. We shall discuss the 
mathematical theory of this point in section 5. The boundary conditions (2.1) 
are generally used in solid state problems. 


l-i,m+1 I+1,me1 
| 
| 
I- I,m I+1,m 
I-1,m-1 I+1,m-1 


Fic. 2. Designation of lattice points in a two dimensional lattice. 


When the particles are displaced from their equilibrium positions the total 
potential energy of interactions can be expanded in a Taylor series about its 
equilibrium value: 


® = — Ci — >. + terms quadratic, cubic, efc., in #1,m and 02,m- 


Since Ci,m and Cj,, represent the forces in the horizontal and vertical directions 
acting on particle (1, m) while it is at equilibrium, they must vanish, equilibrium 
being defined as that state in which no forces act on any particles. By making 
the hypothesis that displacements from equilibrium are small, we are justified 
in neglecting the cubic and higher order terms in the displacements. 

The forces which act in the horizontal and vertical directions respectively 
on the (/,m)-th particle are 


o® 
OU1,m 
o® 
O01,m 


where —4),,%14.,m+, is the contribution of the horizontal force on (/,m) due toa. 
horizontal displacement of a particle \ lattice distances to the right and yu 
above (l,m), etc. 
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The equations of motion of our system of coupled particles are 
(2.3a) = — 
(2.3b) = — 
These equations are linear in the displacements and contain only constant coeffi- 
cients to the “;,m’s and v;,m’s. Hence we can seek solutions of the form 
(2.4a) him = tei (or 
(2.4b) Vm = ve2™ tei (1 
where 4, v, v, $1, and $2 are constants which must be chosen so that (2.4) satisfies 
(2.3), the boundary conditions (2.1), and whatever initial conditions that might 
be imposed. 

The boundary conditions (2.1) are satisfied if ¢; and @ are solutions of 


= 1 and = 1, 


that is, if 

(2.5) = 2xk/N and = 2xj/N 

with k, j=1, 2,---, N or (if N is even) k, j may run through —4N+1, 
—}3N+2,---+,0,1,2,---+,4$N. This gives N? possible choices of the pair (¢1, $2). 


The frequency » is related to ¢; and ¢2 by noting that if (2.4) is substituted into 
(2.2) and (2.3) the constants u and v are solutions of the homogeneous linear 
equations 


(2.6a) $2) — M } + oF 12(¢1, 2) = 0 

(2.6b) uF ($1, 2) + 2) — = 0 

where 

(2.6c) $2) = a $2) = by, yet 


(2.6d)  Far(b1, $2) = ba) = dy pet orto), 


In order for solutions of (2.6) to exist, the determinant of the coefficients of u 
and v must vanish: 


— Fi2(¢1, $2) 
$2) F22(¢1, $2) — 


This equation relates the frequency v to the constants (1, ¢2). There are two 
values of v? for each (¢:, $2) pair. Hence there is a total of 2N? normal modes, or 
solutions (2.4) of our equations of motion (2.3). The most general solution is a 
linear combination of these 2N? normal modes. 


(2.7) 
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3. Curves of constant frequency in (¢;, ¢2)-space and distribution of normal 
mode frequencies. Equation (2.7) is a quadratic equation in the square of the 
circular frequency w=2zv: 


(3.1) — + Foo) + — = 0. 
The solutions of this equation are 
(3.2) [o(d1, bo) ]? = + Fos) + — Foe)? + 4F 21 


Since each F;,; is a doubly periodic function of (¢1, ¢2), w also has this character 
and w(¢1, ¢2) when m and n are positive or negative 
integers. We shall refer to all those frequencies which are generated as (qi, $2) 
runs through the N? values defined by (2.5) as a branch of the frequency spec- 
trum. Our 2D lattice has two branches. We shall refer to the (¢:, ¢2) space over 
which our doubly periodic function w(¢:, $2) is defined as a doubly periodic 
space. 

Henceforth all of our discussions will be concerned with lattices which con- 
tain a very large number, N?, of particles. Our various formulae will refer to 
limit results as N+. In this limit the uniformly distributed points in (¢;, ¢2)- 
space (those defined by 2.5) which correspond to a branch of frequencies of 
normal modes become dense. The number of normal modes associated with a 
closed region in (¢:, $2)-space is proportional to the area of that region. The pro- 
portionality constant is N*/47? since the area of a single period of our doubly 
periodic (¢1, ¢2)-space is associated with N? normal modes. 

It is clear that a doubly periodic function defined on a two dimensional space 
is equivalent to a function defined on a torus. Since the F;;’s are continuous 
functions of (¢1, ¢2), two neighboring points in (¢;, ¢2)-space yield frequencies 
which differ only slightly from each other. Curves of constant frequency can be 
constructed in (¢1, ¢2)-space. They are closed curves on the (gi, ¢2) torus. The 
curve w=w,=constant separates the torus into two regions, one which cor- 
responds to frequencies vy = w/2m which are less than »;=w,/2m and the other to 
frequencies greater than », (there may actually be several closed curves of a 
single frequency of interest; in this case the regions of frequencies <v; may not 
be connected). 

If we let N.(v) be the number of frequencies less than » in the a-th branch, 
we have in the limit as N- 


(3.3) Ne) = ff 


where the integration extends over the set R of all values of (¢:, $2) for which 
$2) (the definition of 2(¢1, $2) is w2 with given by (3.2)) is 
less than v*. Since the number of frequencies between v and »y+dy in the a-th 
branch is 


Za(v)dv = + dv) — N.(v) = (AN./dv)dv, 
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we have 

(3.4) ga(v) = ay op 


An alternative formula, which has been derived by Bowers and Rosenstock 
[6], in terms of 5-functions, allows one to express ga(v) as an integral overa 
complete period of (¢1, ¢2)-space. Now we write 


ga(v)dv = + — Na(v — 4dr) 


= 


where S is the set of values of (¢1, $2) for which 


| 


— vdv < ve($r, ¢2) — < vdy. 
If F,(x) is a function with the property 
1 if -a<x<a 
a(x) = { 
0 if |x| >a 
then 


f — F,(x)dx = 1, 
so that (1/2a)F.(x) approaches the delta function 5(x) as a0. Since 


N? 
= f $2) — v 


we have 


(3.5) = $2) — 


The Fourier integral representation of the 5-function 


1 
= — f 
2rd 


will sometimes be useful in the detailed calculation of (3.5). It leads to 


fex {i ty [va(dr, $2) } | 


Obvious 3D generalizations of (3.4) and (3.6) exist. 


= 
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4. The frequency spectrum of some special examples of two-dimensional 
lattices. We shall now apply the ideas developed in the two previous sections 
to the calculation of the frequency spectrum of 2D lattices whose vibrations 
are transverse (normal to the plane of the lattice) and of several special cases 
with vibrations in the plane. 

A considerable insight into the behavior of curves of constant frequency and 
the FS can be achieved by studying a simple model in which each particle on a 
2D square N XN lattice has only one degree of freedom. This we shall regard 
as the displacement of a particle in a direction normal to the equilibrium plane 
of the lattice. We shall postulate each particle to interact only with its nearest 
neighbors but in an unsymmetrical manner with the force constant for inter- 
actions along a row differing from that for interactions along a column. 

We let 2:,m be the displacement of the (/,m)-th particle from equilibrium, and 
choose the potential energy of interaction between particles to be 


I,m i,m 
The equations of motion of the particles on our lattice are 


+ + 
As usual we use the periodic boundary condition 


Zim = = = 


and find the frequencies of normal modes to be given by 


(4.2) 4n*v?M = 2(a1 + ae) — 2a; Cos — COS 
where 
(4.3) = 2xj/N, $2 = 2xk/N, and j,k =0,1,2,---,N—1. 


If N is even we can choose j and k to run from —$4N+1 to 4N. 
The largest frequency, vz, appears at ¢:=¢:=7. Then 


(4.4a) 4n = 4(a + a). 
Hence, if we let 
(4.4b) and f = v/v 


(4.2) becomes 


(4.5) (1 — 2f*)(1 + 17) = r cos di + cos 


Schematic curves of constant frequency are plotted in Figure 3. 
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We can apply equation (3.4) to the calculation of g(v). A special integral 
must be evaluated in the three cases 


i. +r) 
ii. fs i/(1t+r) 
iii Witnsfsi. 

(-1,1) >, (1,1) 


| 


(1-1) 


Fic. 3. Schematic curves of constant frequency in the reciprocal lattice of a two 
dimensional crystal. 


The typical curves in one quadrant of each of three regimes are given in Figure 
4. We now calculate the integrals of each of the shaded areas in the figures. 


T G 


% 4, 


a b c 


Fic. 4. Typical curves of constant frequency for three ranges of frequency. 


The FS is proportional to the derivatives of these areas with respect to fre- 
quency (equation 3.4). 
Case i. The intercept at ¢:=0 is 


= cos [1 
while that at ¢.=0 is 


cos [1 — 201+ /r] 


0 T 0 0 

| 
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Hence the shaded area in Figure 4a is 


(0) 
cos~! {(1 — 2f)(1+ 1r) — rcos 
0 


so that (remembering that each quadrant yields only one fourth of the frequen- 
cies) 


4N°(1 


dx 


In our interval of interest 0 <f?<r/(1+1r), the roots of the polynomial of the 
denominator of the integrand are ordered with 


Hence a direct application of equation 553 in Peirce’s Table of Integrals yields 


4N°(1 1— 
4.6) neo) it 
r 
where K(k) is the complete elliptic integral 


a dx 
K(k) = [(1 — hx?) ]1/2 


Case ii. Here we compute the shaded area in Figure 4b. This is 


ff {(1 — 2f*?)(1 + 1r) — cos o} dd. 


Hence, if we let x =r cos ¢, and apply (3.4) we find 


4N2 ril2 
(4.7) vig(v) = — t+ rfVi-—f* 


if fs i/(itr). 


Case iii. In the range f?>1/(1+r7), the curves of constant frequency intersect 
the line ¢.=7 at 


= cos {[(2+) +n] 


Hence the area of the shaded portion of Figure 4c is 
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+ cos {(1 — +r) — cos gi} des 


so that 
N*(1 


It is to be noted that (4.8) and (4.6) have exactly the same form. 

In a symmetrical two-dimensional lattice r=1 and the range in which (4.7) 
exists vanishes. The elliptic integral K(k) behaves as log [4/(1—k?)"/?] as k-1. 
This logarithmic singularity occurs in g(v) as v/v, approaches [r/(1+1) ]!/? or 
[1/(1+r) In fact, 


— (N*/n*)(1 + log {v/vz — [r/(1 + 7) 

as (v/yr)? > r/(1+ 17) 
— + log — (1 + 

as 1/(1 + 1). 


As the lattice becomes more and more symmetrical, (7.e., as r—1) the two 
logarithmic singularities approach each other and become one in the limit. In 
other words the singularity is split by reducing the symmetry. We have plotted 
g(v) for several values of r in Figure 5 to exhibit this property. 

The result for a symmetrical lattice agrees with that of Bowers and Rosen- 
stock who showed that the logarithmic singularity persists even when next near- 
est neighbors are taken into account. As long as the nearest neighbor interactions 
predominate, the singularity moves toward higher frequencies with increasing 
next nearest neighbor force constants. Curves based on the formulae of Bowers 
and Rosenstock are given in Figure 6. 

An especially important feature of the results given above (as well as of those 
of Bowers and Rosenstock [6]) is that the logarithmic singularities in the FS 
occur at frequencies which correspond to saddle points in (¢:, ¢2)-space. In 
Figure 3 the saddle point on the f?=(1+7)—! curve occurs at (0, 7) and that on 
the f?=r/(1+7) curve appears at (7, 0). In a symmetrical lattice with r=1, both 
of these saddle points are on the curve with f?=}. 

Since K(k)~47 as k-0, we have as 


(4.9a) 
vig(v) ~ 


(4. 9b) 


vrg(v) ~ 2N%(1 + 


In the Debye continuum theory of a two dimensional lattice g(v) is proportional 
to v for all frequencies. It is not surprising that the continuum and discrete 
models give similar results in the range of wave lengths long compared to lattice 
spacings. 

Equations (4.6)—(4.8) reduce properly to the one dimensional FS as r—40. 
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Fic. 5. Variation of frequency spectrum of two dimensional lattice with one degree of freedom 
per lattice point with increasing anisotropy. The parameter r equals 1 in an isotropic lattice and 
zero in the absence of coupling in one direction. 


Fic. 6. Variation of frequency spectrum with ratio of next nearest to nearest neighbor force 
constants in two dimensional lattice with one degree of freedom per lattice point. Solid line to left 
represents no next nearest neighbor coupling. As this coupling increases the peak moves to the 
right. The curve whose peak falls at »/yz = 1 corresponds to equal nearest and next nearest neighbor 


coupling. 
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In this limit (4.7) describes the entire spectrum. Since K(0) =42, (4.7) becomes 
2N 
nae) = py], 
the one-dimensional g(v). 


Equations 4.6—4.8 could also have been derived from the formula (3.6) which 
involves the Fourier transform of the 6-function. In our problem 


(4.10) fi» f exp {iy[» — vi{(1 +7) 
— cos — cos $2}/2(1 + 


An application of the formula 
1 
(4.11) Jo(x) = cos 099 
and a simple transformation yields 
(4.12) vrg(v) = 2N°(1 + fx“ f y, 


With the aid of the two Bessel function formulae 


mJ o(z)Jo(Z) = + 2? — 2Zz cos 6)d0 


and 


f Jo(at) cos bt dt = 


if 


(4.12) can easily be shown to be equivalent to (4.6)—(4.8). 

We shall now show that the saddle points and logarithmic singularities dis- 
cussed above also exist in 2D lattices of particles which vibrate in the plane of 
the lattice. It was in this type of system that the existence of the logarithmic 
singularities was first noticed [5]. 

The simplest example of a 2D system with vibrations in the plane j is one 
with interactions between nearest neighbors only. One would expect the coupling 
between horizontal and vertical components of lattice vibrations to be very 


weak. On that basis let us consider as a model a system with the potential energy 
function 


d= (t61,m — %141,m)? + (41,m — %t,m41)? 
(4.13) 


+ (P2,m — + $82 — Vr41,m)?. 
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The constants a and #; are force constants associated with displacements 
parallel to the line connecting two particles while a: and f: go with displace- 
ments perpendicular to the connecting line. If all forces between particles are 
central, a2=$2:=0. Since the non-central contribution to forces in a crystal is 
usually small, we assume that asa; and 6.«6;. The 3D analogue of this 
model has been discussed by Newell and Rosenstock [13]. 

The characteristic determinant (2.7) for this system is diagonal. The reader 
can easily verify that the frequencies of normal modes are 


= + Bz) — 28; cos — cos 


where ¢; and @z satisfy (2.5). 

Since these frequencies are generated by exactly the same equation, (4.2), 
as the generating equation for normal mode frequencies in a system with trans- 
verse vibrations, the FS, g(v), must have the form as that plotted in Figure 5. 
In a symmetrical lattice (a1, a2) = (81, and there are two logarithmic singu- 
larities as long as a2~a; (the case a; =a, is quite unrealistic physically). In an 
unsymmetrical lattice with (a, a2) and none of the a’s 
or B's equal to zero, there are four singularities; the removal of a symmetry 
doubles the number of singularities. 

When the forces are central (a2 =82=0), g(v) becomes exactly that of a one- 
dimensional system. In order to restore the two-dimensional character to a lat- 
tice of particles with central forces and vibrations in the lattice plane only, one 
must take into account interactions between more distant neighbors. 

The calculation of g(v) becomes rather complicated when interactions be- 
tween more than nearest neighbors are included. The case of interactions be- 
tween nearest and next nearest neighbors has been discussed at considerable 
length by the author [5]. Here we shall merely give a statement of the problem 
and summarize the results. 

We restrict ourselves to symmetrical systems with central forces (with inter- 
actions only along lines which connect particles). A typical set of interactions 
are those along lines through the vectors given in Figure 7. The total potential 
energy of interaction between all particles is (u:,m and 0;,m represent horizontal 
and vertical displacement of the (J, m)-th particle) 


tad — %r41,m)? + — 


+ (t62,m — — + 


In the usual manner the characteristic equation is found to be 
(1—7) (1-01) +7(1—crce) —2f? TS1S2 
(1—1)(1—c2) +7(1—¢1c2) —2/? 


(4.16) 0 
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where 
(4.17) cj = Cos gj, $s; = sing; and f = 


The largest frequency is given by 


(4.18) vi = + 
and 
(4.19) = = 1/[1 + (a/2y)]. 
e 
I,m 
Fic. 7. Typical directions of interaction between nearest and next nearest neighbors in a 
two dimensional lattice. 


The frequencies of normal modes are generated by 
f = (y/yz)? = 3[27(1 — + (1 — 
+ 4 — 1)(4 1) + (1 — 


as the values of $; run through their usual range (2.5). 

An unfortunate feature of (4.20) is that it involves radicals. The integra- 
tion of (3.4) or (3.6) does not give a closed form result in terms of well known 
functions. However, when r=1/3, the quantity inside the radical is a perfect 
square, and the two branches of frequencies become: 


(4.20) 


(4. 21a) fe = 4(201 aa) — 


(4. 21b) 2. 


The branch (4.21b) has a form that has already been discussed (r =1 in (4.5)). 
The “+” branch can be analyzed easily and also leads to a FS that is a com- 
plete elliptic integral of the second kind. Both branches have logarithmic 


‘ 
3 
i 
j 


FREQUENCY SPECTRUM OF VIBRATIONS OF A CRYSTAL LATTICE 61 


singularities (the “+-” branch at f?=3/4 and the “—” branch at f?=1/3) at 

values of f which correspond to saddle points. The curves of constant frequency 

(in one quadrant of (¢:, ¢2)-space) of the “++” branch are drawn schematically 
in Figure 8. 

T 

%, %, 


0 T 0 0 


%, >, 


a b c 


Fic. 8. Curves of constant frequency in the three ranges. The arrows point in the direction 
of increasing frequency. 

Although one cannot obtain simple expressions for g(v) for other values of r, 
it is not difficult to plot the curves of constant frequency for any fixed r. For all 
values of r<1/2 the “—” branch has a saddle point at f?=r and the *+” 
branch has one at f?=1—7 if r<1/5 and at f?=(1+3r)?/16r if 1/5<7r<1/2. 
It can be shown that logarithmic singularities exist at each of these saddle 
points and therefore that there is a singularity for each branch of the spectrum. 

Smollett [7] has plotted out the curves of constant frequency for a lattice of 
alternate positively and negatively charged particles, all of equal mass, to de- 
termine the effect of long range coulomb forces on their character. These curves 
are somewhat more complicated than those with short range forces, but the 
saddle points and singularities still exist. 

The qualitative similarity between the frequency spectra of all the models 
discussed leads one to ask three questions; (a) do saddle points generally lead 
to logarithmic singularities; (b) is there a physical reason for the existence of 
the saddle points; and (c) would the singularities remain if the periodic boundary 
conditions were replaced by some other ones? We shall now deduce an affirma- 
tive answer to all of these questions. 


5. Topological basis for logarithmic singularities of the 2D FS.* We first 
show that saddle points in (¢;, ¢2)-space always lead to logarithmic singulari- 
ties in g(v) (see also Smollett). If the saddle point is not on the ¢; or ¢: axis, let 
us rotate the coordinate system so that this becomes the case. If it is not sym- 
metrical so that 

2 2 
= — 2a* <0, 


= 


* The topological discussion is too oversimplified for many applications to FS calculations. 
The reader is referred to the work of van Hove for a more detailed and rigorous analysis. The 
author is indebted to Professor van Hove for an informative discussion on which part of this section 
is based. 
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let the scale of one of the variables be changed so that this becomes true at the 
saddle point, (¢%, 0). Since 
of? af? 
rye 
at (6, 0), the Taylor expansion of f? about (¢®, 0) is 


(5.2) 


where f, is the value of v/v, on the curve of constant frequency which goes 
through (¢, 0). The value of the ¢; intercept, ¢@, in the neighborhood of the 
saddle point is related to f through 


(5.3) = (fh 
The contribution of the shaded area of Figure 9 to vzg(v) in a range 


%, 


Fic. 9. Curves of constant frequency in the neighborhood of a saddle point. 


$® $¢1<5¢© in which (5.2) is a good approximation to f? is 


2 2 1/2 
N*f (0) 
— 1)? — (f2 — 
du 
N*f 
(5.4) =— log (f. — f) + bounded function of f. 


8x2a? 
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As the frequency v approaches », (i.e., as f—>f,) this term contributes the 
logarithmic singularity to the FS as we set out to prove. In a similar manner it 
is easy to obtain the same result as y, is approached from higher frequencies. 

Van Hove [10] has recently used the critical point theory of the calculus of 
variations in the large (developed by M. Morse [11]) to show that our saddle 
points are a direct consequence of the periodicity of a crystal lattice and the 
Born-Karman boundary conditions. 

In preparation for the application of Van Hove’s arguments we present a 
heuristic discussion of a special case of one of M. Morse’s theorems. Let G(@y, $2) 
be a continuous doubly periodic function of ¢; and ¢2 with continuous first and 
second derivatives, the derivatives vanishing only at isolated points. Then 
Morse’s theorem states that G has at least two saddle points. As was mentioned 
in section 3, a doubly periodic function is equivalent to a function defined on a 
torus; hence this theorem is equivalent to the statement that a function 
G(¢:, $2) with the properties postulated above and defined on a torus has at 
least two saddle points. 

Since G has been chosen to be continuous it must have at least one maximum 
and one minimum point in each of its periods. Let the location of one of these 
maxima be denoted by an xX in each period in Figure 10 and let the location 
of one set of equivalent minima be represented by heavy dots (points D, E, etc.) 


a 
%, 
xe x 
x x 


Fic. 10. Curves connecting maxima and minima of functions defined on a doubly periodic space. 


If the maxima A and B are connected by a curve such as 1 in Figure 10, 
there is at least one point on the curve where G has a smaller value than at 
neighboring points on the curve. G also achieves a smallest value along any other 
curves, for example on 2 and 3, which connect A and B. The locus of these 
points forms a continuous curve which passes through E and D. G must achieve 
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a largest value somewhere along this locus of smallest values. Let this be ac- 
complished at the point represented by the small triangle in Figure 10. This 
point must be a saddle point because as one passes along 4 from E to D it isa 
relative maximum, and as one passes from A to B it isa relative minimum. 
The same argument could be applied by the examination of paths which 
connect A and C. Hence the function G has at least two saddle points. It is 
easy to trace out the same proof when one period of (¢:, ¢2)-space is wrapped 
into a torus. Those paths which start at x and loop the torus, as in Figure 11a, 
correspond to paths from A to C in Figure 10, while those that go around the 


a b 


Fic. 11. Closed curves through maxima and minima of a function defined on a torus. 


hole (as in Figure 11b) correspond to paths from A to B. No matter how the 
path which goes around the hole is deformed or displaced it cannot be trans- 
formed continuously into the one which goes through the hole. The number of 
saddle points of a function defined at all points on the surface of a 3D body is at 
least equal to the maximum number of closed paths which cannot be trans- 
formed into each other or into a point by a continuous deformation of the paths. 
This number is known as the Betti number of the surface. 

The fact that there might be two saddle points in our reciprocal (¢1, ¢2)- 
space does not guarantee the existence of two logarithmic singularities in g(v). 
The value of »(¢;, ¢2) at the two saddle points may correspond to the same fre- 
quency. This is the case when a square lattice is elastically symmetrical. We see 
in Figure 3 that the two saddle points, one at (¢1, ¢2) = (0, 7) and the other at 
(x, 0) (as well as those at (0, —7) and (—z, 0)) both correspond to the fre- 
quency v= (4)"/*v, when r=1. 

On the basis of the general theorem concerning saddle points, Van Hove 
has pointed out that the logarithmic singularities in the FS of a 2D lattice are 
a direct consequence of the postulate of the harmonicity of the interatomic 
forces and of the doubly periodic structure of the lattice. The number of inter- 
acting neighbors and the detailed nature of the force constants is not critical 
as long as the series F,; in (3.6) converge and as long as the lattice is not so 
weakly coupled that the relaxation of an initial displacement of a given atom 
propagates a wave only along a single row and column (as is the case of a 


square lattice with central force interactions which act only between nearest 
neighbors). 
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The last of our questions about the singularities of the 2D g(v) is whether 
their existence is sensitive to the choice of periodic boundary conditions. To 
show that the detailed choice of boundary conditions is unimportant we apply 
the following argument of Ledermann [15]. It is based on Ledermann’s theorem, 
If, in a Hermitian matrix the elements of r rows and their corresponding columns 
are modified in any way whatsoever, provided only that the matrix remains Her- 
mitian, iten the number of characteristic values which lie in any interval cannot 
increase or decrease by more than 2r. 

The particles in an NXWN square lattice with interactions between nearest 
and next nearest neighbors can be divided into (N—2)? internal particles and 
(4N—4) boundary particles. We might identify the internal ones by the num- 
bers 1, 2, - - -, (N—2)? and the boundary ones by (V—2)?+1,---,N*% The 
equations of motion of the internal particles are all of the same form. 

The matrix of the secular determinant can then be expressed in the form 


(= 2(N 2)? rows. 
My/ 2(4N — 4) rows. 


(5.5) 


where M;; represents the matrix elements of the interactions between internal 
particles, M;, and M,; those between internal and boundary particles, and My» 
those between boundary particles. The corresponding matrix for Born-Karman 
boundary conditions would be of the form 


Mi Be 
(5.6) ( ) 

Bu Bus 
and hence differs from (5.5) in the elements of 2(4N —4) rows. Then the maxi- 
mum number of normal modes which move into or out of a given frequency 
interval is at most 16(N—1) out of the order of 2N*. Hence in the limit as 
N-—»o a negligible fraction make the move, and in this limit g(v)/2N? is inde- 
pendent of the boundary conditions. 

Even in systems with longer (but fixed) range interactions, the boundary 

thickness is independent of N so that the number of rows changed in (5.5) in 
going to periodic boundary conditions is still of O(N). 


6. The nature of the singularities in three-dimensional lattices. Van 
Hove has also discussed the singularities of the 3D FS by topological con- 
siderations. The reciprocal lattice of a 3D one component cubic crystal is 
also 3D. For convenience we shall restrict ourselves to simple cubic lattices. We 
let ¢;=271/N3, with / ranging from 1 to N3, be the third dimension in the RL. 
Then each normal mode frequency is a triply periodic function of the ¢’s, 
v(d1, 2, 3). The 3D RL can be visualized as a set of parallel layers of 2D RL’s 
with $3; constant in each layer. The complete FS of a crystal is the sum of the 
spectra of each layer in the RL. 

The Morse theorem is applicable to each of these layers. The FS associated 
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with each has at least one or two saddle points, the minimum number depending 
on symmetry. If the saddle point remains at the same frequency on all layers, 
the total FS, being a sum of functions each with a logarithmic singularity at the 
same point, has a logarithmic singularity. This, however, is an exceptional case 
which occurs only when the force constants satisfy special conditions. 
Generally the position of the saddle point varies from layer to layer. A 
schematic plot of ».(¢3), the singular frequency, is made in Figure 12. This fre- 
quency, being a continuous periodic function of $3 has at least one maximum 
(say v at ¢%) and one minimum (say v°” at $$”). In the case of a single maxi- 


3 


Fic. 12. Variation of singular frequency with layer. 


mum and minimum all the logarithmic singularities are located between »°” 
and v). Hence we can expect g(v) to have a sharp rise in the neighborhood of 
v and a sharp drop in the neighborhood of v. In the limit as N;— these 
singularities become dense so that if one divides the FS by N; as the limit is 
taken, the final g(v) should vary only slowly in the frequency interval (»%, »®) 
(Figs. 13 and 14). 

We shall now investigate this behavior analytically for frequencies near 
v and v, The singular frequencies on layers of the RL with ¢; close to $$” 
are given by 


(6.1) = ve +0 os) 


where 
1 3? c 
> 0. 
2 
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The FS of a ¢; layer close to ¢{? and at frequencies near v{) is proportional to 
(see (5.4)) 


(6.2) —tlog |»<(s) — 


g (Vv) 


/\A\ 
//A\\ 
\\ 
\\ 
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Fic. 13. Frequency spectra associated with various layers in a thin crystal. 


We choose ¢2 to be an arbitrary number such that (6.1) is valid for all $3 in 
the range ¢° <¢;<¢%. Then the total contribution of the region of the RL 
in the neighborhood of the saddle points near ¢{” to the frequency spectrum is 
proportional to 


Ns 
(6.3) —-2— log | v<(¢s) — v| dos 


2r J 4, 


for those frequencies near vy“. The factor 2 is introduced so that the integration 
can stop at $2 instead of proceeding to a distance $$? —¢9 beyond $$”. 

We shall now show that the FS rises sharply as vv!” (see Fig. 14) from the 
left. When »<v™, (6.3) becomes 


N 
log — ») + — ) 


Ns 2 2 
> f log (a? + u*)du 


N 
~ {c log (a? + c?) — 2c + 2a tan c/a} 


e 
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where 


a= |r —v|}, u = b (vs — andc = (gs — da)b > 0. 


In the limit as (since tan —! 


(6.4) { — — + 2c(1 — log c)}. 


Since the slope of J, is proportional to 1/(v‘ —v)"/?, the FS, g(v), has a vertical 
tangent at v=v" (see Fig. 14). 

It can now be shown that g(v) is continuous at v® but that g’(v) is discon- 
tinuous. When v>v we define ¢¥ so that 


g(v) 


ye ye 


Fic. 14. Typical peak in a 3-D frequency spectrum. 


Then (6.3) becomes 


N 
{ log [b'(¢s. — #3) — a 


(0) 
+ log [a — — 63) 
= = log (u® — a?)du + log (a? — 
6.5) 
Ns; 
= “—— {c log (c? — a*) + a log [(c + a)/(c — a)| — 2c} 
— log) — — )/e + — }: 


as vw”, (6.4) and (6.5) approach each other. Hence g(v) is continuous at 


However, when from the left while a 
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finite number, as vv" from the right. Hence g’(v) is discontinuous. 

One can proceed with a similar analysis to find that g(v) has a rapid drop 
with a vertical tangent at v®, the location of the highest frequency saddle 
point in the RL. This behavior is exhibited in Figure 14. 

The 3D RL can also be built up in layers which are normal to either the ¢, 
or ¢, axes. When the surfaces of constant frequency are not symmetrical func- 
tions of $i, ¢: and ¢; the argument used above can be repeated in each of the 
new alternative layer constructions to locate other singularities. 

Let us consider a case which has symmetrical surfaces of constant fre- 
quency and which has simple peaks of the type given in Figure 14. Then let a 
symmetry be removed by making a small change in the force constants so that 
contours in planes normal to two of the axes are identical but different from 
those on planes perpendicular to the third axis. The total FS is the superposition 
of pairs of peaks with their edges slightly displaced from each other as in Figures 


Fic. 15. Displacement of peaks by slight perturbation in symmetry. 


Fic. 16. Typical peaks in the FS for anisotropic 3-D crystals. 


15a or 15b. The total curve, being the sum of these pairs, generally has peaks of 
the form of those in Figure 16a. Under special conditions only one of the edges 
might be displaced by the perturbation of the force constants, thus yielding 
curves such as those in Figure 16b. When larger perturbations are made one of 
the members of the pair of peaks may be displaced sufficiently to split the 
original peak into two as in Figure 16c. 

Finally when small displacements in force constants are made which remove 


all symmetries with respect to ¢1, ¢2, and $3, three layered peaks such as that in 
Figure 16d may appear. 


) 
| 
\ \ 
/ \ 
a b 
L 


70 SYMPOSIUM ON SPECIAL TOPICS IN APPLIED MATHEMATICS 


We shall now examine some models in which these effects can be studied 
analytically. 


7. Some examples of 3D frequency spectra. As in the case of 2D lattices we 
shall first discuss a model in which one degree of freedom is associated with 
each lattice point ina NXNXVN 3D lattice with interactions between nearest 
neighbors. If the potential energy of interaction is chosen to be 


= > (Zt,m.n — Zt,m.nti)? + — 


7.1 
( ) + das (Z2,m,n 


and periodic boundary conditions are employed, the frequencies of normal 
modes are given by 


(7.2) = 2(a, + a2 + as) — 2a; cos — 2a2 cos — 2a3 cos $3 
where 

= 2nj/N, go: = 2rk/N, $3 = 2ns/N, j, ks =0,1,2,---,N—1. 
As usual we let f=v/vz, so that 
(7. 2a) (1 — 2f?)(a1 + a2 + as) = a COS gi + a2 COS + a3 COS 
with 
(7.2b) 4n = 4(a; + a, + as). 


The saddle points on planes normal to the ¢3 axis occur at the frequencies 


2 a3 + 2ae — COS o3 
7.3 = 
+ a+ a3) 
and 
(7.4) re ) a3 + 2a; — a3 COS o3 


2(a + a2 + as) 
The largest value of f? at a saddle point associated with (7.3) is 


(1), 2 


(7.5a) (fe ) = (a2 + a3)/(a1 + a2 + as) 
while the smallest is 
(7.5b) = + a + as). 
The edges of the peak associated with (7.4) are at 

(2),2 


(7. 6a) (fe ) = (a1 + a3)/(a1 + a + a) 
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on the right and 
(2),2 
(7.6b) (fe) = + + as) 


on the left. A similar analysis of saddle points on planes normal to the ¢,; and @¢2 
axes yield the following four peaks with edges at f? equal to 


(7.7) + a2 +3) and (a + + a2 + as) 
(7.8) as/(a; + 03) and (as + + a2 + as) 
(7.9) + a2 + as) and (a2 + a1)/(a; + a2 + a) 
(7.10) + a2 +3) and (a3 + a)/(a; + a2 + as). 


The qualitative appearance of the FS depends on the particular values assigned 
to the force constants a), a2, and a3. We shall consider several interesting special 


cases. 
When a, =a2:=a; we get a single peak with edges at f?=1/3 and f?=2/3 
as in Figure 17a (and as was first obtained by Bowers and Rosenstock [6]). 


lA LA LA 
a b c d 


Fic. 17. Variation of frequency spectrum of 3-D lattice with one degree of freedom per 
lattice point with the reduction of coupling between layers. 


If we let a;=a.=a and a,=a(1+6) with | 8| <1, the peak edges appear at f? 
equal to (3, 3), ((1+4), 3(2+4)) (4, 3(2+46)), so that the FS has the form of 
Figure 17b. 

If 5=e€—1 with € a small positive number the edges are at (4, 3); (fe, 4(1+¢€)); 
(3, 3(1+¢€)). In the limit as e—0 the system degenerates into a set of inde- 
pendent layers while the peak of the FS narrows (Fig. 17c) until in the limit 
it has the 2D logarithmic singularity (17d). 

The peak in 17a splits into two when 6>1, for then 4(1+6) > 3 and the FS 
has the appearance of Figure 18a. 
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Let us now consider a model which is more representative of a real crystal; 
7 namely, a simple cubic lattice of springs and masses with three degrees of 
LA LA JU 
v/v W/V. Wu 
a b c d 


Fic. 18. Variations of frequency spectrum of simple cubic lattice (with interaction between 
nearest neighbors only) as coupling between layers weakens. 


freedom each and with interactions between nearest neighbors only. Then 
generalizing (4.13) we choose the potential energy function to be 


+ as p> — 

+ 3Bs = — 

+ in — + (wWimin — 


The u, v, and w’s represent components of displacements in the x, y, and z direc- 
tions while a, 6;, and 7; are force constants associated with displacements in the 
direction of the lines which connect particles. If forces are central all other a, 8, 
and y’s vanish. The introduction of periodic boundary conditions yields the 
following normal mode frequencies: 


(7.11) = 2(a, + a2 + a3) — 2a cos di — 2a2 cos — 2a3 cos $3 
(7.12) 42*v?M = 2(6; + Be + Bs) — 283 cos di — 28; cos — 282 cos 
(7.13) = 2(y1 + v2 + ¥3) — 272 cos di — 2y3 Cos — 2y1 Cos $3; 


the ¢’s have the usual range. These equations are of the same form as (7.2). 
Hence the discussion given above is immediately applicable. 

First we consider the symmetrical simple cubic lattice with a:=$:=73 
=a(1+5) and a.=a;=8; =¥2=a. Since the noncentral character of the 
forces acting on atoms in a solid is small, we assume 6>>0. When 6>1 two peaks 
exist as is shown in Figure 18a. Detailed calculations have been made on this 
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model by Newell and Rosenstock [13]. This type of FS is typical for single 
cubic lattices with interactions between more distant neighbors. Newell [12] 
has shown this through several exact numerical calculations. 

It is interesting to see how the FS in Figure 18a degenerates into that of a 
two dimensional lattice as the interactions between layers are made weaker. 
We let a1=ys=a(1+6), B2=ea(1+6), Bi =a2=B;s=y2=a and y:=a;=ea and 
let e—0. Then the peaks become progressively narrower until in the limit they 
become logarithmic. 

If one wishes to find the FS of a real crystal, several courses are available. 
He should first locate the various singularities and the low frequency behavior 
by analytical methods. Then either by numerical integrations based on the 3D 
analogues of (3.4) or (3.6) or by actual calculation and sorting of normal mode 
frequencies on a high speed computing machine, a FS curve can be sketched 
which is consistent with the location of singularities. Another scheme is to 
choose an empirical function which is a linear combination of one with properly 
located singularities and a polynomial (or the product of these two types of 
functions) and find the values of the constants which lead to independently 
calculated moments of the FS. These moments are easily obtained from the 
characteristic equations which define the normal mode frequencies [16]. 
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